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This review deals with the structure of hadrons, strongly interacting many-body systems consisting of quarks and 
gluons. These systems have a size of about 1 fm, which shows up in scattering experiments at low momentum 
transfers Q in the GeV region. At this scale the running coupling constant of Quantum Chromodynamics (QCD), 
the established theory of the strong interactions, becomes divergent. It is therefore highly intriguing to explore 
this theory in the realm of its strong interaction regime. However, the quarks and gluons can not be resolved at 
the GeV scale but have to be studied through their manifestations in the bound many-body systems, for instance 
pions, nucleons and their resonances. The review starts with a short overview of QCD at low momentum transfer 
and a summary of the theoretical apparatus describing the interaction of hadrons with electrons and photons. 
In the following sections we present the experimental results for the most significant observables studied with 
the electromagnetic probe: form factors, polarizabilities, excitation spectra, and sum rules. These experimental 
findings are compared and interpreted with various theoretical approaches to QCD, such as phenomenological 
models with quarks and pions, dispersion relations as a means to connect observables from different experiments, 
and, directly based on the QCD lagrangian, chiral perturbation theory and lattice gauge theory. 
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"This review is dedicated to the memory of the late Peter Brix (MPI Heidel- 
berg), the Nestor of low Q 2 physics. 



[. INTRODUCTION 

Hadrons are composite systems with many internal degrees 
of freedom. The strongly interacting constituents of these 
systems, the quarks and gluons are described by Quantum 
Chromodynamics (QCD). This theory is asymptotically 
free, that is, it can be treated in a perturbative way for very 
large values of the four-m o mentum transf er squared, Q 2 
(iGross and Wilczekl Il973alb1: iPolitzerl Il973l) . However, the 
binding forces become increasingly strong if the momentum 
transfer decreases towards the region of about 1 GeV, which 
is the natural habitat of nucleons and pions. In particular, 
the "running" coupling constant of the strong interaction, 
a s (Q 2 ), is expected to diverge if Q 2 decreases to values 
near Aq CD w (250 MeV) 2 , which defines the "Landau pole" 
of QCD. This behavior is totally different from Quantum 
Electrodynamics (QED), for which the coupling constant 
Q: e m(C} 2 ) diverges for huge momentum transfers at the Planck 
scale, corresponding to Q 2 w 10 38 GeV 2 or 10~ 35 m, much 
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below any distance ever to be resolved by experiment. On the 
contrary, the Landau pole of QCD corresponds to a resolution 
of the nucleon's size, somewhat below 1 fm or 1CP 15 m. This 
is the realm of non-perturbative QCD, in which quarks and 
gluons appear as clusters confined in the form of color-neutral 
hadrons. As of today it is an open question whether this 
confinement can be derived directly from QCD or whether it 
is a peculiarity of a strongly interacting many-body system or 
based on some deeper grounds. Therefore, the study of QCD 
in the non-perturbative domain serves less as a check of QCD 
per se, but is concerned with the highly correlated many-body 
system "hadron" and its effective degrees of freedom. 

Quantum Chromodynamics is a non-abelia n gauge theory 
developed on the basis of quarks and gluons (iFritzsch et al. , 
11973b iGross and WilczekL Il973bt IWeinbergl Il973l) . The non- 
abelian nature of this theory gives rise to a direct interac- 
tion among the gluons, and the forces among the quarks 
are mediated by the exchange of gluons whose chromody- 
namic vector potential couples to the vector current of the 
quarks. If massless particles interact via their vector cur- 
rent, the helicity (handedness or chirality) of the particles 
is conserved. The nucleon is essentially made of the light 
u and d quarks plus a small admixture of s quarks, with 
masses m u = 1.5 to 3.0 MeV, nid = 3 to 7 MeV, and 
m s = (95 ± 25) MeV dYao et all 120061) . In the zero mass 
limit, these light quarks can be classified according to their 
chirality by the group SU(3).r<£> SU(3)l. Several empirical 
facts give rise to the assumption that this symmetry is spon- 
taneously broken down to its vectorial subgroup, and in addi- 
tion the finite quark masses cause an explicit symmetry break- 
ing. The spontaneously broken symmetry is a most remark- 
able feature of QCD, because it can not be derived from the 
Lagrangian. This is quite different from the explicit symme- 
try breaking, which is put in by design through the finite quark 
masses in QCD and appears in a similar way in the Higgs sec- 
tor. As a result one obtains the conserved vector currents J" 
and the only partially conserved axial vector currents J^, 

J l. = 97m Y q ' J ^ = 97^75 y9 , (1) 

where q are Dirac spinors of point-like (light) quarks and 
7 M ,75 the appropriate Dirac matrices. The quantities A°, 
a = 1 ... 8 denote the Gell-Mann matrices of SU(3) describing 
the flavor structure of the 3 light quarks, and A is the unit ma- 
trix. The photon couples to the quarks by the electromagnetic 
vector current J™ ~ + jjj , corresponding to isovec- 
tor and isoscalar interactions respectively. The weak neutral 
current mediated by the Z° boson couples to the 3 rd , 8 th , and 
th components of both vector and axial currents. While the 
electromagnetic current is always conserved, 9 M JP™ = 0, the 
axial current is exactly conserved only for massless quarks. 
In this limit there exist conserved charges Q a and axial 
charges Q5?, which are connected by commutation relations. 
The corresponding "current algebra" predated QCD and 
was the basis of various low-energy theorems (LETs), which 
govern the low-energy behavior of (nearly) massless particles. 



The puzzle we encounter in the physics of hadrons is the 
following: The massless quarks appearing in the QCD La- 
grangian must conserve the axial currents. The nucleons 
should eventually emerge from the same Lagrangian as mas- 
sive many-body systems of quarks and gluons. However, the 
conservation of the axial current in the Wigner-Weyl mode 
would require a vanishing axial coupling constant for these 
massive nucleons, which is ruled out by the observed (3 decay. 
A solution of this puzzle was given by Goldstone's theorem. 
At the same time as the "three-quark system" nucleon be- 
comes massive by means of the QCD interaction, the vacuum 
develops a nontrivial structure due to finite expectation val- 
ues of quark-antiquark pairs (condensates (qq)), and so-called 
Goldstone bosons are created, qq pairs with the quantum num- 
bers of pseudoscalar mesons. These Goldstone bosons are 
massless, and together with the massive nucleons they warrant 
the conservation of the axial current. Because the quarks are 
not really massless, the chiral symmetry is slightly broken in 
nature. As a consequence also the physical Goldstone bosons 
acquire a finite mass, in particular the pion mass follows 
to lowest order from the Gell-Mann-Oakes-Renner relation, 

m lfl = -{m u + m d )(qq) H , (2) 

with the condensate (qq) » -(225 MeV) 3 and f„ ps 93 MeV 
the pion decay constant. Since the pions are now massive, the 
corresponding axial currents are no longer conserved and the 
4-divergence of the axial current becomes finite, 

d^J^ « -Umltt , (3) 

where (j>^ describes the local pion field. In other words the 
charged pion decay 7r + — * /i + + and 7r~ — * + 
proceeds via coupling to the axial current of Eq. (0. While the 
charged pions decay weakly with a life-time of 2.6 x 10~ 8 s, 
the neutral pion decays much faster, in 8.4 x 10 -17 s, by means 
of the electromagnetic interaction, ir° — > 7 + 7. This provides 
an additional source for the neutral pion or the axial current 
with index 3, 

d»JI=^E.B, (4) 

where E and B are the electromagnetic fields. We note that 
the scalar product of the two electromagnetic fields is a pseu- 
doscalar. This decay can be mediated by a triangle of inter- 
mediate quark lines, and therefore it is often called the tri- 
angle anomaly. It is "anomalous" because such processes do 
not appear in classical theories but only in quantum field the- 
ories through the renormalization procedure (Wess-Zumino- 
Witten term). The analogous anomaly in QCD is obtained 
from Eq. (01 by replacing the electromagnetic by the corre- 
sponding color fields, E c and B c , a em by the strong coupling 
constant a s , and with an additional factor 3 for u, d, and s 
quarks, 

9"J 5 ° = 3 — E c -B c . (5) 

As a consequence also J®^ is not conserved, not even for 
massless quarks ("J7a(1) anomaly"). 
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Unfortunately, no ab-initio calculation can yet describe the 
intriguing but complicated world of the confinement region. 
In principle, lattice gauge theory should have the potential to 
describe QCD directly from the underlying Lagrangian. This 
theory discretizes QCD on a four-dimensional space-time lat- 
tice and approaches the physical world in the continuum limit 
of vanishing lattice constants dKogut and Susskincfl, 1 1 975b 
IWilsonl 1 1 974b . However, these calculations can as yet only 
be performed with u and d quark masses much larger than the 
"current quark masses" mentioned above, and therefore also 
the pion mass turns out much too large. As a consequence the 
Goldstone mechanism, the abundant production of sea quarks 
is much suppressed. Lattice gauge theory has progressed 
considerably over the past decade and further progress is 
foreseen by both improved algorithms and increased com- 
puting power. For recent developments see the fo ll owing 
references: dAlexandroui 2007: lAlexandrou et all 2006 ; 
Boinepalli et al. 2006; Edwards et all 120061: iGockeler et a 



20061) . Semi-quantitative agreement has been reached for 
ratios of masses and magnetic moments for the hadrons, 
there also exist predictions for nucleon resonances and 
electromagnetic form factors in qualitative agreement with 
the data. However, some doubt may be in order whether 
such procedure can ever fully describe the pionic degrees of 
freedom in hadronic physics, particularly in the context of 
pion production and similar reactions. 

A further ab-initio calculation is chiral p erturbation theory 
(ChPT), which has been established by IWeinbergl d 19791) 
in the framework of effective Lagran gians and put into a 
systematic perturbation theory by iGasser and Leutwvlerl 
(11984 119851) . This theory is based on the chiral symmetry 
of QCD, which is however expressed by effective degrees 
of freedom, notably the Goldstone bosons. Because of the 
Goldstone mechanism, the threshold interaction of pions and 
other Goldstone bosons is weak not only among themselves 
but also with the nucleons. Furthermore the pion mass 
is small and related to the small quark masses m u and 
nid according to Eq. @. Based on these grounds, ChPT 
has been set up as a perturbation theory in the parame- 
ters p :— (pi,P2, Tn u , rrid), where pi are the external 
4-momenta in a particular (Feynman) diagram. Chiral 
perturbation t heory has been applied t o man y photoinduced 
reactions by iBernard et all dl991al 1 19951 1 19931) in the 
1990's. As a result several puzzles have been solved and 
considerable insight has been gained. However, ChPT can 
not be renormalized as QED by adjusting a few parameters 
to the observables. Instead, the appearing infinities must 
be removed order by order in the perturbation series. This 
renormalization procedure gives rise to a growing number of 
low-energy constants (LECs) describing the strength of all 
possible effective Lagrangians consistent with the symmetries 
of QCD, at any given order of the perturbation series. These 
LECs, however, can not (yet) be derived from QCD but must 
be fitted to the data, which leads to a considerable loss of 
predictive power with increasing order of perturbation. A 
further problem arises in the nucleonic sector because of 
the large nucleon mass M, which is of course not a small 



expansion parameter. The latter problem was solved by heavy 
baryon ChPT, a kind of Foldy-Wouthuysen expansion in 
1/M. This solution was however achieved at the expense of 
approximating the relativistic description by a non-relativistic 
one. Over the past decade new schemes have been developed, 
which provide a consistent e xpansion within a ma n ifestly 
Lorentz invariant formalism dBecher and Leutwvlerl 1 19991: 
Fuchs et aZ.l 120031; iKubis and Meifineiil200ltlSchindler et all 
2004). Fo r recent reviews of ChPT see the work of ISchereri 
(120031) and lBernardl d2007l) . 

If quarks and gluons are resolved at high momentum 
transfer, they are asymptotically free and their momentum 
distribution can be described by evolution functions as 
derived from perturbative expansions ("higher twists"). 
This domain has been studied in great detail ever since the 
discovery of parton scaling at the end of the 1960's. Such 
investigations have given confidence in the validity of the 
QCD Lagrangian. Systems of heavy quarks (charm, bottom, 
top) can be well described by effective field theories based 
on QCD. However, these approaches are less effective for 
systems of light quarks (up, down, strange), for which the sea 
quarks and notably pionic degrees of freedom become very 
important. In order to incorporate the consequences of chiral 
symmetry, a plethora of hybrid models with quarks and pions 
has been developed. Quark models have been quite successful 
in predicting the resonance spectrum of the nucleon as well as 
the electromagnetic decay and excitation of these resonances. 
However, they have problems to describe the spectrum and 
the size of the nucleon at the same time. We do not dwell on 
these models in the review but occasionally refer to them at 
later stages. 

On the following pages we concentrate on hadronic struc- 
ture investigations with the electromagnetic probe, that is 
electron and photon scattering as well as electro- and pho- 
toproductio n of mesons. A broader ac count can be found in 
the book of Thomas and Weise (200l]). In section^ we give 
a brief introduction to the formalism relevant for these stud- 
ies. The following section [III] summarizes the information on 
the form factors of nucleons and pions. Another bulk prop- 
erty of the hadrons is their polarizability, which can be de- 
termined by Compton scattering as discussed in section [TV] 
These global properties are related to the excitation spectrum 
of the particles, meson production at threshold, resonances 
and continuum backgrounds as detailed in section [V] Finally, 
we examine the origin and relevance of several sum rules in 
section |VI] In all of these fields, there has been a rapid evo- 
lution over the past years triggered by high-precision experi- 
ments. These have been made possible by a new generation of 
electron accelerators with high current, high duty-factor, and 
highly polarized beams in combination with improved target 
and detection techniques, notably for polarized particles. In 
many of the presented phenomena we recover the role of the 
pion as an effective degree of freedom at low energy and mo- 
mentum transfer to the nucleon. It is therefore the leitmotiv of 
this review to look at the hadrons as interesting and compli- 
cated many-body systems whose direct description by QCD 
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proper is a major challenge for particle physics over the years 
to come. 



II. THE ELECTROMAGNETIC INTERACTION WITH 
HADRONS 

A. Kinematics 

Let us consider the kinematics of the reaction 

e(k 1 )+N{p 1 )^e(k 2 ) + N{p2) , (6) 

with ki — (uji,ki) and pi = (Ei,pi) denoting the four- 
momenta of an electron e with mass m and a nucleon N with 
mass M. The 4-momenta are constrained by the on-shell con- 
ditions v\ — v\ — ^l 2 i k\ = k 2 — to2 ' an d by the conser- 
vation of total energy and momentum, k\ + pi = k 2 + p 2 - In 
order to make Lorentz-invariance manifest, it is useful to ex- 
press the amplitudes in terms of the 3 Mandelstam variables 

s = {k 1 +p 1 f , <=(fc 2 -fci) 2 , u={ P2 -k 1 ) 2 . (7) 

Due to the mentioned constraints, these variables fulfill the 
relation s + t + u = 2(m 2 + M 2 ), and therefore we may 
choose s and t as the two independent Lorentz scalars. For 
reasons of symmetry, the center-of-mass (cm) system is used 
in the following. The 3-momenta of the particles cancel in 
this system, and therefore s = (uj cm + E cm ) 2 — W 2 , i.e., the 
Mandelstam variable s is the square of the total cm energy 
W. Furthermore, the initial and final energies of each particle 
are equal, and hence t = — (k 2 — ki) 2 m is related to the 3- 
momentum transfer in the cm system. From these definitions 
it follows that physical processes occur at s > (m + M) 2 and 
t < 0. Because of the smallness of the fine-structure constant 
a em w 1/137, it is usually sufficient to treat electron scatter- 
ing in the approximation that a single photon with momentum 
q = k\ — k 2 = (ui, q) is transferred to the hadronic system. 
We call this particle a space-like virtual photon 7*, because 
t = q 2 < 0, i.e., the space-like component of the 4-vector q 
prevails. Since t is negative in the physical region of electron 
scattering, it is common use to describe electron scattering by 
the positive number Q 2 = —q 2 . This contrasts the situation 
in pair annihilation, e + e~ — > 7*, which produces a time-like 
virtual photon with q 2 = m 2 , > 0. The above considerations 
can be applied to real Compton scattering (RCS), 

-f(k 1 ) + N(p 1 )~, 1 (k 2 )+N(p 2 ) , (8) 

by replacing m\ 2 = k\ 2 — > and to virtual Compton scat- 
tering (VCS), 

7 *(fc 1 )+iV(p 1 ) -> 7 (fc 2 )+7V(p 2 ) , (9) 

by replacing m 2 = k\ — > q 2 < and m 2 = k 2 — > 0. 

Let us now turn to the spin degrees of freedom. The 
virtual photon with momentum q carries a polarization 
described by the vector potential A, which has both a 
transverse component, At -L q, as in the case of a real 



photon, and a longitudinal component q ■ A, which is related 
to the time-like component Aq by current conservation, 
q ■ A = ujAq — q ■ A = 0. Since the electron is assumed to be 
highly relativistic, its spin degrees of freedom are described 
by the helicity h = s ■ k = ±5, the projection of the spin s 

on the momentum unit vector k. In the following we denote 
the polarization of the incident electron by P e = 2h = ±1, 
for example, P e = 1 describes a beam of fully polarized 
right-handed electrons. The polarization vector P of a target 
or recoil nucleon is represented in a coordinate system with 
the z-axis pointing in the direction of the virtual photon, 
e z = q, the y-axis perpendicular to the electron scatter- 
ing plane, e y ~ k\ x k 2 , and the x-axis "sideways", i.e., 
in the scattering plane and on the side of the outgoing electron. 

The scattered electron probes the charge and magnetization 
distributions of the hadronic system via the interaction of the 
electromagnetic currents, which leads to a transition matrix el- 
ement M. = J2 t j,j^( e )J >1 W- If we neglect higher order QED 
corrections, the electron is a Dirac point particle with its cur- 
rent given by j M = — e«27 Ai Wi, where 7 M are Dirac matrices 
and Ui Dirac spinors characterized by the quantum numbers 
i = {ki,hi}. In the one-photon exchange approximation, the 
cross section is then obtained by the square of the transition 
matrix element multiplied by phase space factors, 

da~^|M| 2 ~^^(e)^(ft), (10) 

spins spins 

where r]^ u = j^j* can be calculated straightforwardly. By 
varying the incident electron energy and the scattering angle 
as well as the polarizations of the respective particles, it is 
then possible to enhance or suppress specific components of 
the hadronic tensor W^ v = J^J^, and thus to study differ- 
ent aspects of the hadronic structure in a model-independent 
way. For further details and a general introduction to the struc- 
ture of hadrons a nd nuclei, we refe r to the book of lBoffi et al.\ 
(119961) (see also dBoffi et alllm^ ). 



B. Elastic electron scattering 

The hadronic current for elastic electron scattering off the 
nucleon is given by the most general form for the vector cur- 
rent with the same spin-i particle before and after the colli- 
sion, 

J» = u P2 (7^i(Q 2 ) + *^f F 2 (Q 2 )^u pl , (11) 

where u Pl and u P2 are the 4-spinors of the nucleon in the 
initial and final states, respectively. The first structure on 
the rhs of Eq. (fTTT > is the Dirac current, which describes the 
finite size of the nucleon by the Dirac form factor Fi(Q 2 ). 
The second term reflects the fact that the internal degrees 
of freedom also produce an anomalous magnetic moment k 
whose spatial distribution is described by the Pauli form fac- 
tor F 2 (Q 2 ). These form factors are normalized to Ff (0) = 1, 
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F 2 P (0) = k p = 1.79 and F™(0) = 0, F 2 "(0) 
for proton and neutron, respectively. 
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From the analogy with non-relativistic physics, it is seduc- 
ing to associate the form factors with the Fourier transforms 
of the charge and magnetization densities. The problem is 
that the charge distribution p(r) has to be calculated by a 
3 -dimensional Fourier transform of the form factor as func- 
tion of q, whereas the form factors are generally functions 



of Q z 



q 2 



However, there exists a special Lorentz 



frame, the Breit or brick-wall frame, in which the energy of 
the (space-like) virtual photon vanishes. This can be realized 
by choosing, for example, pi — — \q and p 2 = +\q leading 
to E 1 = E 2 = (M 2 + y 2 )V2, u = o, and Q 2 = q 2 
Equa tion ( fTTT i takes the following form in this frame (ISachsl 
1962): 



Gb(q2) ' ?; w Gm(q2; 



(12) 



where Ge(Q 2 ) stands for the time-like component of J M and 
hence is identified with the Fourier transform of the electric 
charge distribution, while Gm(Q 2 ) appears with a structure 
typical for a static magnetic moment and hence is interpreted 
as Fourier transform of the magnetization density. The Sachs 
form factors G e and G m are related to the Dirac form factors 
by 



Ge{Q ) = Fi(Q 2 ) 
G M {Q 2 ) = F 1 (Q 2 ) 



tF 2 (Q 2 ) 
F 2 (Q 2 ) , 



(13) 



where r = Q 2 /AM 2 is a measure of relativistic (recoil) ef- 
fects. Although Eq. ( fT3T > is a covariant definition, the Sachs 
form factors can only be Fourier transformed in a special 
frame, namely the Breit frame, with the result 



Because the hadronic current is completely defined by 
Eq. (fTTT) . any observable for elastic electron scattering can be 
uniquely expressed in terms of the two form factors. In par- 
ticular the unpolar ized differential cross section is given by 
(RosenblutlQIIi 



da 

dh 



da 



1+T 



2rtan 2 -G 2 vf ) , (15) 



with (da/dCl)o the cross section for electron scattering off a 
point particle and 8 the scattering angle of the electron in the 
laboratory (lab) system. Equation ( fT~5T > gives us the possibility 
to separate the form factors by variation of tan 2 | while keep- 
ing Q 2 constant. In fact the data should lie on a straight line 
("Rosenbluth plot") with a slope that determines the magnetic 
form factor Gm- However, there are limits to this method, in 
particular if one of the form factors is very much smaller than 
the other. In such cases a double-polarization experiment can 
help to get independent and more precise information. Such 
an experiment requires a polarized electron beam and a polar- 
ized target, or equivalently the measurement of the nucleon's 
polarizatio n in the final state. The measured asymmetry takes 
the form (lArnoldef a/Ill 98 ll) 



A= -P t 



v /2re(l - e) G E G M P X 



tVT^G 2 m P z 



eG% 



T G 2 M 



(16) 

where e = 1/[1 + 2(1 + r) tan 2 |] is the transverse polar- 
ization of the virtual photon. In particular we find that the 
longitudinal-transverse interference term, appearing if the nu- 
cleon is polarized perpendicularly (sideways) to q, is propor- 
tional to GeGm, while the transverse-transverse interference 
term, appearing for polarization in the q direction, is propor- 
tional to G 2 M . The ratio of both measurements then deter- 
mines Ge /Gm with high precision, because most normaliza- 
tion and efficiency factors cancel. 



G E (q 2 ) = p(^ f d 3 f 



(14) 



p(r)d?r- — I p{r)r 2 d?r + 
6 



where the first integral yields the total charge in units of e, 
i.e., 1 for the proton and for the neutron, and the second 
integral defines the square of the electric root-mean-square 
(rms) radius, (r 2 )E- We note that each value of Q 2 requires 
a particular Breit frame, i.e., the information on the charge 
distribution is taken from an infinity of different frames, 
which is then used as input for the Fourier integral in terms 
of GE{q 2 )- Therefore, the density p{f) is not an observable 
that we can "see" in any particular Lorentz frame but only 
a mathematical construct in analogy to a "classical" charge 
distribution. The problem is that an "elementary" particle has 
a small mass such that recoil effects, measured by r, and size 
effects, measured by (r 2 ), become comparable and can not be 
separated in a unique way. The situation is numerically quite 
different for a heavy nucleus, in which case the size effects 
dominate the recoil effects by orders of magnitude. 



C. Parity violating electron scattering 

In the previous section we have tacitly assumed that the 
interaction between electron and hadron is mediated by the 
virtual photon and therefore parity conserving. With this as- 
sumption the polarization of only one particle does not yield 
any observable effect. However, it is also possible to exchange 
a Z° gauge boson, although this is much suppressed in the 
low-energy region because of the large mass Mz» — 91 GeV. 
This boson couples to electrons and nucleons with a mixture 
of vector and axial vector currents typical for the weak inter- 
action. If the Z° is emitted from one of the particles by the 
vector coupling and absorbed by the other one by the axial 
vector coupling, it produces a parity-violating asymmetry that 
can be observed if one of the particles (typically the incident 
electron) is polarized. The coupling of the Z° to the electron 
involves the current 



Jit = ffVabvC 1 - 4 sin 2 Q w ) - 7^75]^! 



(17) 



where Qw is the Weinberg angle and g' a weak coupling con- 
stant. Because sin 2 Oyy ~ 0.23, the vector current in Eq. ( TTTb 
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is largely suppressed compared to the axial vector part con- 
taining the 75 factor. The corresponding weak hadronic cur- 
rent can be parameterized as follows: 

J, = u P Al»F l {Q 2 ) +l ^fF 2 {Q 2 ) 

+l^ 5 G A (Q 2 )]u Pl , (18) 

where the tilde signifies the coupling to the Z°. The weak 
Sachs form factors are defined as in Eq. dT"3l ). and the cross 
sections and asymmetries are calculated as in the previous 
section. However, the contribution of the weak current to 
the differential cross section is well below the experimen- 
tal error bars, and information can only be obtained from 
the interference between the electromagnetic and the weak 
transition amplitudes. The parity-violating asymmetry A — 
(da + — da~) / '(da + + da~), where d<r + and da~ are the dif- 
ferential cross sections for incident electrons with positive and 
negative helicities, takes the form 

A = __Gt_ 

( sGeGe + tGmGm - e'(l - 4sin 2 G W )G M G A \ 
X V eG 2 E +rG 2 M J 

= A e (Ge)+A m (G m ) + A a (Ga), (19) 

with e' = y/(l - e 2 )r(l + r). 

D. Pseudoscalar meson electroproduction 

The reaction 

7*(«)+JV(pi) -^n(k)+N( P2 ) (20) 

is described by the transition matrix element e M J M , with e M 
the polarization of the (virtual) photon and J M the transition 
current leading from the nucleon's ground state to a meson- 
nucleon continuum. This current can be expressed by 6 differ- 
ent Lorentz structures constructed from the independent mo- 
menta and appropriate Dirac matrices. Since the photon cou- 
ples to an electromagnetic vector current and the pion is of 
pseudoscalar nature, the transition current appears as an axial 
4-vector in the nucleon sector. The space-like (J) and time- 
like (p) components of the transition operator take the follow- 
ing form in the hadronic cm frame: 

J = aFi + i(q x a) {a ■ k)F 2 + k{a ■ q)F 3 

+ ~k(<7 ■ k)F 4 + q(a ■ q)F 5 + q(<7 ■ k)F e , (21) 
p = (<j-k)F 7 + (<j-q)F 8l (22) 

with q and k the 3-momentum unit vectors of virtual photon 
and pion, respectiv ely, and F\ to Fg the CGLN amplitudes 
dChew et al ., 1957). The structures in front of the Fi are all 
the independent axial vectors and pseudoscalars that can be 
constructed from the Pauli spin matrix a and the independent 
cm momenta k and q. We further note that a and k are the 




N N* N' 



FIG. 1 Multipole notation for pion photoproduction. See text for 
further explanation. 

transverse components of a and k with regard to q. With these 
definitions F\ to F4 describe the transverse, F5 and F$ the 
longitudinal, and F7 and F$ the time-like or Coulomb com- 
ponents of the current. The latter ones are related by current 
conservation, q ■ J — uup = 0, leading to \q\F 5 = uF s and 
\q}Fe = loFj. The CGLN amplitudes depend on the virtu- 
ality of the photon, Q 2 , as well as the total hadronic energy 
W and the pion-nucleon scattering angle 6>* in the hadronic 
cm system. These amplitudes are complex functions, because 
the transition leads to a continuum state with a complex phase 
factor. They can be decom posed in a series of multipoles (see 
iDrechsel and Tiatori d 1992b for further details), 

M l ±={E l ±,M l ±,L l ±,S l ±}, (23) 

where E%±, Mi±, Li±, and Si± denote the transverse electric, 
transverse magnetic, longitudinal, and scalar (time-like or 
Coulomb) multipoles, in order. The latter two are related 
by gauge invariance, |<t|L/± = wSi±, and therefore we may 
drop the longitudinal multipoles in the following without loss 
of generality. The CGLN multipoles are complex functions 
of 2 variables, Me± — Me±(Q 2 ,W). The notation of the 
multipoles is clarified by Fig.Q] The incoming photon carries 
the multipolarity L, which is obtained by adding the spin 
1 and the orbital angular momentum of the photon. The 
parity of the multipole is V = (— l) L for E, L, and S, and 
V = (— for M, The photon couples to the nucleon 
with spin i and V = +1, which leads to hadronic states of 
spin J = I L ± \ I and with the parity of the incoming photon. 
The outgoing pion has negative intrinsic parity and orbital 
angular momentum I, from which we can reconstruct the spin 
J = I I ± \ I and the parity V = of the excited 

hadronic state. This explains the notation of the multipoles, 
Eq. d23l . by the symbols E, M, and S referring to the type of 
the photon, and by the index £± with I standing for the pion 
angular momentum and the ± sign for the two possibilities to 
construct the total spin J = \ I ± i | in the intermediate states. 

We complete the formalism of pion photoproduction by 
discussing the isospin. Since the incoming photon has both 
isoscalar and isovector components and the produced pion is 
isovector, the matrix elements take the form 

M? ± = \[r a MM^ + 1 -{r a ,r Q }M\t ) +r a M^ , (24) 

where r Q are the isospin Pauli matrices in a spherical basis, 
i.e., a = {+,0,—}. It follows that the intermediate state 
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\ or I = §. The 4 



in Fig. Q] can only have isospin I 

physical amplitudes with final states {p7r°, nir°, mr + ,pir~} 
are given by linear combinations of the 3 isospin amplitudes. 
We should however keep in mind that the isospin symmetry 
is not exact but broken by the mass differences between 
the nucleons (n,p) and pions (tt , 7T°) as well as explicit 
Coulomb effects, in particular near threshold. 



The calculation of the observables is straightforward but 
somewhat tedious, and therefore we choose pion photopro- 
duction at threshold as an illustrative example. Near threshold 
the partial wave series may be truncated to s and p waves, i.e., 
the transverse multipoles E a +, M\+, Ei+, and M\-. With 
Pi = 3E 1+ + M 1+ - M x -, P 2 = 3E 1+ - M 1+ + Mi-, and 
P3 = 2Mi+ + Mi- the differential cross section takes the 
following form in the cm frame: 



(25) 



with A 
and C 



±|P 2 | 2 + ±|P 3 | 



B = 2Re(E* + P 1 ), 



1| ftl 2 . As is to be expected, 



l^o+ 1 

I-P1I 2 ~ - 5K3 

the s-wave multipole yields a constant angular distribution, 
whereas the forward-backward asymmetry is given by the in- 
terference between the s wave and the p-wave combination 
Pi. The terms in cos 2 #* determine a further p-wave combi- 
nation, I P2 1 2 + I P3 | 2 . A complete experiment requires to 
measure one further observable, the photon asymmetry 



E(0J) 



do- 1 - daW 



da^ 

-( 

2q V 



■daW 



P 2 | 2 -|P 3 | 2 )sin 2 i 



'da 

dn { 



(26) 



where _L and || stand for photon polarizations perpendicular 
and parallel to the reaction plane. 

The theory of meson electroproduction is more in- 
volved and we re f er the reader to the literature, see 
iDrechsel and Tiatorl (1 1992b and references quoted therein. 
The scattered electron serves as a source of virtual photons 
whose flux IV and transverse polarization e can be controlled 
by varying the electron kinematics. Moreover, we assume that 
the electron beam is polarized. The five-fold differential cross 
section for meson electroproduction is written as the product 
of a virtual photon flux factor T v and a virtual photon cross 
section, 



da 



dVL2d€2d£l%- 



da 



(27) 



The electron kinematics is commonly given in the lab system, 
whereas the hadrons are described in the hadronic cm system 
as indicated by an asterisk. The reaction plane and the elec- 
tron scattering plane have the same z-axis, but the former is 
tilted against the latter by the azimuthal angle ^ 7T . With these 
definitions, the virtual photon cross section takes the follow- 



ing form for polarized electrons but unpolarized hadrons: 



da 
dttt 



daj 



df_L 



—r da LT , 
el — — — cos® 
' dfli 



daTT „ / — ; r da t.t> 



(28) 
sin $„■ . 



Denoting the initial and final electron lab energies by ei and 
€2, respectively, the photon lab energy is u>l — ei — £2, and in 
the same notation the photon lab three-momentum is given by 
<7l. With these definitions the transverse electron polarization 
and the virtual photon flux take the form 



1 



1 



nil 



tan 2 



2tt 2 ei Q 2 



1 



1 



(29) 



with K = (W 2 — M 2 )/2M the photon "equivalent energy" 
in the lab frame. The partial cross sections in Eq. d28l ) are 
functions of the virtuality Q 2 , the pion scattering angle 0*, 
and the total hadronic cm energy W. The first two terms 
on the rhs of this equation contain the transverse (jjt) and 
longitudinal (ox) cross sections. The third and fifth terms 
yield the longitudinal-transverse interferences a^x and a^T' ■ 
These terms contain the explicit factors cos^Tr and sin*^, 
respectively, and an implicit factor sin #* in the partial cross 
sections, and therefore they vanish in the direction of the vir- 
tual photon. The latter is also true for the fourth term, which 
contains the transverse-transverse interference (ottX which 
is proportional to sin 2 #* and appears with the explicit factor 
cos2$ w . These 5 partial cross sections can be expressed in 
terms of the 6 independent CGLN amplitudes Pi to F$, or in 
terms of 6 helicity amplitudes Hi to H§ given by linear com- 
binations of the CGLN amplitudes. The particular form of the 
t&jr dependence in Eq. d28l ) is of course related to the fact that 
the virtual photon transfers one unit of spin. A close inspec- 
tion shows that the 5 responses provided by the polarization 
of the electron can be separated in a "super Rosenbluth plot". 
This requires measuring the polarization e of the virtual pho- 
ton, the beam polarization P e , and the angular distribution of 
the pion with at least one non-coplanar angle <lv. A double- 
polarization experiment measuring also the target or recoil po- 
larizations of the nucleon yields 18 different response func- 
tions alt ogether. The relevant expre s sions can be found in th e 
work of IDrechsel and Tiatori dl992b ; iKnochlein et (119951) . 



E. Resonance excitation 

As shown in the previous section, each partial wave is char- 
acterized by 3 quantum numbers, orbital angular momentum 
I, total angular momentum J, and isospin I. Most of these 
pion-nucleon partial waves show distinct resonance structures 
at one or more values of the hadronic cm energy. Furthermore, 
there are generally 3 (independent) electromagnetic transi- 
tions between the nucleon and a particular partial wave, an 
electric, a magnetic, and a Coulomb transition. Let us con- 
sider as an example the most important resonance of the nu- 
cleon, the A(1232) with the spectroscopic notation P33, which 
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decays with a life-time of about 0.5 x 10 s into a pion- 
nucleon state, except for a small electromagnetic decay branch 
of about 0.5 %. This intermediate state contains a pion in a p 
wave, i.e., t = 1 and V = +1. The indices 33 refer to twice 
the isospin and spin quantum numbers, J = J = S. The elec- 
troexcitation of this resonance takes place by magnetic dipole 
(Ml), electric quadrupole (E2), and Coulomb quadrupole 
(C2) radiation, which is denoted by the 3 complex func- 



3/2 3/2 

tions Mi| , EZ{ , and S^. If one neglects the small pho- 
ton decay branch, unitarity requires that all 3 electroproduc- 

tion multip oles carr y the p hase 5^(W) of the pion-nucleon 
final state (Watsonl 11954 . For a st able particle with the 
quantu m numbers of a P33 resonance, Rarita and Schwinger 
(Il94ll) have developed a consistent relativistic theory involv- 
ing 3 (real) form factors G* (Q 2 )- However, the physical pion- 
nucleon state has a complex phase factor, the resonance phe- 
nomenon spreads over more than 100 Me V in excitation en- 
ergy, and there is no model-independent way to extract the 
"bare" resonance parameters from the observables. It is there- 
fore common practice to relate the form factors to the transi- 
tion multipoles taken at the resonance position, W — Ma = 
1232 MeV. Corresponding to the independent transition multi- 
poles, the following 3 form factors for the N A transition have 
been defined: 



,3/2 



<(M A ,Q 2 ) 
<(Ma,Q 2 ) 



= -iN 



-iN 



M 
9a(0 2 ) 



M 
1a(Q 2 ) 2 
2 MM A 



G* e (Q 2 ), (30) 



G* C (Q 2 ). 



with <7a and &a the 3 -momenta of photon and pion at the 
resonance, and N = \f"ia em /{%k\T \) a kinematic factor 
relating pion photoproduction to total photoabsorption at 
the resonance. We note that these definitions divide out 
the q dependence of the multipoles at pseudothreshold 
(q — * 0) such that the form factors are finite at this point. 
Equati on (f30b corresponds to t he definition of lAsh et al\ 
(1 1967b . the form factors of IJones and Scadronl d 1973b 
are obtained by multiplication with an additional factor, 
G JS = ^l + Q 2 /(M + M A y G Ash . 



Because the background becomes more and more important 
as the energy increases, the concept of transition form fac- 
tors is usually abandoned for the higher resonances. Instead 
it is common to introduce the helicity amplitudes, which are 
uniquely defined for each resonance by matrix elements of the 
transition current between hadronic states of total spin J and 
projection M. With the photon momentum q as axis of quan- 
tization, the virtual photon can only transfer intrinsic spin 1 to 
the hadronic system, with projections ±1 for right- and left- 
handed transverse photons (current components J±i) and 
for the Coulomb interaction (time-like component p). Using 
parity and angular momentum conservation, we find 3 inde- 



pendent helicity amplitudes, 
1 



/2 



.4 



3/2 



Si/2 



(N*(J,±)\J+\N&-±)), 
(N*(J,I)\ J+|iV(|,±)), (31) 
(N*(J,±)\p\N(l,±)). 



In particular we note that the amplitude A 3 / 2 exists only 
for resonances with J > |, and neither does this amplitude 
exist for a free quark. Hence asymptotic QCD predicts that 
A 3 / 2 should vanish in the limit of large momentum transfer, 
Q 2 — ► 00. The electromagnetic multipoles can be expressed 
by combinations of the helicity amplitudes. For the A (1232) 
these relation take the following form: 



M 



E 



S 



3/2 
1+ 

3/2 
1+ 

3/2 
1+ 



1 



2^3 
1 

2\/3 



(V3Ai /2 
(V3Ai /2 



3 A 3/2 ) . 
A3/2) , 



(32) 



"7i 5l/2 



It is interesting to observe that asymptotic QCD predicts the 
following multipole ratios in the limit Q 2 — > 00: 



Rem — 



Im E 



3/2 



ImM 



3/2 



Rs 



Im S 



A I 



3/2 
1+ 



Im Ml 



3/2 



W=M A 



W=Ma 



const. 



(33) 



In these relations, the multipoles are evaluated at resonance, 
defined by the energy for which the real part passes through 
zero (K-matrix pole). This should happen at the same energy 
for all 3 multipoles as long as the Fermi- Watson theorem is 
valid. 



F. Dispersion relations 

Dispersion relations (DRs) play an important role in the 
following sections. They are based on unitarity and analyt- 
icity and, by proper definitions of the respective amplitudes, 
fulfil gauge and Lorentz invariance as well as other symme- 
tries. The analytic continuation in the kinematic variables 
allows one to connect the information from different physical 
processes and thus to check the consistency of different sets 
of experiments. As demonstrated in section [IV] DRs are 
prerequisite to determine the polarizabilities of the hadrons 
from Compton scattering, and several sum rules are derived 
in section [Vl] by combining DRs with low-energy theorems. 
Most of these techniques are very involved and we have to 
refer to the literature. Therefore we only give an overview of 
the dispersive approach for the nucleon form factors, which 
are discussed in more detail in the following section llHl 



Let G(t) be a generic (electromagnetic) form factor describ- 
ing the ground state of the nucleon. The real and imaginary 
parts of G(t) are then related by DRs. Assuming further an 
appropriate high-energy behavior, these amplitudes fulfill an 
unsubtracted DR in the Mandelstam variable t, 



0.06 



Re G{t) 



1 



d£ 



Im G(f) 



if - t 



(34) 



where to is the lowest threshold for the electroproduction of 
pions by e + e - pair annihilation. These form factors can 
be measured by electron scattering for space-like momentum 
transfer (t = — Q 2 < 0) and by collider experiments for time- 
like momentum transfer (t > AM 2 ). The imaginary part or 
spectral function, Im G(t), vanishes in the space-like region, 
and therefore the ie in Eq. (|34l can be dropped for elastic elec- 
tron scattering. However we note that Eq. (|34| > defines the 
real part of the form factor in both the space-like and time- 
like regions, provided that the spectral function is sufficiently 
well known. The dispersive formalism also yields informa- 
tion on proton and neutron at the same time as is evident 
from the following reasoning. The spectral function can be 
obtained from the two-step process 7* — * X — ► N N, with 
X a hadronic state with the quantum numbers of the photon. 
In the usual notation of these quantum numbers with isospin 
I, G-parity, spin J, parity P, and C-parity, the isoscalar 
photon has I G {J PC ) = CP(1 ) and the isovector photon 
I G (J PC ) = 1 + (1 ). The lightest hadronic system X in 
the intermediate state is a pion pair, which has even G-parity 
and therefore contributes only to the isovector current. This 
part of the spectral function is therefore composed of (I) the 
vertex 7* — > it tt given by the pion form factor F^if), here in 
the time-like region and therefore a complex function, and (II) 
the process it tt — > N N. The latter process is needed in the 
unphysical region, which can however be reached by analytic 
contin uation of the p -wave amplitudes for pion-nucleon scat- 
tering dHohlerl [l983). As a result, the two-pion contribution 
to the spectral function can be constructed from to = Am 2 up 
to about 1 GeV 2 as 

lmG%{t) = A^F^tffHt), 



AlVt 



Im Gl.it) 



9t 



FAt)*f-(t), 



(35) 



with q t = y/t/A — m 2 the pion momentum in the interme- 
diate state and f±(t) the p-wave im — * N N amplitude. 
The spectral functions for the Sachs form factors are plot- 
ted in Fig. [2] The figure shows a rapid rise of the spectral 
functions at the two-pion threshold (t = A m 2 ), because the 
projection of the nucleon Born graphs to the p wave yields 
a singularity on the second Riemann sheet at t = 3.98 to 2 , 
just below threshold. Quite similar results for the two-pion 
continuu m have also be en obtained by a two-loop calculation 
in ChPT dKaiserll2003l) . Furthermore, we observe the strong 
peak near t w 28 m 2 , which is due to the p meson with mass 
770 MeV and a large width. The spatial distribution of charge 
and current can be obtained by the respective form factors in 
the Breit frame. The starting point is Eq. (l34l for space-like 




FIG. 2 The isovector spectral functions in units of m n . Solid 
lines: 2lmG v E (t)/t 2 , dashed lines: 21mG V M (t)/t 2 . The thin lines 
are from the work of Hohler and Pietarine n M 19751) . the thick lines 
include modern data f or the pion form factor. The figure is from 
teelushkin et aZ.Ll2006h . 



t = —Q 2 — > — q 2 , which is Fourier transformed to r-space 
with the result 



p{r) = i 



dtlmGit)- 



-Vtr 



(36) 



in 



The mean squared radius for a particular region of the spec- 
tral function at t = /1 2 , where fi is the mass of the intermediate 
state, is given by < r 2 >= 6//1 2 . For instance, the onset of 
the spectral function corresponds to an rms radius of about 
1.7 fm, the p meson to about 0.6 fm, and so on. We conclude 
that the density at large distances is dominated by the lightest 
intermediate states. As a consequence, the tail of the density 
distribution at large radii should take a Yukawa form, e _A ""/r, 
with p v = 2™^ and /i s = 3771^ for the isovector and isoscalar 
form factors, respectively. It is therefore natural to identify 
the "pion cloud" with the two-pion contribution to the spectral 
function, which remains after subtraction of the p peak from 
the spectral function of Fig. [2] Whereas the isovector spectral 
function can be constructed from the available experimental 
information up to t w 1 GeV 2 , the higher part of the spectrum 
has be modeled from information about resonances and con- 
tinua. Because the isoscalar spectral function contains at least 
3 pions in the intermediate state, it can not be obtained directly 
from experimental data. In the region below t m 1 GeV 2 it is 
dominated by the to and $ mesons, the three -pion continuum 
has be en shown to couple only weakly, see Belushk in et al.\ 
(2007) and references to earlier work. 



FORM FACTORS 



Ever since iHof stadterl d 1 956L 1 1 957h first determined the size 
of the nucleon, it has been taken for granted that the nu- 
cleon's electromagnetic form factors follow the shape of a 
dipole form, with some minor deviations and, of course, mod- 
ified for the vanishing charge of the neutron. This form was 
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FIG. 3 The world data for the 4 nucleon form factors according to iFriedrich and Walched l2003h . The solid line is the phenomenological 
fit given by Eqs. d38t and l |39l >. The other lines sho w the two dipole contribution s and the bump/dip term of the fit separately. Note that the 
absolute value is plotted for negative quantities. See Friedrich and Walcher ( 2003) for a detailed listing of the data. 



conveniently parameterized as 

G(Q2) ° a +Q Wr 1371 

with Ajj w 0.84 GeV a universal parameter. Because this 
parameter is close to the mass of the p meson, it was assumed 
that the nucleon structure is dominated by a vector meson 
cloud which was described by the "vector dominance model". 
This idea was of course in conflict with the quark model after 
its establishment in the 1970's and many attempts were made 
to reconcile these conflicting models. 

In order to set the scene, let us recall the following proper- 
ties of proton, neutron, and heavier baryons: 

• The complexity of these strongly interacting many- 
body system reflects itself in the finite size in space, 
the anomalous magnetic moment, and the continuum of 
excited states with strong resonance structures. These 
three aspects are of course closely related, and can be 
connected quantitatively in some cases by sum rules as 
detailed in sectionlVTl 

• Because of the approximate SU(3) symmetry of u, d, 
and s quarks, the nucleon forms a doublet in isospin 



with strangeness S = in an octet of states. The other 
partners are: the A and the triplet with 
S = -1, and the with S = -2. The strange 

baryons decay weakly with a typical mean life of about 
10~ 10 s. Because the E° can also decay into the A by 
photoemission, its mean life is only about 10 -19 s. The 
most important resonance of the nucleon, the A(1232) 
appears as an isospin quadruplet with strangeness S = 
in a decuplet of states. The partners of the A(1232) 
are: an isospin triplet £+(1385), E°(1385), £"(1385) 
with S = -1, a doublet 3°(1530), H"(1530) with S = 
-2, and the ft" (1672) with S = -3. The latter lives 
about 10~ 10 s, because it can only decay weakly. All 
the other particles in the decuplet decay by the strong 
interaction with a mean life of order 10~ 23 s. 

• The size effect reflects itself in the form factors of elas- 
tic electron scattering. Because of the life time, only the 
proton and, with some caveat, the neutron can be stud- 
ied as a target. With the dipole form of Eq. d37] > and 
the definitio n of an rms radiu s according to Eq. ( fT4l . 
the result of iHofstadterl d 1 956b was r p E = 0.81 fm for 
the charge distribution of the proton. In the mean time 
several new experiments have led to the larger radius of 
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r E ps 0.88 fm. The form factors of the strange baryons 
can in principle be measured by scattering an intense 
beam of these particles off the atomic electrons of some 
nuclear target. Such an experiment was performed 
by the SELEX Collaboration at the F ermi Lab Teva- 
tron w ith a high-energetic £~ beam (Es chrich et all 
120011) . The result is a first datum on a hyperon radius, 
r|" = (0.78 ± 0.08 ± 0.05) fm, distinctly smaller than 
the accepted value for the proton. 

• The "normal" magnetic moment of a particle i expected 
for a pointlike Fermion is given by eQi/(2Mi), with 
Mj the mass and e Qi the charge of the particle. If the 
magnetic moments /ij are given in units of the nuclear 
magneton [/ijv] = e/(2M p ), the values for proton and 
neutron, \i p = 2.79 \fijf] and fi n = —1.91 [/xjv], sig- 
nal a large isovector anomalous magnetic moment of 
the nucleon. From electron scattering we also know 
that the rms radius of the magnetization distribution 
is quite similar to the radius of the charge distribu- 
tion. Because of their long mean life, the magnetic 
moments of 5 other octet baryons and, in addition, the 
one of the f2~(1672) in the decuplet are known from 
spin precession experiments. Without going in the de- 
tails, also these particles have large anomalous mo- 
ments. As an example, even the f2~, a configuration 
of 3 s quarks with the large mass of 1.672 GeV, has 
fin = —2.02 [/j,n] compared to a "normal" magnetic 
moment of —0.56 [jUjv]- In order to get more infor- 
mation about the decuplet, several experiments were 
performed to measure the magnetic moment of the 
A ++ as a subprocess of radiative pion-nucleon scat- 
tering dBosshard et aZ.Ul99ll:lNefkens et a?lll978l) and 
of the A + as a subpro cess of radiativ e pion photopro- 
duction on the proton dKotullai 120031) . The results for 
the magnetic moments are in qualitative agreement with 
quark model predictions but still with l arge model errors 
dPascalutsa and Vanderhaegheni 120071) . 

At the turn of the century new surprising results put the nu- 
cleon form factors into focus once more. These new results 
became possible through the new generation of cw electron 
accelerators with sources of high-intensity polarized beams 
combined with progress in target and recoil polarimetry. As 
summarized in section ITXBl the measurement of asymmetries 
allows one to determine both form factors even if they are 
of very different size. This situation occurs in two cases: 
(I) Because of its vanishing total charge but large anomalous 
magnetic moment, the neutron's electric form factor is very 
much smaller than the magnetic one, at least for small and 
moderate momentum transfer. (II) As shown by Eq. ( fl~5b . 
the magnetic form factor Gm(Q 2 ) appears with a factor 
t = Q 2 /4M 2 , whereas Ge{Q 2 ) is suppressed by a factor 
1/(1 + t). As a consequence, Ge(Q 2 ) is less well deter- 
mined by the Rosenbluth plot if Q 2 becomes large. Even 
though this was known, it was to the great surprise of ev- 
erybody when asymmetry measurements showed a dramatic 
deviation from previous results based on the Rosenbluth sep- 
aration and, at the same time, from the dipole shape of the 



proton form factors dGavou et al. I l2002t IJones et all [2000). 
Another open question concerns the behavior of the form fac- 
tors at low 4-momentum transfer, such as oscillations at very 
small Q 2 and conflicting results for the rms radius of the pro- 
ton. All these experimental findings have caused an intense 
theoretical investigation t hat has been summarized by sev- 
eral recent review papers (jArr ington efadl2007bt lGaoT.12 003 ; 
iHvde- Wright and de Jagerl l2004t iPerdrisat et all 120071) . In 
the present work we concentrate on the low momentum trans- 
fers and, therefore, the phenomena for Q 2 > 1 GeV 2 will only 
be discussed briefly. 



A. Space-like electromagnetic form factors of the 
nucleon 

The new results in the field of space-like form factors 
have been obtained at basically 3 facilities, the cw e l ectron 
accelerators CEBAF at the Jefferson Lab (ICardmanL [2006f) 
and Mainz Microtron dJankowiakl 12006). and the electron 
stretcher ring at MIT/Bates dMilneil 120061) . All these facilities 
provide an intense beam of polarized electrons. The second 
essential for measuring the asymmetries was the develop- 
ment of polarized targets and polarimeters to determine the 
polarization of the recoiling particles. For details of the 
new accelerators, targets, and particle detectors as well as 
pe rtinent references, we refer the re ader to a recent review 
of IHvde- Wright and de Jagerl d2004 . In this context it is 
however important to realize that there exist no free neutron 
targets, but only targets with neutrons bound in a nucleus. 
Therefore, any analysis of the data requires theoretical models 
to correct for the binding effects, which include initial-state 
correlations, meson exchange and other two-body currents 
with intermediate resonance excitation of the nucleons, and 
final-state interactions while the struck nucleon leaves the 
target. In this situation, we infer that the deuteron provides 
the most trustworthy neutron target, because its theoretical 
description is far more advanced than in the case of heavier 
nuclei. Of course, measurements with heavier nuclei, in 
particular polarized 3 He targets, provide complementary 
information and are interesting for their own sake. 

In Fig. [3] we display the nucleon form factors as functions 
of Q. We choose this somewhat unusual presentation in or- 
der to e mphasize the small Q 2 region . The data base shown 
is from iFriedrich and Walcheri d2003l) . which meanwhile has 
been complemented by results from the following references : 



been complemented by results rrom tne rollowing rererences 
dAnderson et all 2007tlBermuth et all 120031: Crawford et al. 



120071: iGlazier et all 120051: iPlaster et all 120061: IZiskinl 120051) 



The phenomenological fit shown by the solid line in Fig. [3] 
is composed of two dipoles and a bump/dip structure. The 
dipole form is given by 



G S (Q 2 ) = 



aw 



»20 



(l + Q 2 /a n ) 2 (l + Q 2 /a 21 ) 2 ' 



(38) 



and the bump/dip structure, seen at Q 2 — Q 2 ss 0.2 GeV 2 on 
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4 The nucleon form factors G P E , G v M j \i v , and Gjjf / fi„ as func- 



tion of Q, divided by the stan dard dipole form factor of Eg. fl37b . The 
full line represents the fit of iFriedrich a nd Walchej] J2003h . the bro- 
ken line shows the smooth contribution of the two dipoles. 



top of the smooth dipoles, is parameterized as 



G b (Q 2 ) = a b Q 2 



.If' 
2 V 



h\2 



1( Q+Qb \2 
2\ a b > 



(39) 



We note that this ansatz provides an even function of Q as 
required by g eneral argu ments. A similar form has been 
introduced by ISickl dl974l) in r-space in order to obtain a 
non-singular function for his model-independent analysis of 
nuclear charge distributions. 

Figure [3] shows that the bump/dip structure is of the or- 
der 3% and only visible for Q 2 < 1 GeV 2 . Therefore we 
have magnified the data and their structure by the linear plot 
of Fig. 2] which shows the form factors divided by the stan- 
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squar es: d(e, e'n)p 



FIG. 5 The world data for the electric 
neutron, G r, as functio n of Q 2 . Solid 
at MAMI, iGlazieref al\ d2005h; lOstrickl d2006h . solid trian- 
gles: d{e, e'n)p at Jlab, Ma dev et al\ 1(2003), solid diamond: 
d(e, e'n)p at BATES . lEden et al\ il994l) . solid circle: d(e, e'n)p at 
NIKHEF, Passchi er et al\ jl999h . open triangles: d(e, e'n)p at Jlab, 
Warre n et al. I J2004h , open diamonds: die, e'n)v at BATES JZiskinl 
j2005h . open squares: 3 He(e,e'n)pp at MAMI. lRohl i2006h . The 
open circles are der i ved fr om the deuteron quadrupole form factor 
(Schiavilla and Sick, 200lh. Solid line: new fit of the phenomeno- 
logical model of Fried rich and Wa lcher (2003), da shed line: the re - 
sult of dispersion relation {gelushkin et ai, 2007; MeiBner. l2007h . 
dotted line: the original Galster fit of Eq. < |40t . The new fit includes 
all the shown data except for the preliminary o flZiskinl J2005) and the 
values derived by Schiavilla and Sick (2001). The data in the figure 
were updated by J. Friedrich and the theory curves complemented by 
L. Tiator. 



dard dipole form factor of Eq. d37] i. The electric form fac- 
tor of the neutron is quite special because of its vanishing 
charge, which results in an overall small value of G%. There- 
fore we have plotted this form factor in a different way in 
Fig. |5J which displays the published world data as measured 
with polarized elect r ons. W e note, however, that the results of 
ISchiavilla and Sickl d200ll) have been deduced from an anal- 
ysis of the deuteron quadrupole form factor Fc2, which re- 
quires a careful investigation of the model dependence due to 
the nucleon-nucleon potential. The error bars of these data are 
therefore not statistical but indicate the (systematic) model er- 
ror. The combined data shown in Fig. [5] clearly support the 
existence of the bump structure at Q 2 w 0.2 GeV 2 as in the 
previous cases. The solid line in this figure is the result of a 
new fit with the phenomenological model given by Eqs. d38l ) 
and ( f39b . Th e dashed line in this f igure is the parameterization 
first given bv lGalster et al\ dl97ll) . 



G n E {Q 2 ) - 



a G t 



1 



(l + b G T) (l + QVA^ 



(40) 



with a G = — Hn an d b G = 5-6. The result from dispersion 
theory is displayed by the dotted line. Neither dispersion 
theory nor the Galster fit reproduce the data at low Q 2 . 
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NEUTRON CHARGE DENSITY 
visualization of different ideas 
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FIG. 6 The charge density distribution for the neutron, 4%r p%, as 
function of r. The notation for the lines corresponds to that of the 
form factors shown in Fig. [5] Figure by courtesy of J. Friedrich and 
L. Tiator. 



The Fourier transform of G E is of particular interest, be- 
cause the overall charge of the neutron must vanish. A finite 
charge distribution is therefore a definite sign of correlations 
among the charged constituents, for example, between the u 
quark and the two d quarks of the constituent quark model. 
The charge distribution p r E is displayed in Fig.|6]for the 3 fits 
to the neutron factor shown in the previous figure. We recall 
the arguments of subsection III.BI that the Fourier transform 
can only be obtained in the Breit or "brick wall" frame, which 
however is a different Lorentz frame for different values of 
Q 2 . The visualization of the Fourier transforms as charge 
and magnetization distributions in r-space is therefore only 
approximately correct if the momentum transfer is small 
compared to \Q 2 \ — AM 2 , which defines the threshold of 
nucleon pair production at time-like momentum transfer. 
With this caveat in mind, we may interpret the charge 
distribution of Fig. [6] by the dissociation of the neutron in 
a proton and a pion, i.e., as a negative "pion cloud" around 
a positiv e core. As we see from the figure, the pion cloud 
found bv lFriedrich and Walcherl d2003l) extends to large radii. 
It is important to realize that this result does not depend 
on a model assumption but is borne out by a statistically 
satisfactory reproduction of the data. We also note that the 
signal of the pion cloud is empirically present in all 4 form 
factors. 



In a recent paper. lMillerl d2007l) has found a negative density 
in the center of the neutron from an analysis of generalized 
parton distributions as function of the impact parameter b. The 
charge distribution p\ (b) is then defined by the 2-dimensional 
Fourier transform of the Dirac form factor F[ L (Q 2 ). We note 
that this does not contradict the results shown in Fig. [6] In fact 
our results for the 3 -dimensional Fou rier transform of ^"(Q 2 ) 
agree very much with the findings of Miller ( 2007): a negative 
density in the center, positive values for 0.4 fm< r < 1.2 fm, 
and a negative tail for the larger distances. 
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FIG. 7 The world data for the time-like form factor of the 
proton as function of the beam energy Ekeam, extracted with 
the assumption \Ge\ = I Cm I- Open symbols : data from 
e+e" -» vv jAblikim et all 120051: lAntonelli et all 1 19961, fl998 



Aubert ef a/.L 120061; Bisello et all 119831 Il99d ; ICastellano et al. 
19731 : iDelcourt et allll979l;|Pedfar a/.ll2005t). solid symbols: data 
from pp — > e~e + IXm brogiani et al., 1999; Andreotti et al j [20031; 
Armstrong et all 1 19931 ; iBardin et all 1 19941 : fBassompier re et all 
19771) . The figure is from fRossi et aZ.ld2006ri . 



B. Time-like electromagnetic form factors of the nucleon 

Let us next discuss the form factors for time-like momen- 
tum transfer, i.e., positive values of the Mandelstam variable 
t = q 2 = —Q 2 > 0, see subsection III. Al for definitions. By 
inspection we find that the previously defined dipole form fac- 



tors have poles t 



My, that is, the phenomenological 



fits "predict" the existence of vector mesons in the time-like 
region. The form factors in the space-like and time-like re- 
gions are connected by analyticity and unitarity, and therefore 
also the knowledge of the time-like form factor s is mandatory 
for a c omplete understa n ding o f the nucleo n (Baldin i et al. 



19991: iBelushkin et all 120071: iGeshkenbein et all 1 1974 



Hammed 120061: iMergell et all 1 19961) . The time-like pho- 
tons are obtained in collider experiments by the reaction 
e+er -> N N for t = (2E be . dm ) 2 > 4Af 2 . Whereas the 
space-like form factors are real, the time-like form factors 
are complex functions because of the strong interaction 
between the produced hadrons. However, the unpolarized 
cross section in the time-like region only depends on the 
absolute values of the two form factors, | Ge (t ) | and | Gm (t) \ ■ 
In order to get information on the relative phase between 
the form factors, polarization experiments are required. 
Unfortunately, the present data basis does not even allow for 
a Rosenbluth separation. Therefore the data are analyzed 
with the assumption |G£:(i)| = |Gm(^)|, which follows from 
Eq. ( fT3l at threshold, t — AM 2 , but is of course not expected 
to hold for higher beam energies. 
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FIG. 8 The world data for the time-like form factor of the neutron 
as function of the beam energy Eteam, extracted with the assumption 
\G E \ = \G M \- The data are obtained by the FENICE collabora- 



tion from the reaction e + e 



The s olid line is a perturba- 



tiveQCD extrapolation usi ng the proton dat a jAntonellief~a?lll994 
Il998h . The figure is from Ros si et al\ d2006r) . 



Figure [7] displays a compilation of the proton data known 
so far. These data cover the range of 4M 2 < t < 6.8 GeV 2 . 
The figure shows an overall falloff with the beam energy for 
G P M , somewhat faster than 1/t 2 , and some structure near 
E heam = 1.05 GeV or t = 4.4 GeV 2 , which may indicate 
a resonance in that region. We note that the decrease of 
the form factor is in qualitative agreement with perturbative 
QCD, which requires a falloff like \t\~ 2 for both space-like 
and time-like photons. However, a comparison shows that 
the space-like form factor at t = —Q 2 = —16 GeV 2 is al- 
ready about a factor 3 smaller than the time-like form factor 
at t = +16 GeV 2 , that is, asymptotia is still far away. A look 
at Fig.[8]tells us that our knowledge about the neutron's time- 
like form factors is still far from satisfactory. We hope that 
the currently planned experiments will improve on the preci- 
sion of the time-like form factors and, in particular, also de- 
termine their relative phases, which is absolutely necessary in 
order to get the full information on the structure of the nucleon 
dRossi et ■a/.Ll200d) . 



C. Theoretical considerations 

The electromagnetic form factors encode information on 
the wave functions of the charged constituents in a bound sys- 
tem. However, in the case of the hadrons we face severe obsta- 
cles to get a real grip on the elementary quarks. As has been 
mentioned in section [TJ only two ab-initio approaches exist to 
describe QCD in the confinement phase, chiral perturbation 
theory and lattice gauge theory. Chiral perturbation theory is 



restricted to small values of the momenta. Moreover, if ex- 
tended to higher order in the perturbation series, ChPT loses 
predictive power, because the number of unknown low energy 
constants increases. Lattice gauge theory, on the other hand, 
is still hampered by the use of large quark masses. This has 
the consequence that the pionic effects appearing at low mo- 
mentum transfer are underestimated. Beyond these two ap- 
proaches, which are in principle exact realizations of QCD, a 
plethora of "QCD inspired" models with quarks and pions has 
been developed. The problems are twofold: 

• Starting directly from QCD, one would have to use 
the small u and d quark masses of order 10 MeV. The 
many-body system is therefore highly relativistic from 
the very beginning. However, a typical constituent 
quark model (CQM) has quarks with masses of sev- 
eral hundred MeV. It is therefore obvious that these en- 
tities are many-body systems of quarks and gluons by 
themselves. In any case, the constituent quarks wave 
functions have to "boosted" if hit by the virtual pho- 
ton. However, there exists no unique scheme to boost a 
strongly interacting relativistic many-body system. 

• In view of the small current mass of the quarks, the in- 
teraction as mediated by gluon exchange inevitably pro- 
duces a considerable amount of quark-antiquark admix- 
ture. These effects have to be modeled by properties 
of the constituent quarks, such as mass and form fac- 
tor, see lDe Sanctis et aZl d2005al) . or by explicitly intro- 
ducing a meson clou d of the "bare" constituent quarks 
dFaessler et adl2006l) . 

Since this article is dedicated to the low-Q 2 domain, it 
suffices to consider some models useful in this region. A 
more det ailed discussion of the wide range of models is 
given by iPerdrisat et al.\ d2007l) . The traditional model of 
the nucleon is the CQM with quark masses m q w M/3. 
Except for the smallest momentum transfers, the quark 
wave functions have to be relativized, which is usually done 
by relativistic boosts of the single-quark wave functions. 
Figure [9] shows the result of several calculations for the ratio 
G P E /G P M compared to the recent data from the Jefferson Lab 
obtained with a double-polarization experiment. The rapid 
falloff of this ratio was a real surprise, because previous 
experiments without polarization did not find big deviations 
from the dipole fit for both form factors. The solution of 
problem was explained by two-photon effects, which are 
usu ally small but turn out to be c ome l arge i n special cases, 
see Carlson and Vand erhaeghen (2007) and Arringto n et al.\ 
(l2007al) for recent reviews. Even though the figure shows 
only a small selection of models, quite different results are 
obtained by similar models, depending on the properties of 
the constituent quarks, their interaction, and the boosting 
mechanism. We are also not aware of many predictions for the 
rapid drop of the ratio, before the JLab double-polarization 
data were obtained. In any case, the shown models describe 
the data qualitatively well, and they support the expected 
zero crossing of the electric form factor at Q 2 w 7 GeV 2 . 
The zero crossing is also predicted from a Poincare covariant 
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Q a (GeV 2 ) 

FIG. 9 The ratio of the electric and magnetic form factors 
of the proton, fi p G p E /G p K[ , compared to relativistic calculations 
in the framework of the CQM . Dotted line: light-cone front 
CQM (Chung and CoesteJ 1 19911) . thick solid line: light-cone front 
CQM ( Fran k et aZl Il996h . dashed-dotted line: light-cone front 
CQM with point-like consti tuent quarks (Cardarelli et~ail 1 19951 : 
Cardar elli and Sim ula. 2000), dashed line: G oldstone bo s on ex - 
change between point-like constituent quarks ( Boffi et al., 2002), 
thin solid line: covariant spectator CQM (Gross and Agbakpe. l2006l) . 
The data hav e been taken a t Jeffe rson Lab bv lPuniabi etali (2003) 
(circles) and iGavou et al\ (20021) (squares). The figure is from 
iPerdrisatef al.\ (20071) . 



Faddeev calculation describing the nucleon as a correlated 
quark-diquark system, gen erated with an interaction fitted 
to the structure of mesons (lAlkoferef all l2005t iHoll et all 
2005). On the other hand, the models without explicit 
Goldstone bosons, in particular pions, can not describe the 
region of Q 2 J$ 0.5 GeV 2 , in which the pionic degrees of 
freedom play a decisive role. 

The inclusion of a pion cloud into quark models of the 
nucle on started with the "little bag model" dBrown and Rhol 
119791) and was first applied to the nucleo n form factors 
in the form of the "cloudy bag model" dLu et all Il998t 
iMilleiL 120021: iMiller and Frank! 120021: iThomasl 1 1984b . More 
recently, iPasquini and Boffil d2007l) studied a system of 
valence quarks surrounded by a meson cloud with light-cone 
wave functions. They found distinct features of such a cloud 
below Q 2 = 0.5 GeV 2 , however without a pronounced 
bump-dip structure. Another model incorporating quarks 
and Goldstone boson s is th e "chiral quark soliton model" 
dDiakonov and Petrovl Il986t iDiakonov eiail Il988l) . The 
nu cleon form factors we re calculated within this framework 
by IChristov et al. (1995). The linear decrease of the ratio 
G V eIG v m with Q 2 was shown to follow f rom this model quite 
naturally dHolzwarfhl 1 19961 120021. 120051) . However, in all of 
these models the bump/dip stru cture of the form fac tors is not 
in focus. On the other hand, iFaessler et al. I d2006l) showed 
that this structure can be reproduced within a chiral quark 
model by a cloud of pseudoscalar mesons. One should also 
keep in mind that usually the authors do not give a complete 



description of the data. For example, a cross section ratio may 
be obtained in agreement with the data, whereas the model 
fails to describe the indivi dual cross sections. On the c ontrary, 
the parameterization of iFriedrich and Walcherl d2003l) covers 
the full Q 2 domain up to 7 GeV 2 , and therefore the charge 
distribution derived from this work is directly based on the 
experimental data. As is obvious from Fig. [6] the neutron 
charge distribution (7-), as defined by the Fourier transform 
of the electric Sachs form factor C^(Q 2 ), is positive in the 
interior region and negative for radii larger than about 0.7 fm. 
However, it certainly takes a model to quantify the separation 
into components, say a core and a "pion cloud". As an 
example, there is no unique way to break the spectral function 
of Fig.|2]into parts belonging to the two-pion continuum and 
heavier intermediate states like the p meson. On the other 
side, it is also evident that the tail of the density at large radii 
is determined by the lightest hadron, the pion. 

Since the size o f the bump/dip signal found by 
IFriedrich and Walcherl d2003l) (FW) is in conflict wi t h cal- 
culations using dispersion relations ( Belushki n et al .1 120071 
MeiBner, 2007) (BHM), it is worthwhile to discuss the dif- 
ferences more closely. 

• The fit of FW describes the data in the space-like re- 
gion with a very good x 2 /dof w 0.9 for about 160 de- 
grees of freedom (dof), because the fitting function is 
designed for the space-like data. The dispersion rela- 
tions try to reproduce both the space-like and the time- 
like form factors by use of all the available spectral data 
for the involved hadrons, and therefore the fits are much 
more constrained. As a result, the fits of BHM have 
the much larger x 2 /dof ps 1.8 with dof ps 200. Such 
large values of x 2 have to be taken with great caution 
since they are somewhat outside the range of the valid- 
ity of statistics. In particular, the statistical probability 
P(x 2 /dof > 1.8, dof « 200) is smaller than 10~ 10 . 
This leaves the usual suspects: the probl em is in the data 
dBelushkin et all 120071: iMeiBneii 120071) . the dispersion 
relation have still an incomplete input, or both data and 
theory have problems. In any case, the l- a bands of 
iBelushkin etaH d2007l) and lMeifiner! d2007l) derived by 
increasing the absolute x 2 by 1 are not meaningful if 
the % 2 /dof is as much off as 1.8. 

• In order to obtain the bump/dip structure at Q 2 w 
0.2 GeV 2 , BHM would have to include two more "ef- 
fective" poles: an additional isoscalar pole near the 
(isoscalar) u> meson, but with the opposite sign and 
twice the strength of the u>, and a weaker isovector 
structure close to the mass of 3 pions, which is the 
threshold of the isoscalar channel. With these modifi- 
cations, also BHM could obtain a x 2 /dof w 0.9. There 
is however no evidence for such structures in e + — e~ 
collisions nor are such objects known to interact with 
the nucleon, and therefore BHM discard these fits. 

• The electric rms radius of the proton, r p E , is another 
piece of evidence showing some peculiarity around 
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0.2 GeV 2 . From a fit to all the available low Q 2 data, 
ISickl d2003l) finds the radius r p E = (0.95 ± 0.018) fm . 
On the other hand, FW obtain r E = 0.794 fm without 
and r p E — 0.858 fm w i th the bump/dip structure. Ac- 
cording to Rosenfeld er (2000), Coulomb and recoil cor- 
rections have to be added to these results, which leads 
to r p E = 0.876 ± 0.015 f m in accord with re sults from 
Lamb shift measurements dUdem et a/.Ul997l) . (For an 
overview of the results from atomic ph ysics and their 
interpr etati on, we refer to the work of iKarshen boim 
(1998) and lCarlson and Vanderhaeghenld2007l) .) How- 
ever, as in all previous work based on dispersion rela- 
tions, the electric rms radius of the proton also turns out 
to be small in the work of BHM, r p E = 0.844 fm or even 
smaller. 

• The dip structure reported by FW for the proton corre- 
sponds to a bump structure obtained for the neutron at a 
similar value of Q 2 . This change of sign makes sense, 
because the pion clou d couples to the isovec tor photon. 
As mentioned before, Kopec kv et al\ (Q997) obtained a 
mean square radius (r 2 ) E = -(0.115 ± 0.004) fm 2 
from low-energy neutron scattering off 208 Pb, however 
this extraction is certainly model depen dent. BHM get 
(r 2 )% = —0.118 fm 2 , in agreement with lKopeckv et al\ 
dl997b . FW take (r 2 )| = -0.115fm as a fixed param- 



eter or obtain 



„2\n 



-0.147fm in the new fit of the 



analytical form of the phenomenological model. 

• It follows from dispersion relations that the tail of the 
charge distributions at large radii has a Yukawa shape 
with the mass of the lightest intermediate state, that is 
the 2 pion masses for the isovector and 3 pion masses 
for the isoscalar densities. Hence it would take a consid- 
erable cancellation of positive and negative structures 
in the lower part of the spectral function if one wants to 
shift the pion cloud to rms radii above 1.7 fm. This is in 
conflict with the bump/dip structure of Eq. ( [391 , which 
results in a considerable amount of charge in the "pion 
cloud" above 1.7 fm, as seen in Fig. [6] 

• The fit of FW is restricted to the space-like form fac- 
tors. This approach can not be extended to the time-like 
region, and another purely empirical fit would make lit- 
tle sense in view of the restricted data basis in the this 
region. The dispersion relations, on the other hand, are 
built just to make the connection between the two re- 
gions, and the results of BHM give a good overall de- 
scription in both domains. However, they miss a struc- 
ture at £ beam « 1.05 GeV or t » 4.4 GeV 2 . 

In concluding these arguments, we mention that dispersion 
theory and FW agree on the dip seen for the magnetic form 
factors of both proton and neutron at Q 2 » 0.2 GeV 2 . There 
is also qualitative agreement that the charge and the magneti- 
zation in the surface region of the nucleon, r ^ 1 fm, are dom- 
inated by the pion cloud, which reaches much beyond the rms 
radius of the proton. It remains a challenge for both experi- 
ment and theory to answer the raised questions concerning the 



distributions of charge and magnetization inside an nucleon, 
which we consider a key aspect of the nucleon structure. 



D. Weak form factors of the nucleon 

1 . Axial form factor of the nucleon 

The axial current of the nucleon can be studied by 
anti-neutrino and neutrino scatterin g, pion electroproduc - 
tion, and radiative muon cap ture, see iBernard et al.\ d2002al) : 
iGorringe and Fearing! d2004l) for recent reviews. The (isovec- 
tor) axial current between nucleon states takes the form 



A, 



U P2 1^Ga{Q ) + i 



(P2 

2M 



G P (Q 2 ) 75^ 



(41) 

As for the vector current, Eq. (fTTT i, there appear two form fac- 
tors, the axial form factor Ga and the induced pseudoscalar 
form factor Gp. A linear combination of the form fac- 
tors Ga and Gp is related to the pion-nucleon form fac- 
tor G-xN by the PCAC relation. The experimental infor- 
mation about the induced pseudoscalar form factor is lim- 
ited. The data are mostly obtained from muon capture by 
the proton, fi~ + p — > n + v^. This determines the value 
of G P at Q 2 = 0.88 m 2 » 0.01 GeV 2 , which is usually de- 
scribed by the induced pseudoscalar coupling constant, gp — 
jj(jGp(Q 2 = 0.88m 2 ). A recent experiment at P SI yielded 
the value g P = 7.3 ± 1.1 (Andr eev et at „ 2007), in agree- 
ment w ith the result from heavy baryon ChPT dBernard et all 
l2002al) . gp = 8.26 ± 0.16, a nd manifestly Lorentz-invariant 
ChPT dSchindleref a/.Ll2007h . gp = 8.29 ± 0.7, with an esti- 
mated error stemming mostly from the truncation of the chiral 
expansion. The axial form factor Ga is usually parameterized 
in the dipole form of Eq. d37| i, with a parameter called the 
"axial mass", 



G A {Q 2 



9A 



(1 + Q 2 /Ai) 2 



(42) 



with g A = 1.2695(29) dYao et a/1120061) . A recent (corrected) 
global average of the axial mass as determin ed by neutrino 
scattering has been given by Bud d et al\ d2003l) . 



A", 



(1.001 ± 0.020) GeV. 



(43) 



However, quite a different value, A A = (1.20 ± 0.12) GeV 
has been derived by the K2K C ollaboration from q uasi-elastic 
v^n — > n~p in oxygen nuclei dGran et al. , 2006). The axial 
form facto r has a lso been studied by pion electroproduction 
(Baumann, 2005). The Rosenbluth separation of these data 
is shown in Fig. |T0j The results are i n agre ement with an 
earlier experiment by iLiesenfeld etal (1999), but are com- 
plemented by a data point at the very low momentum transfer 
of Q 2 — 0.058 GeV 2 . An exact Rosenbluth separation is 
prerequisite, because the transverse cross section ot is sen- 
sitive to Ga and the longitudinal cross section to the pion 
form factor F n , which is discussed in subsection IIII.EI These 
electroproduction data have been analyzed with MAID2007 
as follows: (I) The cross sections <jt and er^ of MAID were 
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FIG. 10 The transverse and longitudinal cross sections f or charged 
pion p roduction as function of Q 2 . The data are from iBaumand 
20051) . Dashed line s: cross sections predicted by MAID2007 



Drechsel et all l2007h . solid lines: fit to the data with MAID2007 
(See text for details.) Figure by courtesy of L. Tiator. 



normalized to the data by factors 0.825 and 0.809, respec- 
tively, and (II) the axial dipole mass and the corresponding 
"mass" for the monopole form of the pion form factor were fit- 
ted. The result is A^ = (1.028 ± 0.025) GeV. However, this 
value has to be corrected for the "axial mass discrepancy", 
A^ — ss 0.055 MeV, which is due to loop corrections 
dBernard et all 19921). With th is correction, the electropro- 
duction data o f lBaumannld2005l) yield 



a coir 
A A 



(0.973 ± 0.025) GeV, 



(44) 



which agrees with the corrected value from neutrino scattering 
gi ven by Eq. d43l), but dis agrees with both the pre vious result 
of Liesenfeld et al. I (fl999h and the measurement of lGran et al\ 
(2006). In view of the relatively large normalization factor 
applied to the electro-production data, it would be helpful to 
check the normalization at Q 2 = by pion photoproduction. 



2. Strangeness content of the nucleon 



As outlined in subsection III. CI the parity violating com- 
ponent of electron scattering provides access to the weak 
form factors Ge and Gm- These form factors are related to 
the strangeness content of the nucleon by the universality of 
the electroweak interaction with the quarks. For a detailed 
derivation of the strange form factors and their experimen- 
tal det ermination see, e.g., the review of iBeck and McKeownl 
(120011) . Because the strangeness in the nucleon appears only 
through the presence of the heavy ss pairs, these observables 
are of great importance for our understanding of the nucleon 
in terms of large vs. small scales. The strangeness content 
is related to the a term, which has been derived from pion- 
nucleon scattering at the (unphysical) Cheng-Dashen point , 
s = u = M 2 , t = 2ml dSainiol [20021: iThomas and Weisel 
120011) . This term is a direct measure of the chiral symmetry 
breaking in QCD, the chiral properties of the strong interac- 
tions, and the impact of sea quarks on the nucleon's properties. 



Its relation to the strangeness contribution is given by 

(N\m(uu + dd- 2ss)\N) 
a = 1 ' 



(45) 



where m = (m u + m<j) /2 is the average of the u and d quark 
masses, and y is a measure for the scalar strange quark content 
of the nucleon, 



2(N\ss\N) 
(N\uu + dd\N) 



(46) 



From a recent detailed analysis, [Pavan et all d2002l) found the 
value y w 0.46, indeed a surprisingly large strangeness con- 
tent in the nucleon, w hereas a muc h smaller value was ob- 
tained in earlier work dSainioll2002h . These inconsistencies 
were a very strong motive to study the strangeness content 
with the electromagnetic probe. At large momentum trans- 
fer, 1 GeV 2 < Q 2 < 100 GeV 2 , the strangeness contribution 
has been derived from unpolarized deep-inelastic lepto n scat- 
tering at the Fermi Lab Tevatron (Baza rko et al.l\l99$) . The 
momentum fraction of the sea quarks carried by the strange 
quarks extracted is 



(x(s(x) + s(x))) 
(x(u(x) + d(x))) 



0.5, 



(47) 



or about 3% of the total nucleon momentum. If this con- 
tribution is extrapolated to large spatial scales by the quark 
evolution, a rather small value is obtained. 

On the theoretical side a plethora of nucleon models have 
usually predicted strangeness for m factors of consi derable 
size, see for example the review of IBeck and Holsteinl dlOOl). 
These considerations have initiated an intense experimen- 
tal program at several laboratories. These activities started 
at the Bates/MIT laboratory with the SAMPLE experiment, 
which first proved that it is feasible to measure the small 
asymm etries of order 10~ 6 in parity-violatin g electron scat- 
tering dKowalskil 120061 ISpavde et all |2004 . This experi- 
ment was based on a particular technique using Cherenkov 
detectors developed previ ously for a parity -violation experi- 
ment at the Mainz linac (Hei l et all 1 19890 an d an improve- 
ment of the SLAC polarized electron source (ISouder et all 
1990). At the Mainz Microtron MAMI, the A4 collabora- 
tion built a Cherenkov detector consisting of 1022 PbF2 crys- 
tals, which in conjunction with electronics allowing for on- 
line identification of electr omagnetic clu sters, made it possi- 
ble to count single events (Maas, 2006). Furthermore, two 
experiments were performed at the Jefferson Lab. The first 
of these experiments (HAPPEX) used the two-spectrometer 
set-up of Hall A taking advantage of a pair of septum mag- 
nets for measurements at very small scattering angles and 
low momentum transfers. This project was passing through 
different phases of improvement. HA PPEX-I measured on 
a hydrogen target dAniol et all |2004 at Q 2 = 0.48 GeV 2 
only. In this geometry the combination G S E + 0.392 G S M = 
0.014 ± 0.020 ± 0.010 was determine d. In the next step 
HAPPEX-II measured on both hydrogen dAniol et aZ.ll2006d) 




FIG. 11 The strangeness form factors G% and G S M at Q = 
0.1 GeV 2 as obtained by the SAMPLE, A4, HAPPEX, and GO exper- 
iments. The bands represent the one-sigma (statistical plus systemat- 
ical) error of the individual experiments, the ellipse is the combined 
two-si gma area for all measurements. The figure is from lKowalskil 
(2006). 



FIG. 12 The one-sigma band of HAP PEX-II toge ther with the one- 
and two-sigma ellipses derived by Acha et al ( 2007). The theoretical 
predictions with their e stimate errors are identified by the numbers 
given in the figure: [181 Park and Weiael ( 1992), [19] Hammer et al 
1996), [20] Hammer and Ramsey-Musolf (1999 Q, [21]"lSilva et al 
f 



20021) . [22] iLewis et al\ d2003l) . r23l lLeinweber et al\ d2005l.l2006l) . 



and helium targets dAniol et all 12006 b). The nucleus 4 He is 
quite special as a target, because only the electric form fac- 
tor can contribute due to its zero total spin. The results of 
HAPPEX are compared to those of other collaborations in 
Fig-El Each measurement gives an error band in the plot of 
G% versus G S M . The common error ellipse indicates values for 
Gs(0.1 GeV 2 ) and Gm(0.1 GeV 2 ) that are consistent with 
zero but at variance with most theoretical predictions, how- 
ever, not incompatible with the experiments obtained with the 
other methods mentioned above. Recently the second phase 
of HAPPEX-II w as completed, with the result of a much im- 
proved precision dAcha et al. 1 12007b . These results are com- 
pared with several theoretical predictions in Fig. [12] As is 
obvious from the figure, the strangeness form factors are cen- 
tered about zero, whereas most of the models predict large 
values. The only theoretical results compatible with these ex- 
periments are from lattice gauge calc ulations with chiral ex- 
trapo l ation to the physica l pion mass dLeinweber et al.l 120051 
120061: ILewis etal\. 120031) . The second JLab experiment was 
performed by the GO collaboration. This collaboration has 
built an e ight-sector superconduc ting toroidal magnetic spec- 
trometer dArmstrong et al.l 120051) . Figure [T3] displays the Q 2 
dependence of the world data including the GO results. From 
this figure we get the impression of a small but finite value for 
that particular combination of the two strangeness form fac- 
tors. 
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FIG. 13 The form factor combination G 3 E (Q 2 ) + vGliiQ 2 ) as ob- 
tained by kinematical extrapolation from A4, HAPPEX II, and GO 
experiments. The point "HAPPEX III" i ndicates the erro r bar for a 
planned measurement. The figure is from[Kpwalski (20061). 



E. Form factor of mesons 

Because mesons are unstable particles, their form factors 
can not be measured directly by lepton scattering but have 
to be obtained by more indirect methods. In the following 
we concentrate on the form factor of the charged pion, how- 
ever similar methods can also be used t o measure the form 
factors of heavier me s ons or rare decays dGuidal et al. 1 119971: 
IVanderhaeghen eiall 1 1998b . At present, only the data for 
the pion are precise enough to allow for a reliable extraction 
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of the form factors over a large region of momentum trans- 
fer. There exist two experimental methods to overcome the 
missing target problem. The first method is the scattering 
of relativistic mesons with dilated lifetime on atomic elec- 
trons, which are then identified by measuring the recoil of the 
struck electron. This method is limited to relatively small mo- 
mentum transfer, Q 2 < 0.5 GeV 2 . As a consequence, this 
method is essentially sensitive to the rms radius of the free 
pion, r n = -J (r 2 )^, which is related to the mass parameter 
A 2 in the usual monopole form as follows: 



i + Q 2 /H 



O (Q 4 ) • (48) 



From an exp eriment at the CERN SPS, Ame ndolia et al\ 
(11984 [l986b) derived the rms charge radius of the pion, 
TV = (0.663 ± 0.006) fm. At the same time also th e kaon 
form factor was measured (lAmendolia et all [1986a). The 
rms charge radius of the charged kaon was found to be 
fx = (0.58 ± 0.05) fm, somewhat smaller than the pion 
radius, which is to be expected because of the heavier strange 
quark in the kaon. Because of the small momentum transfer 
involved, these results depend only little on the monopole 
form of the ansatz. 

The second possibility to study the pion form factor is given 
by electron scattering on the pion cloud of the nucleon, which 
is part of the reaction p(e, e'Tr + )n. The obvious problem is 
that the initially bound pion is off its mass-shell and that many 
other diagrams contribute as well. The idea is therefore to 
study this process in kinematic regions for which the t-channel 
pion exchange is dominant. In principle one should extrapo- 
late the cross section to the pion pole, which however lies at 



the unphysical 4-momentum transfer t 



At very high 



momentum transfer, Q 2 ^> 1 GeV 2 , this extrapolation can 
be pe rformed within the Regge model (jVanderhaeghen et all 
1998). Following this approach. iHorn et al.\ (120061) have re- 
cently determined by a Rosenbluth separation of the lon- 
gitudinal and transverse cross sections at the Jefferson Lab. In 
this way they extracted two precise values of the pion form 
factor at Q 2 = 1.60 GeV 2 and 2.45 GeV 2 , which are shown 
together with previous data and several model calculations in 
Fig- 03 Beside the model dependence on the Regge analysis, 
another problem arises because the data at the higher momen- 
tum transfer do not follow the monopole form, i.e., the value 
of A 2 diffe rs by as much as 7% fr om A 2 = 0.53 GeV 2 as 
obtained by Amendolia et aZ.I(ll986bl) . This leads to an incon- 
sistency of about 1 standard deviation. At small values of Q 2 
one can also try to derive the rms radius from pion electro- 
production as pointed out in subsection IIII.D.ll A previous 
result yielded A 2 = 0. 425 GeV 2 , equivalent to an rms radius 
7v = fO.74 ± 0.03) fm dLiesenfeld et a/.lll999l). This result is 
clearly at variance with the value o flAmendolia et al. I(fl986b). 
A reason f or this discrepancy was given by iBernard et al.\ 
d2002all2000h in terms of the loop corrections for the longitu- 
dinal s-wave multipole Lq + \ which dominates the cross sec- 
tion <jl at small Q 2 . In fact, these loop corrections increase 



the downward slope of Lq + ' substantially, such that 

(r 2 > 7r [p(e, e' 7r+)n] = (r 2 )^ + 0.266 fm 2 . 



(49) 



With t he pion radius according to Amendoli a et al. I (11984 



1986b), 



0.440 fm , the electroproduction experiment 



should therefore measure (r 2 ) w [p(e, e' 7r + )n] = 0.706 fm 2 
or an "effective" rms radius of 0.84 fm, which is surprisingly 
cl ose to the nucleo n radius. The fit with MAID07 to the data 
of lBaumannl(l2005l) yields A 2 = (0.386±0.042) GeV 2 , which 
is smaller than the result of iLiesenfeld et al. (1999). If we 
include the loop correction, we obtain rv = (0.78 ± 0.04) fm, 
which is quite close to the prediction of ChPT Of course, 
there are weak points in our reasoning. In the first place, there 
is the already discussed overall reduction of the MAID model 
to fit the data. Second, the s-wave yields only about half of the 
measured cross section, and much smaller loop corrections 
are expected for the higher partial waves. Nevertheless we 
may conclude that the virtual constituent pion looks quite 
different from the free pion and appears, in this particular 
experiment, nearly as large as the nucleon. 

Because the neutral pion is its own antiparticle, its form fac- 
tor vanishes identically. However, the reaction 7* (Q 2 ) + 7 — > 
7T° can be studied as function of Q 2 . This provides informa- 
tion on the transition form factor .F~~~_o of the Wess-Zumino- 



Witten anomaly defined by Eq. (O of section U As shown 
in section [IV] this anomaly is quite important in Compton 
scattering and in particular for the spin polarizability of the 



nucleon. The transition form factor F, 



7*77r L 



was measured 

at the e + e~ co llider at Cornell by the CLEO collaboration 
dGronberg et aZ.lll998l) and at the PETRA storage ring by the 
CELLO collaboration dBehrend et adll99l1) . Analyzed with a 
monopole form factor as given by Eq. d48l . these experiments 
yielded the parameter A 7 * 77r a = (776 ± 10 ± 12 ± 16) MeV, 
and similarly for the corresponding transition form factor of 
the 77 meson A T , 77) = (774 ± 11 ± 16 ± 22) MeV. Both results 
confirm the prediction of the vector dominance model assum- 
ing that the virtual photon coexists with the neutral p meson, 
which decays in a pion or an 77 and a real photon. The Q 2 de- 
pendence of the transition is then simply given by the propaga- 
tor of the p, that is, A « m p = 775.5 MeV for both reactions. 
With the same arguments as for the ordinary form factors, we 
can turn this value in a transition radius ryvzw = 0.62 fm. 



IV. POLARIZABILITIES 

The polarizability measures the response of a system to a 
quasi-static electromagnetic field. In particular the energy of 
a homogeneous and isotropic system is described by the elec- 
tric (ctEi) and magnetic ((3mi) dipole polarizabilities. In the 
case of a macroscopic system these polarizabilities are related 
to the dielectric constant e and the magnetic permeability p. 
The classical theory of Lorentz describes the dispersion in 
a medium in terms of electrons bound by a harmonic force. 
In the presence of a static and uniform electric field Eq, the 
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FIG. 14 The pi o n form factor F n multiplied with Q 2 . Crosses: 
lAmendolia eIo/1 dl984l. Il986bl) (it — e scattering at CERN), as- 
terisks: I Ackermann et al\ d 19781) (DESY), diamonds: iBrauel et al\ 
dl977l 1 19791) (DES Y), circles: i Tadeyosvan et al\ d2007h (Jefferson 
Lab), and squares: iHorn et al\ d2006f) (Jefferson Lab). The experi- 
mental results are comp ared to four model calculations as indicated 
by the lines. Figure from lHorn et al. (2006). 



Hamiltonian of the harmonic oscillator takes the form 



H 



p mio, 
2~^ + ~ 



+ e r ■ E 



(50) 



We note that Gaussian units are used in this subsection in or- 
der to concur with the standard notation of classical electrody- 
namics, i.e., the fine structure constant is related to the charge 
by a cm = e 2 ~ 1/137. In all other parts of this review we 
use Heaviside-Lorentz units, i.e., a cm = e 2 /(4ir) ps 1/137. 
Substituting r = r' + Ar, where Ar = — eEo / (mui 2 ) is the 
displacement due to the electric field, we may rewrite Eq. 
as 



H 



Ira 



AE , 



(51) 



that is, the applied electric field has induced a dipole moment 
d and an energy shift AE, 



d = -eAr = T En 



AE 



' 2niLol 



(52) 



The electric dipole polarizability olei is obtained by varying 
the induced dipole moment or the energy shift with regard to 
the electric field, 



Sd 
SE 



5 2 AE a 



(SE y 



(53) 



This result is quite general and also valid for quantum me- 
chanical systems. The energy shift to first order (linear 



Stark effect) vanishes for a system with good parity, and if 
the ground state is spherically symmetric, the second order 
(quadratic Stark effect) yields 



AE 



-at c 



n>0 



\(n\z\0)\ 2 
- eo 



(54) 



with the electric field pointing along the z-axis, and e n the 
energies of the eigenstates \n >. Combining Eqs. d53l and 
|, we obtain the (quasi-static) electric dipole polarizability, 



OLEl 



2a cm 



n>0 



\(n\z\0)\ 2 



(55) 



As an example for a classical extended object we quote the po- 
larizabiliti es of a small dielectric and permeable sphere with 
radius a (|JacksonLfl975l) . 



OLEl 



e + 2 



Ml 



1 



/' 



(56) 



The polarizabilities for a perfectly conducting sphere are 
obtained from Eq. (TS6b in the limits e — > oo and /j, — > 0, 
otEi = a 3 and Pmi = — ^ a 3 - Up to a factor 47r/3, the 
electric polarizability of a conducting sphere is the volume 
of the sphere. Because of the different boundary conditions 
for the magnetic field, the magnetic polarizability turns out 
negative. The induced currents in the conductor lead to 
a magnetization opposite to the applied field according to 
Lenz's law, i.e., diamagnetism. A permeable sphere can 
be diamagnetic (fi < 1) or paramagnetic (/i > 1), in the 
latter case the magnetic moments are already preformed and 
become aligned in the presence of the external field. Whereas 
the magnetic polarizabilities of atoms and molecules are usu- 
ally very small, \fi — 1| < 10 -2 , the electric polarizabilities 
may be quite large compared to the volume. For example, 
with a static dielectric constant of s = 81, water is a nearly 
perfect conductor, although in the visible range this constant 
is down to e — 1.8, corresponding to a refraction index 
n = 1.34. A further, quantum mechanical example is the 
hydrogen atom. Its ground state has good parity and spherical 
symmetry and therefore Eq. d55l ) applies. It is even possible 
to sum over the excited states and to obtain the closed form 
olei C^H) — I a^, where as is the Bohr radius dMerzbacheil 
Il970h . With an rms radius given by (r 2 ) = 3a 2 B , the equiva- 
lent hard sphere has the radius R = V5as, and as a result 
the hydrogen atom is a pretty good conductor, a^i/volume 
re 1/10. 

In the following we report on the polarizabilities of the nu- 
cleon (subsections lIV.Al and lIV.BI ) and pion (subsection llV.CI ). 
Both particles are very rigid objects. They are held together 
by strong interactions, and the applied electromagnetic field 
can not easily deform the charge distribution. If compared 
to macroscopic matter, the nucleon is a dielectric medium 
with e ps 1.002, that is a very good insulator. Furthermore, 
magnetic effects are a priori of the same order as the elec- 
tric ones, because the charged constituents, the quarks and 



mesons, move close to the speed of light. However, the dia- 
magnetic effects of the pion cloud and the paramagnetic ef- 
fects of the quark core of the nucleon tend to cancel, with the 
result of a relatively small net value of (3mi- The polariz- 
ability of the nucleon can be measured by Compton scatter- 
ing: The incoming photon deforms the nucleon, and by mea- 
suring the energy and angular distributions of the outgoing 
photon one can determine the induced current and magnetiza- 
tion densities. Particularly interesting is the case of "virtual 
Compton scattering" (VCS), which yields information on the 
spatial distribution of the polarization densities. Furthermore, 
the nucleon has a spin and therefore polarized nucleons ap- 
pear as anisotropic objects. This leads to the spin or vector 
polarizabilities whose closest parallel in classical physics is 
the Faraday effect. 



A. Real Compton scattering 

The reaction 7(3, e) + N(p,X) -> 7(9', e') + N(p',X') 
involves the absorption of an incident real photon with 
4-momentum q and polarization e on a nucleon with 4- 
momentum p and polarization A, leading to a spectrum of in- 
termediate hadronic states, which finally decay by the emis- 
sion of a real photon leaving the nucleon back in its ground 
state. Typical intermediate states are shown diagrammatically 
in Fig. [15] and the following Fig. [16] shows the contribu- 
tions of these diagrams to the differential cross section. For 
a point Dirac particle only the diagrams (a) and (b), with a 
nucleon in the intermediate state would contribute. These two 
"nucleon Born terms" yield singularities for the (unphysical) 
kinematics s = M 2 and u = M 2 , respectively. The differen- 
tial _cross_S£cdonforsuch a point nucleon was first calculated 
by iKlein and Nishrnal d 19291) . The predicted cross section in- 
creases surprisingly much by adding the Pauli current, i.e., the 
anomalous magnetic moment of the nucleon. The result is the 
Powell cross section (Powell, 1949). If we further add the 
"pion pole term" of Fig.[l5](f), the cross section falls back to- 
wards the Klein-Nishina values. The pion pole term has a sin- 
gularity at t = rn 2 , it results from the decay ir° — ► 7 + 7 as a 
co nsequence of the axial ano mal y derived on ge neral grounds 
bv lWess and Zuminol (1 197 lb and I Wittenl dl983l) . This term is 
often referred to as triangle anomaly, because the vertex ^77 
can be microscopically described by a triangular quark loop, 
a diagram not allowed in any classical theory and only ap- 
pearing due to the renormalization process of quantum field 
theory. As we see from Fig. [16] the pion pole term yields a 
large contribution for backward angles. All further contribu- 
tions in Fig. [15] do not have pole structures, but correspond 
to excited states in s-, u- or t-channel processes. As such 
they yield dispersive contributions that determine the polariz- 
abilities of the nucleon. If we include only the electric and 
magnetic dipole polarizabilities, we obtain the low energy ex- 
pansion (LEX) in Fig. [16] This expansion describes the data 
only up to a photon lab energy of about 80 MeV, over a region 
in which the polarization effects are small and the data scat- 
ter. Therefore, the analysis of the modern data has been based 
on dispersion relations whose results are shown by the solid 
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(e) (f) (g) 



FIG. 15 Typical intermediate states contributing to Compton scatter- 
ing off the nucleon. Upper row: The direct (a) and crossed (b) Born 
diagrams with intermediate nucleons, a typical resonance excitation 
in the s-channel (c) and its crossed version (d). Lower row: mesonic 
contributions with photon scattering off an intermediate pion (e), the 
pion pole diagram (f), and a correlated two-pion exchange such as 
the "a meson" (g). 
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FIG. 16 Differential cross section for Compton scattering off the 
proton as function of the photon lab energy i? 7 and at scattering angle 
#kb ~ 135°. The lines show the results of fixed-f subtracted disper- 
sion relations (solid), Klein-Nishina (small dots), Powell (dashed), 
Powell plus 7T° pole (large dots), and LEX including the leading or- 
der contributions of q^i and , (dash e d-dotted). The experimen - 
tal data are fromlFederspiel et al] dl99ll) . iMacGibbon et al.\ |1995 | ) . 
and lOlmos de Leon et q/.l J200lh . Figure from lDrechsel et q/.H2003l) . 



line in the figure. Clearly the higher order terms become more 
and more important with increasing photon energy, particu- 
larly after crossing the pion threshold (seen as a kink at about 
150 MeV), from thereon the energy increases sharply towards 
the A(1232) resonance. 

1 . Compton amplitudes and polarizabilities 

Assuming invariance under parity, charge conjugation and 
time reversal symmetry, the general amplitude for real Comp- 
ton scattering (RCS) can be expressed by 6 independent struc- 
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ture functions Ai(y, t) depending on the two Lorentz invariant 
variables v = (s — u)/(AM) and t, see Eq. (0 for the defini- 
tions. These variables are related to the initial (E 7 ) and final 
(E') photon lab energies, and to the lab scattering angle 9\. d \, 
as follows: 

t = -AE 1 E! y sin 2 ^ = -2M(E^-E' J ), 

* = E i+m = V E ^ + ^ ) - (57) 

Note that for 9\. d \, = (forward scattering) the Mandelstam 
variable t vanishes and the crossing-odd variable v is simply 
the incident photon lab energy E~ . 

The general RCS scattering matrix takes the form 

Tfi = e M e 1 : u f 0', X' N )H^ Ui (p, X N ) , (58) 

where u and u are the nucleon spinors. The Compton tensor 
H^ v contains the hadronic transition currents J M and J v as 
well as the propagation of the intermediate hadronic state. It 
can be decomposed into a complete basis of 6 tensor structures 
constructed from the indep endent momentum 4-vecto rs and 
appropriate Dirac matrices dL vovl \l9S ll iPrangel 1 1 958l) . 

H»» = M>t v Ai{v,t) . (59) 

i=l,6 

For further details we refer to Drechse l et al. (2003). 

Let us now consider the forward scattering of a real pho- 
ton by a nucleon. The incident photon is characterized by the 
Lorentz vectors of momentum, q = (\q\ ,q), and transverse 
polarization, s\ — (0, E\), with q ■ q — for real photons and 
£\ ■ q = 0. The corresponding quantum numbers of the out- 
going photon are denoted by primed quantities. If the incident 
photon moves in the direction of the z-axis, q = \q \ e z , the 
polarization vectors 

e± = =F-^= (e x ± ie y ) (60) 

correspond to circularly polarized light with helicities A = 
+1 (right-handed) and A = — 1 (left-handed). The forward 
Compton amplitude takes the form 

T(v, 6 = 0) = s'* • ef(u) + i 3 ■ {t * x e ) g{v) . (61) 

Because T is invariant under the crossing transformation, 
e'* «-* e and v — ► —u, / must be even and g odd as 
function of v. These fo rward scattering amplitude s have 
the low-ene rgy expansion dGell-Mann and G oldberger, 1954; 
EowL IT954h 

= ~fi + ^ + ^i) v 2 + O(^) , (62) 
9{v) = -^^ + ^ 3 + 0(^), (63) 

with 6 ]v the charge of the nucleon in units of e and kn the 
anomalous magnetic moment in units of nuclear magnetons. 



The leading contribution to f(y) is the Thomson term famil- 
iar from non-relativistic theory. The term linear in v vanishes 
due to the crossing symmetry, and the term 0(v 2 ) contains 
the sum of the scalar polarizabilities giving information on 
the internal structure. Being odd under crossing, the spin-flip 
amplitude g(y) starts with the term 0(y) proportional to the 
square of the anomalous magnetic moment, and its next order 
term is described by the forward spin polarizability 70. The 
leading terms for both amplitudes are obtained from the pole 
terms typical for a point-like particle, whereas the polarizabil- 
ities are contained in the sub-leading terms. As is evident from 
the above equations, the scalar and spin polarizabilities have 
different units. In the following all scalar polarizabilities are 
given in units of 10~ 4 fm 3 , while the vector or spin polariz- 
abilites have units of 10~ 4 fm 4 . As will be detailed in sec- 
tion |VTA] the forward scalar (olex + [3 mi) an d forward spin 
(70) polarizabilities of Eqs. d62b and d63l can be determined 
by energy-weighted integrals over the photoabsorption cross 
sections. In particular, Baldin' s sum rule yields the f ollowing 
results for proton and neutron (Ba busci et q/|[l9~98bl) : 

(«bi+/3mi) p = 13.69 ±0.14, 
(asi+/3 M i)n = 14.40 ±0.66. (64) 

The T-matrix for general scattering angles is described by 
the 6 Lvov amplitudes Aj(z/, t). These amplitudes have no 
kinematical constraints, are symmetrical under crossing, and 
contain both the pole terms of (Fig.[15]a, b, and f) and an inte- 
gral over the excitation spectrum, which we call the dispersive 
amplitude, 

Ai(y,t) = Af e (v,t) + Af sp (v,t). (65) 

The polarizabilities are determined by the dispersive ampli- 
tudes at v = t = 0, that is, at the threshold for RCS. This 
defines 6 real numbers a; = Af sp (0,0), from which we can 
derive 2 scalar and 4 vector (or spin) polarizabilities by lin- 
ear combinations. In the scalar sector, we find the familiar 
electric (a^i) and magnetic (/3jwi) polarizabilities, which ap- 
pear as qei + (3m 1 for forward and a,E\ — Pmi f° r backward 
Compton scattering. The physical content of the 4 vector po- 
larizabilities is best described in a multipole notation. Since 
the initial and the final states contain a nucleon in its ground 
state with total spin J = |, the transition operator must have 
even parity and angular momentum or 1 . The electric polar- 
izability describes the absorption of an electric dipole photon 
followed by the emission of a photon with the same multipo- 
larity, that is ole\ ~ [E\ x Ei] ^ , and in the same way we find 
Pmi ~ [Mi x Mi][°J with the multipoles coupled to (scalar 
polarizabilities). In the spin-dependent sector there are 4 po- 
larizabilities at the lowest order: jeiei, Imimi, Imiei, and 
1e\mi- In this case the multipolarities are coupled to 1 (vec- 
tor polarizabilities). As an example, Jmie2 ~ [Mi x E^]^ 
defines a spin polarizability with an electric quadrupole ab- 
sorption followed by a magnetic dipole emission. It is useful 
to define the forward (9 = 0) and backward (9 — 71") polariz- 
abilities also in the spin-dependent sector, 

70 = — 7-B1E1 - 7A/1A/1 — 1MIE2 ~ lEXMI , (66) 
7?r = — 7-E1E1 + 1MIMI + 1MIE2 ~ 1EIM2 ■ (67) 



23 



It is of course possible to define higher polarizabilities related 
to the v and t derivatives of Af sp (v, t) taken at v = t = 0. An 
often discussed example is the electric quadrupole polarizabil- 
ity, which appears among the terms of 0(v A ) in Eq. 



2. Theoretical developments 

In a nonrelativistic approach like the constituent quark 
model (CQM), the scalar dipole polarizabilities can be ex- 
pressed by 

a ei = 2a em ~ ™ 1" & a Ei , (68) 



Pmi = 2a em 2J 



E n — Eq 



A/3 M i , (69) 



where d = d q — e q^q an( ^ P = 2 fiq are sums over the 
electric and magnetic dipole operators of the constituents. For 
simplicity the quark masses may be taken as m q — \M, and 
the quark charges e q are in units of e. Clearly the first terms on 
the rhs of the above equations are positive, because the exci- 
tation energy E n — Eq is positive. The second terms describe 
recoil and retardation, Auei = a em (0| Yl e q ?\ |0)/(3M) 
and A0 M1 = -a em (0| d 2 + £ &\ |0}/(2M). These are 
small corrections in atomic physics but quite sizeable for the 
quark dynamics of the nucleon. They turn out positive for 
a ei but negative for Pmi- The leading term of the magnetic 
polarizability describes the paramagnetism, mainly by a quark 
spin-flip transition from the nucleon to the A (1232), while the 
sub-leading term represents Langevin's diamagnetism. The 
simple CQM with an oscillator potential connects the rms ra- 
dius (r 2 ) 1 / 2 with the oscillator frequency, mo = 3/(M(r 2 )), 
and yields oei = 2a em /(Mcc> 2 ) + 0(M~ 2 ). However, this 
model is not able to describe both size and excitation energy. 
If we use the proper size, cxei is grossly overestimated, 
whereas a fit to the excitation energy of the dominant dipole 
mode N* (1520) leads to a value much below the experiment. 
For the magnetic polarizability, the Ml transition to the 
A (1232) yields a large paramagnetic value, /3 M1 w 12, 
which is somewhat reduced by the subleading diamagnetic 
terms. It was therefore early recognized that a complete 
picture of the nucleon must also include the pion cloud 
dWeiner and Weiselll985l) . 

Systematic calculations of pion cloud effects became pos- 
sible with the development of chiral perturbation theory 
(ChPT), an expansion in the external momenta and the pion 
or quark mass (p expansion). The first calcula tion of Comp- 
ton sca ttering in that scheme was performed bv lBernard et al\ 
( 1991b). Keeping only the leading term in l/m n , they found 
the following remarkable relation at O (p 3 ) : 



&E1 



10/3 



Ml 



12.2. 



(70) 



with f n w 93 MeV the pion decay constant and gA ~ 1-26 
the axial coupling constant. The calculation was later re- 



peated in heavy-baryon ChPT, which allows for a consistent 
chiral power counting, and extended to Q(p A ) yie lding 
o&i = 10.5 ± 2.0 and P M1 = 3.5 ± 3.6 dBernard et all 
1994, 1993). The error bars for these values indicate that 
several low-energy constants appear at this order, which 
were determined by resonance saturation, that is by use 
of phenomenological information about resonances and 
vector mesons. Since the A(1232) is close in energy and 
very important for photoabsorption, it has been proposed 
to include this resonance dynamically. This leads to an 
additional expansion parameter, the NA mass splitting (e 
expansion). Unfortunately, the "dynamical" A increases the 
polarizabiliti es to values far above the data, a p E1 = 16.4 and 
(3 P M1 = 9.1 dHemmertef a/.LIl998l) . This can be changed by 
introducing large low-energy constants within a higher-order 
calculation, however at the expense of losing the predictive 
power. 

The spin polarizabilities have been calculat ed to 0(p 3 ) 
in both relativis tic ChPT dBernard et all 1 19951) and heavy- 
baryon ChPT dHemmert et all 1998 ). As an example we give 
the predictions of the latter reference: 
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4.6 
4.6 



2.4 
2.4 



0.2 
0.2 



= 2.0, 
43.5 = -36.7. 



(71) 
(72) 



1 



the 4 separate contributions referring to N7r-loops, A-poles, 
A7r-loops, and the pion pole, in order. As is obvious from 
these results, the tt° pole dominates the backward spin 
polarizability but does not contribute in the forward direction 
Independent calculations of the forward spin polarizabilit 
to 0(p 4 ) resulted in 70 = -3.9 dBirse et all 1200 it Weta 
2000), which indicates a slow convergence of the expansion. 

Because a reliable data analysis is based on dispersion re- 
lations (DRs), we recall some pertinent features of this tech- 
nique in the following. The invariant amplitudes Ai are free 
of kinematical singularities and constraints, they also obey the 
crossing symmetry and gauge invariance. Assuming further 
analyticity and an appropriate high-energy behavior, these 
amplitudes fulfill unsubtracted DRs at fixed t, 



Re Ai{v, t) = Af e (v, t) + -V I dv 



, v'lm Ai(i/,t) 



(73) 

where A v ° le is the nucleon pole term and V denotes the 
principal value integral. The latter can be calculated if 
the absorptive part of the amplitude, Im Ai, is known to a 
sufficient accuracy. Because of the energy-weighting, the 
pion production near threshold and the mesonic decay of the 
low-lying resonances yield the biggest contributions to the 
integral. With the existing information on these processes 
and some reasonable assumptions on the lesser known higher 
part of the spectrum, the integrand can be constructed up to 
cm energies W « 2 GeV. However, Regge theory predicts 
that the amplitudes A\ and A2 do not drop sufficiently fast to 
warrant a convergence of the integral. This behavior is mainly 
due to fixed poles in the i-channel. In particular the f-channel 
exchange of pions and a mesons leads to the bad convergence 
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for A-2 and Ay, respectively. The latter meson has a mass of 
about 600 MeV and a very large width, it models correlations 
in the two-pion channel with spin zero and positive parity. 
In order to obtain useful results for these two amplitudes, 
iL'vov et all d 19971) proposed to close the contour integral in 
the complex plane by a semi-circle of finite radius v ma x, and 
to replace the contribution from the semi-circle by a number 
of energy independent poles in the t channel. This procedure 
is relatively safe for A^ because the ir° pole or triangle 
anomaly is well established by both experiment and theory. 
However, it introduces a considerable model-dependence for 
A x . 

In order to avoid the convergence problem and the phe- 
nomenology necessary to determine the asymptotic contri- 
butions, it was sugge sted to subtract the DRs at v = 
dDrechsel et all |2000|) . This subtraction improves the con- 
vergence by two additional powers of v' in the denominator 
of the dispersion integrals, Eq. d73l >. The subtraction func- 
tions Aiiy = 0, t) can be obtained from subtracted DRs in 
t with the imaginary part of the amplitude 77 — * tttt — > 
NN as input. In a first step, a unitarized amplitude for the 
77 — > tttt subprocess is constructed from the available ex- 
perimental data. This information is then combined with the 
tttt — > NN amplitudes determin ed by analytic al continua- 
tion of nN scattering amplitudes ( IHohleii U983). Once the 
t dependence of the subtraction functions Ai(0, t) is known, 
the subtraction constants <ij = Ai(0, 0) have to be fixed. Al- 
though all 6 subtraction constants a\ to a$ could be used 
as fit parameters, it is sufficient to fit a\ and 122, or equiv- 
alently (ccei — Pmi) and 7^ to the data. The remaining 
4 subtraction constants can be calculated through an unsub- 
tracted dispersion integral. Yet another method are hyperbolic 
(fixed-angle) DRs, which improve the conve rgence for large 
value s of t or backward s c attering angles ( Bernabe u et all 
19741 iHolstein and Nathan) Il994t t7vov and Nathanl 1 1 999b . 



Holstein an d Nathan ( 1994) investigated backward DRs in or- 
der to get rigorous bounds for the backward scalar polariz- 
ability of the proton. The important finding was that the phe- 
nomenological a meson can be replaced by experimental in- 
formation on the nn continuum. The results of a more recent 
analysis are (a_Bi~/3 A /i) s = -5.6 and (a E i-f3Mi)t = 16.5, 
leading to a total value of abou t 10.9 in good agreement with 
the data dDrechsel et all 120031) . The importance of the t- 
channel contribution has also been found in a careful analy- 
sis of the new experimental data by ISchumacheij d2007l) who 
obtained (otEi — /?A/i)t = 15.2. We conclude that the polariz- 
ability of the nucleon is largely determined by the subprocess 
7 + 7 — * 7r + 7r, and therefore intertwined with correlations 
of the two-pion system and the polarizability of the pion. 



3. RCS data and extraction of the proton polarizabilities 

The pioneering experim ent in Compton scatter ing off the 
proton was performed by iGol'danski et all (Q~960). They ob- 
tained an electric polarizability a p E1 — (9±2), with a large un- 
certainty in the normalization of the cross section giving rise 



to an additional syste matical error of ±5. In a later experiment 
iBaranov et all (1 1975b used bremsstrahlung providing photons 
with energies up to 100 MeV. The data obtained by these au- 
thors were later reevaluated by DRs with the result a E1 w 12 
and P M1 ps —6. This outcome was much to the surprise 
of everybody, because one expected a large paramagnetic ef- 
fect of at least (3 p . dri sa 10 from the quark spin alignment 
in the N — > A(1232) transitio n. The first modern ex peri- 
ments were performed at Illinois dFederspiel et all 199 lb. fol- 
lowed by the work of the Saskatoon group (Halli n et al. . 1993; 
iMacGibbon et all 1 1995b . With tagged photons at 70 MeV< 
Ej < 100 MeV and untagged photons for the higher ener- 
gies, the latter group obtained a E1 = (12.1 ±0.8 ±0.5) 
and (3 P M1 = (2. 1 =F 0.8 =F 0-5). New pre cision measure- 
ments at MAMI dOlmos de Leon et all [2001 ) have been per- 
formed with tagged photons and the photon detector TAPS. 
The measured differential cross sections from various labo- 
ratories are shown in Fig. [17] as a function of the photon lab 
energy and at different scattering angles. The data have been 
compared to the results from four different types of DRs. The 
figure shows that the differences among the predicted results 
are hardly visible, except for the unsubtracted hyperbolic DR 
at ftab = 107°, because this angle is too much forward for 
this DR. A fit to all modern low-energ y data constrained by 
Baldi n's sum rule of Eq. d64l l yields dOlmos de Leon et all 
12001b 



a 



12.1 ± 0.3 T 0.4 ±0.3, 



0: 



El 

P M1 = 1.6 ±0.4 ±0.4 ±0.4. 



(74) 
(75) 



in units of 10~ 4 fm 3 and with errors denoting the statistical, 
systematical, and model-dependent errors, in order. This new 
global average confirms, beyond any doubt, the dominance 
of the electric polarizability a p E1 and the tiny value of the 
magnetic polarizability ii p Ml , which has to come about by 
a cancelation between the large paramagnetic s-channel 
contribution of the iVA spin-flip transition and the somewhat 
smaller diamagnetic t-channel contribution of the "pion 
cloud". The huge improvement by the new data is seen in 
Fig. [18] which displays the error ellipses in the a p E1 - f3 p M1 
plane as obtained from recent experiments. For further details 
of the experiments and their interpretation, see the review of 
ISchumache3d2005b . 

Much less is known about the spin sector, except for the 
forward and backward spin polarizabilities. The most recent 
values are 



To 



7tt 



-0.94 ±0.15, (76) 

-36.1 ± 2.2 (Olmos de Leon et al. 2001) 
-37.9 ± 3.6 fGaller et al. 20011 }(77) 
-38.7 ± 1.8 rSchumacher. 20051 



in units of 10~ 4 fm 4 . The small value for 70 in Eq. d76b 
was not measured by Compton scattering but has been de- 
termined by a sum rule based on the helicity-dependent ab- 
sorption cross sections, see section IVI.AI The upper line 
in Eq. ( 1771 ) gives 7^ as determined from low-energy data 
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FIG. 17 Differential cross section for Compton scattering off 
the proton as a function of the lab photon energy _E 7 and 
at different scattering angles f?iab. Solid lines: fixed-t sub- 
tracted DRs, dashed lines: fixed-t unsubtracted DRs, dotted 
lines: hyperbolic subtracted DRs, dashed-dotted lines: hyper- 
bolic unsubtracted DRs. All results are shown for fixed val- 
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37. The data a r e from loimos de Leo n et al. (2001) (full circles), 
iFederspiel ela/l il99lh (s q uares) , iMacGibbon et al\ i 19951) (tri- 
j 19921) (open circles). Figure from 



angles), and Zieger et al\ 



Drechsel et al. (2003). 



doi mos de Leon et al .,2001), in which case the error is domi- 
nated by the statistical plus systemat ical uncertainty, wh ereas 
the middle line refers to the work o f lGallerer aZ.I (f2001) who 
found that the model error prevails in the A region, and the 
lower line gives the weig hted average of several MAMI re- 
sults dSchumacherl [2005 ) . For all the other spin and higher 
order polarizabilities, both ChPT and DR predict small val- 
ues that can not be determined without dedicated polariza- 
tion studies. A new generation of experiments with polar- 
ized beams, polarized targets, and recoil polarimetry holds 
the promise to disentangle all the scalar and vector polariz- 
abilities of the nucleon and to quantify the pr oton's full spin 
respon s e to an external ele c tromagnetic field dBabusci et all 



1998a; iBeane et all 120051; 



Hilde brandt et all |2004 . The 



HIGS project dWellerl 120071) of a high-intensity beam with 
circularly polarized photons in an energy range up to 140- 
160 MeV is ideally suited to perform such experiments in 
the pion threshold region. Complementary investigations are 
planned in the first resonance region u s ing th e Crystal Ball 
detector at MAMI dArends and Schererl 120071) . We strongly 




+ 1 for different measure- 
3 dashed ellipses are ob tained 



FIG. 18 Contour plot of x' 
ments of and {3\j,^ . The 

from lOlmos de Leon et al] J200lh . Ifederspiel et al\ ( 1 199 lh . and 
IMacGibbon et al\ il995|) as indicated. Th e dashed area corresponds 
to the measurement of Zieg er et al\ j 19921) at 8 — 180°, the area be- 
tween the other two straight lines to the Baldin sum rule. The small 
ellipse draw n with the solid line is a com mon fit to all the above data. 
Figure from Olm os de Leon et oZl d200fh . 



believe that only a combination of these experimental projects 
will provide the "sharp knife " to extract the spin p olarizabili- 
ties in an unambiguous way (Pasquini et al. .1 120071) . 



4. RCS data and extraction of the neutron polarizabilities 

The experimental situation concerning the polarizabilities 
of the neutron is still quite unsatisfactory. The electric 
polarizability a] si can in principle be measured by scatter- 
ing low-energy neutrons on the Coulomb field of a heavy 
nucleus, whereas the magnetic polarizability (3 1 ^ I1 remains 
essentially unconstrained by such an experiment. This 
technique seemed to be very promising until t he beginning 
of the 1990's, when Schmiedmave r et al\ dl99ll) obtained the 
value a 1 ^ l = 12.6 ± 1.5 (stat) ± 2.0 (syst) by scattering 
neutrons with energies 50 eV < E n < 50 keV off a 208 Pb 
target. Shortly later, however, Nikolenko a nd Popovl d 19921) 
argued that the errors were underestimated by a factor 5. 
These finding s were confirmed by a similar experiment 
dKoester et al. , 19951) resulting in a% x = ± 5, and by 
lEnikeraZ.ld 19971) who obtained 7 < a n E1 < 19 after a further 
analysis of the systematic errors. 



The neutron polarizabilities can also be measured by quasi- 
free Compton scattering off a bound neutron and elastic scat- 
te ring on a deuteron . The former experiment was performed 
bv lRose et al.\ (Q990). Interpreted in conjunction with Baldin's 
sum rule, the result was < a r E1 < 14 with a mean value 
a Ei ^ 10-7- The large error bar arises from the fact that the 
Thomson amplitude vanishes for a neutral particle, and there- 
fore also the interference between this term and the leading 
non-Born amplitude is absent. It was therefore proposed to 
repeat such an experiment at higher energies and backward 
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angles for which the sensitivity to a 1 ^ 1 — fj^ [1 is highest. Be- 
cause the data analysis is very sensitive to final-state interac- 
tions and two-body currents, it was suggested to measure the 
polarizabilities of the bo und proton at the same time. The 
proton values obtained bv lWissmann et all d 19991) were quite 
promising, Q?e\—$m\ ~ 10.3±1.7 (stat + syst)±l.l (mod). 
The experiment was then exte nded to the neutron by the 
CATS/SENECA Collaboration dKossert et all 120021) . Data 
were collect ed for both deuterium and hydrogen targets and 
analyzed by iLevchuk and L'vo vl (120001) . The agreement be- 
tween the polarizabilities of free and bound protons was still 
quite satisfactory, and the final result for the (bound) neutron 
was 

=9-8±3.6(stat)i"(syst)±2.2(mod). (78) 

This value i s compatibl e with an earlier datum of the Saska- 
toon group dKolb et all 2000), obtained at similar energies 
and angles but with a much larger error bar, a\ x — (3 Ml « 
12. The comparison between proton and neutron demon- 
strates that there is no significant isovector contribution to the 
scalar polarizabilities of the nucleon. Unfortunately, the ex- 
perimental data from elastic photon scattering off a deuteron 
are much more prone to mo del errors. Such expe riments 
have been performed at SAL dHornidge et all H000) and at 
MAX-lab jL undin et al, 2003), and within the formalism of 
ILevchuk and L'vovl d2000l) the following results have been ob- 
tained: 



ft Ml — 



-4.8 ± 3.9 (Hornidge efaL 2000) 
-2.3 ± 3.4 (Lundin et al., 2003) 



•(79) 



Altogether these numbers speak for a very small value of the 
backward scalar polarizability, which is difficult to understand 
on theoretical grounds. The quasi-free Compton scattering 
experiments off a bound neutron have also provided a first 
glimpse at the backward spin polarizability of the neutron. 
Whereas the large pion-pole contribution is negative for the 
proton, it carries a positive sign for the neutron. The disper- 
sive contributions, on the other hand, are positive for both nu- 
cleons. As a result, we expect a large positive number for 7™. 
This is consistent with the value 



7£ = 58.6 ± 4.0 , 



(80) 



obtained from a fit to quasi-free Compton scattering off a 
bound neutron dSchumachen 120051) . 



B. Generalized polarizability of the nucleon at Q 2 > 

Virtual Compton scattering is formally obtained from real 
Compton scattering by replacing the incident real photon with 
a virtual photon 7*. It is realized by a subprocess of the reac- 
tion e + p — * e' + p' + 7. As displayed in Fig. [19] the real 
photon can be emitted by either the electron or the proton. The 
former process is called Bethe-Heitler (BH) scattering that can 
be calculated from QCD, whereas the latter process is referred 




(a) 




FIG. 19 Contributions to the reaction e + p — > e' + p' + 7, (a): 
virtual Compton scattering on the proton, (b): Bethe-Heitler process. 
Fhe blob in diagram (a) represents both nucleon intermediate states 
(Bom terms) and excited states of the nucleon (non-Born terms). 



to as virtual Compton scattering (VCS). Because the two pro- 
cesses lead to the same final state, the amplitudes add coher- 
ently, 



T 



vcs 



(81) 



The VCS amplitude 7~ vcs can be further decomposed into a 
Born and a non-Born contribution. For the Born contribution, 
the nucleon remains always in its ground state, and therefore 
this amplitude can be calculated once the (ground state) form 
factors of the nucleon are known. The non-Born term con- 
tains all the contributions with excited intermediate states, that 
is nucleon resonances, pion-nucleon scattering states, and so 
on. The physics interest is, of course, in the non-Born am- 
plitude, because it contains the information on the nucleon's 
internal structure in the form of generalized polarizabilities 
(GPs). These GPs depend on the virtuality Q 2 transferred by 
the virtual photon. The physics of VCS is visualized best if we 
consider the time-reversed version: As in RCS the real pho- 
ton plays the role of a quasi-static electromagnetic field that 
induces a polarization of the charges, currents, and magneti- 
zations whose spatial distributions are resolved by the virtual 
photon through variation of Q 2 . To lowest order in the en- 
ergy, VCS is determined by 6 independent GPs, which can be 
determined by measuring the interference between the Bethe- 
Heitle r and VCS amplitud es by means of angular distribu- 
tions dGuichon et all] 1995b and double-polarization asymme- 

m 



tries d Vanderhaegheri u 997l) . A word of caution for the reader 



familiar with the formalism of meson electroproduction: the 
cross section for the reaction e + N—>e' + N' + 'y does not 
take the form of Eq. d28b , which is based on particle produc- 
tion from the nucleon only. Instead, the final-state photon can 
be emitted from both the electron and the nucleo n. The pio- 
neerin g VCS experiment was made at MAMI bv lRoche et al\ 
(2000), and a first double-polarization experiment is under- 
way. Further experiments ha ve recently been perfor med at 
MIT/Bates for very small Q 2 (Bourgeo is et all 12006) and at 
JLab for large Q 2 dLaveissiere et a ll 12004 . From such ex- 
periments we find different spatial distributions for the dia- 
magnetism and paramagnetism in the nucleon. Furthermore, 
the planned double-polarization experiments will give a direct 
comparison with the spin polarizabilities predicted by ChPT, 
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which are free of low-energy constants at the leading order. 

1 . Kinematics and invariant amplitudes 
In the following we only consider the subprocess 

1 *(q) + N(p)^ 1 (q') + N(p'). (82) 

Because the space-like virtual photon has a "mass" q 2 = 
—Q 2 , the kinematic relations change with regard to the real 
photon case, Eq. 07] ). in particular 



s + t + u = 2M Z - Q z , 
s-u „ t-Q 2 



v — 



AM ' 



AM 



(83) 
(84) 



t = 2E' 1 (cos 6 lah ^E 2 +Q 2 -E 1 )-Q 2 , (85) 

with 6>i a b the lab scattering angle, and E' and E 1 the lab 
energies of the real and virtual photon, respectively. In 
the following we choose v, t, and Q 2 as the independent 
variables. 

The VCS Compton tensor is constructed as for RCS, 
Eq. d58l . except that the polarization four-vector of the vir- 
tual photon has 3 independent components, that is the helici- 
ties A = ±1 (transverse polarization) and A = (longitudinal 
polarization). The VCS tensor H^ v can be expanded in a ba- 
sis of 12 independent tensors with amplitudes depending on 3 
variables, 



12 



(86) 



The number 12 is given by the possible choices for the helic- 
ities in the initial and final states, namely 3 x 2 x 2 x 2 di- 
vided by two because of parity invariance. The same consider- 
ation yields 8 for RCS, but this number is further reduced to 6 
independent combinations by time-reversal, which of course 
does not apply for VCS. It is possible to find a special ten- 
sor basis such that each term is gauge invariant, even under 
crossing and free of k inematical singularities and constraints 
dDrechsel et all [T998I) . Furthermore, only 6 amplitudes con- 
tribute for Q 2 — > 0, because 4 tensor structures and 2 ampli- 
tudes vanish in this limit. The result are 6 relations between 
the VCS amplitudes Fi and the RCS amplitudes A». 



2. Generalized polarizabilities 

If the emitted photons have small energies, the VCS exper- 
iments can be analyze d in terms of a low-e nergy expansion 
(LEX) as proposed bv lGuichon et all (11995b . In this approx- 
imation the non-Born part of the amplitudes is expanded in 
E!y, and only the linear term is kept. This reduces the multipo- 
larities of the emitted photon to electric and magnetic dipole 
radiation. Furthermore, the GPs are given by linear combi- 
nations of the amplitudes at threshold (u = 0, t = —Q 2 ), 



which contains the definition of the polarizabilities for RCS 
in the limit Q 2 — > 0. Let us next discuss the multipole de- 
composition of the non-Born VCS tensor of Eq. ( |86i > 
at small real photon energy, q' — > 0, but for arbitrary three- 
momentum q = | q cm | of the virtual photon. For this purpose 
we denote the GPs by P( M ' £',MC)s dQuichon et a/.lll995l) . 
In this notation, C refers to the angular momentum and M. 
to the electric (E), magnetic (M), or longitudinal (L) na- 
ture of the virtual photon, with the primed variables denoting 
the real photons. Furthermore, the quantum number S dif- 
ferentiates between the spin-flip (S = 1) and non spin-flip 
(S = 0) character of the hadronic transition. Within the LEX 
we may use the dipole approximation, £' = 1. With this as- 
sumption, the conservation of ang ular momentum and parity 
restricts the number of GPs to 10 dGuichon et q/.l[l995l) . Four 
more constraints are provided by nucleon crossing combined 
with charge conjuga tion symm etry, which leaves 6 indepen- 
dent GPs dDrechsel et 0/1,119981) . 



p(Ll,Ll)0 ((?); p(Ml,Ml)0 {q) 
p(Ll,Lin (qh p{Ml,Ml)l {q) 
p(Ml,£2)l/ ? \ p(il,M2)l/A 



(87) 



We note that the transverse electric multipoles have been elim- 
inated from the above equations, because they differ from the 
longitudinal multipoles only by terms of higher order in q. In 
the limit q — > one finds the f ollowing relations be tween the 
VCS and RCS polarizabilities dDrechsel et all 1 19981) : 



p{Ll,LX)0 



p(Ml,Ml)0 



V2 



V8 



a ei 



Ml 



%/3a e m 

p(Ll,Ll)l ^ Q) p(Ml,Ml)l _ >Q) (g8) 

v/2 



p(Ll,M2)l 
p(Ml,L2)l 



3V3 a t 



73 ■ 



(72 + 74) 



In order to connect the scalar VCS and RCS polarizabilities, 
we introduce the definitions 

vei(Q 2 ) = -^a^P^^^Q 2 ), (89) 

/WQ 2 ) = -yfa em P (M1 ' M1) °(Q 2 ), (90) 

withaEi(O) = a E i and/?Mi(0) = (3mi as measured by RCS 
according to Eqs. d74i > and d75l >. 



3. Theoretical developments 

For the given tensor basis, the associated non-Born VCS 
amplitudes F" B , (i = 1,...,12) are free of kinematical singular- 
ities and constraints, and even under crossing. Assuming fur- 
ther an appropriate analytic and high-energy behavior, these 
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amplitudes fulfill unsubtracted dispersion relations in the vari- 
able v and at fixed t and Q 2 , 



ReiT (Q 2 ,",*) = ^T le {Q\v,t)-F?{Q\v,t) 



°° dv , u'lmF t (Q 2 lV ',t) 



(91) 



We recall that the Born amplitudes Ff are given by diagrams 
with nucleons in the intermediate state, whereas the pole 
amplitudes Ff° Xe are obtained from the Born amplitudes at 
the pole position, that is, with all numerators evaluated at 
the pole. Furthermore, Im Fi are the discontinuities across 
the s-channel cuts, starting at the pion production threshold 
u = ma + (2ml + t + Q 2 )/(^M). 

Besides the absorptive singularities due to physical inter- 
mediate states, one might wonder whether additional singu- 
larities like anomalous thresholds can contribute to the dis- 
persion integrals. The latter arise when a hadron is a loosely 
bound system of other hadronic constituents which can go on- 
shell, thus leading to so-called triangular singularities. How- 
ever, it was shown that within the strong confinement of QCD, 
the quark-gluon structure of hadrons does not give rise to addi- 
tional anomalous thresholds (Jaffe and Mendel fl 9921: lOehmel 
11995b . and that possible quark singularities turn into hadron 
singularities as described through an effective field theory. 
Therefore, the only anomalous thresholds arise for hadrons 
which are loosely bound systems of other hadrons, as for ex- 
ample the £ particle in terms of a A-ir system. Such anoma- 
lous thresholds are absent for the nucleon, and therefore the 
imaginary parts in Eq. ( |9TT l are only given by absorptive ef- 
fects due to ttN, ttttN, and heavier hadronic states. Of course, 
Eq. d9~n > is only valid if the amplitudes drop fast enough such 
that the integrals converge. The hi gh-energy behavior of the 
amplitudes Fi was investigated by Pas quini et all (1200 1 ab in 
the Regge limit (y — * oo, t and Q 2 fixed). As in subsec- 
tion IIV.A.2I the dispersion integrals diverge for two ampli- 
tudes, Fi and F$ in our notation. These amplitudes are domi- 
nated by the t-channel exchange of a and 7r° mesons, respec- 
tively. As long as we are interested in the energy region up 
to the A(1232), we may saturate the s-channel dispersion in- 
tegral by the ttN contribution, choosing v max w 1.5 GeV as 
upper limit of integration. The asymptotic contribution to F5 
is saturated by the pion pole and therefore independent of v, 



F™(v,t,Q 2 ) 



(Q 2 ' 



Me 2 t 



(92) 



with a monopole form factor from the ir° — > 7 + 7 decay. As 
a result a^Q 2 ) ~ F 5 as (0, -Q 2 , Q 2 ) has a dipole form, with 
ag 1 (0) known from RCS. Although the pion pole contribu- 
tion is certainly dominant, there may be other effects such as 
more-pion and heavier intermediate states. In view of our still 
limited knowledge on these reactions, we simply parameterize 
the Q 2 dependence in a dipole form with a parameter A a , 



a 



li(Q 2 



km 



(l + QW) 5 



(93) 



In the same spirit we also estimate the contribution of the a 
meson by a dispersion relation in t at v = 0, 



Ff(Q 2 ) 



1 



dt 



,ImtFi(0,i / ,Q 2 



t' + Q 2 



(94) 



with Imt-Fi as determined from the t-channel reaction 7 + 
7* — > 7r + 7r — > NN. The result of this calculation can a lso 
be parameterized by a dipole form dPasquini et a/.luOOlal) : 



Pmi{Q 2 ) 



(1 



r/A 2 ) 2 



(95) 



with « 0.4 GeV and /3 as (0) known from RCS. We note 
that A/3 is small compared to the parameter A^ = 0.84 GeV 
of Eq. ( |37| |. which gives the scale of the nucleon's magnetic 
form factor, that is, the asymptotic diamagnetic polarization is 
related to surface phenomena as expected from the pion cloud. 



4. Experiments and data analysis 

At small 3 -momentum q' of the emitted real photon, the 
measured cros s secti on can be analyzed through the LEX of 
iGuichon et all dl995h . This expansion is based on a low- 
energy theorem (LET) stating that the radiative amplitude for 
point-like particles diverges like 1/q' for q' — > 0, whereas the 
dispersive amplitude vanishes like q' in that limit. As a conse- 
quence the spin-averaged (unpolarized) square of the matrix 
element takes the form 



\M\ 2 = 



A- 



t\2 



■Ao + O(q'). 



(96) 



with coefficients _4_ 2 an d A^i that are fully described by 
the Bethe-Heitler (BH) and Born terms, which can be cal- 
culated from QED once the proton (ground state) form fac- 
tors are known. The next order term _4 contains contribu- 
tions from the BH and Born terms but also an interference 
between the 0(1/ q') contribution of BH plus Born ampli- 
tudes and the leading term of the non-Born amplitude 0(q'), 
which is proportional to the GPs. This interference term can 
be expre ssed by the structure functions Pll(q), Ptt(q), and 
P L T(q) (IGuichon et all 1 1995b . 



yf xp - ^C +B - 2^ 2 ^ vi [eP LL (q) - P TT (q)} 



—V3W2e{l + e)P LT (q)\, (97) 



U>2 



with K2, Vi ,«2> an d W3 kinematical functions depending on e, 
q, and the polar and azimuthal cm angles, 9 cm and <f> cm , respec- 
tively. Furthermore, go is the cm energy of the virtual photon 
in the limit q' — * 0. The 3 structure functions of Eq. ([97]) 
can be expressed by the GPs as follows dGuichon et ail 1995; 



29 



iGuichon and Vanderhaeghen[ |l998): 

Pel = -~2V6MG E P iL1 ' L1)0 , (98) 



P ll -P tt /e (GeVO 



P LL -P TT /e (GeV" z ) 



?! 

P ^ = \/I^g £ p( m w + ^GmP (l1 - l1)1 , 

with and Gm the electric and magnetic nucleon form 
factors. 

In Fig. [20] we compare the measured response functions 
to the predictions of DR (left column) and HBChPT (right 
column). The response function Pel — Ptt/s is displayed 
in the upper panels of this figure. According to Eqs. (l89l and 
j98l , Pll is directly proportional to the scalar GP aEi(Q 2 ), 
whereas Ptt contains only spin GPs. As discussed in 
subsection IIV.A1 the dispersive and asymptotic contributions 
to a ei have the same sign at the real photon point, which 
leads to a large total value. However, o>ei(Q 2 ) drops rapidly 
as function of Q 2 . The difference between the solid and 
the dashed lines is due to the spin GPs whose importance 
rises with Q 2 . T here is a gene r al agr eement between the 
results from PR JPasauini et all 1200 lab and the HBChPT 
dHemmert et all l2000l 1 1997a). however the spin GPs turn 
out much larger in the latter approach. The lower row of 
Fig- E01 gives the same comparison for the response function 
Plt, which contains both the scalar magnetic polarizability 

p (A/i,Mi)o [3 M1 (Q 2 ) and the spin GP p\ L ^)\ As 

shown in subsection IIV.AI (3 mi (0) is the sum of a large 
dispersive (paramagnetic) contribution, which is dominated 
by A(1232) excitation, and a somewhat smaller asymptotic 
(diamagnetic) contribution with opposite sign. Moreover, 
we expect that the diamagnetic contribution is largely due to 
pionic degrees of freedom, and therefore of longer range in 
r-space than the paramagnetic component. This expectation 
is corroborated by the minimum of Pet at Q 2 « 0.05 GeV 2 : 
the (positive) diamagnetic component of Pet decreases 
faster in Q 2 -space than the (negative) paramagnetic term, 
and therefore Plt decreases over a small Q 2 region to the 
minimum, from whereon the form factor effects lead to a 
rapid approach towards zero. Although the full results of 
DR and HBChPT agree qualitatively, there is again a large 
difference in the spin-dependent sector. 

In the region between pion threshold and A-resonance 
peak, the sensitivity to the GPs is much enhanced, because 
the contributions of the GPs interfere with the rapidly rising 
amplitude of the A-resonance excitation. It is of course 
not possible to extend the LEX to these energies, but the 
dispersive approach is expected to give a reasonable frame 
to extract the GPs. When crossing the pion threshold, the 
VCS amplitude also acquires an imaginary part due to the 
opening of the irN channel. As an interesting result, single 
polarization observables appear above pion threshold. A 
particularly relevant observable is the electron single-spin 
asymmetry (SS A), which is obtained by flipping the electron 
beam helicity dGuichon and VanderhaegherJ, 1 1998b . The 
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FIG. 20 Comparison between the unpo larized VCS structure func- 
tions calculated by dis persion relations dPasauini et all 12001 ah (left 
column) and HBChPT Memmert et q/.Ll2000ll 1997hF at Q(v 3 ) (right 
column). Upper row: Result for Pll — Ptt/s with e = 0.62 (solid 
lines) and e = 0.9 (dashed-dotted lines) compared to the result for 
otEi only and e = 0.62 (dashed lines). The dispersive results for 
e = 0.62 and s = 0.9 are obtained with A a = 1.79 GeV and 
A Q = 0.7 GeV, respectively. Lower row: Results for Plt (solid 
line) compared to the result for /3m i only (dashed li ne). The data are 
from | 01mos de Leon et al\ d200ll) (open t riangl es), iBourgeois et al\ 
(2006) (full triangles), an d iRoche et al\ <2OO0l) (squares). Figure 
from lDrechsel et al\ J2003h updated by B. Pasquini. 



main source of the SSA is an interference between the 
(real) Bethe-Heitler and Born amplitudes and the imaginary 
part of the VCS amplitude. Because the SSA vanishes 
in-plane, its measurement requires an out-of-plane experi- 
ment. Such an experiment h as recently been performed at 
MAMI (Bensaf a et all 12007b . The measured asymmetry at 
W=1.19 GeV and Q 2 = 0.35 GeV 2 is displayed in Fig. ED 
and compared to the predictions of dispersion theory. The 
figure shows a rather weak dependence of the asymmetry on 
variations of the GPs. Therefore, a measurement of the SSA 
provides an excellent cross-check of the dispersive input, i.e., 
the imaginary parts of the tt N multipoles, in particular by 
studies of the A region by VCS and pion electroproduction in 
parallel. 

At larger virtuality, the VCS process has been investi- 
gated by the Hall A Collaboration at JLab, and data have 
been obtained at Q 2 = 0.92 GeV 2 and Q 2 = 1.76 GeV 2 
dLaveissiere et all 120041) . A reasonable description of 
these data is obtained by the values A a — 0.71 GeV and 
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FIG. 21 Beam single spin asymmetry (SSA) for VCS as function 
of the photon scattering angle. The dispersive predictions are shown 
for different pairs of (A Q , Kg) given in GeV. Solid line: (1.4, 0.7), 
dashed-dotted line: (1.4, 0.4), lon g-dashed line: (0.8, 0.7), short- 
dashed line: (0.8, 0.4). Figure from lBensafa et al\ 120071) . 



hp = 0.51 GeV shown by the solid lines in Fig. [22] We 
note that A Q corresponds to a spatial range comparable 
to the nucleon's charge distribution. However, the best 
fit value for is substantially lower, indicating that the 
diamagnetism is related to pion cloud effects at distances 
above 1 fm. Subtracting the spin-dependent terms according 
to the dispersion predictions, we obtain the Q 2 dependence of 
the scalar GPs shown in Fig. [23] It is obvious that the electric 
GP aEi(Q 2 ) is dominated by the asymptotic term, which 
however can not be described by a single dipole form over 
the full Q 2 range. The magnetic GP (3mi(Q 2 ) clearly shows 
a characteristic maximum at Q 2 w 0.05 GeV 2 , which comes 
about by cancelation between the positive paramagnetic A 
contribution and the negative diamagnetic contribution of 
the ^-channel 7T7t exchange. By Fourier transforming the 
GPs aEi(Q 2 ) and (3mi(Q 2 ) in the Breit frame, one obtains 
the spatial distribution of the ind uced electric polar ization 
and magnetization of the nucleon (iL'vov et al. 1 120011) . The 
emerging picture is as expected from a classical interpretation 
of diamagnetism. Due to the external magnetic field, pionic 
currents start circulating in the nucleon and give rise to an 
induced magnetization opposite to the applied field. At 
distances r > l/m^, the diamagnetic effect dominates 
and the Fourier transform /3mi(0 takes negative values, 
whereas the paramagnetic contributions prevail at the smaller 
distances giving rise to positive values of (3mi(t) in the 
interior of the nucleon. As the momentum transfer Q 2 
increases, the negative contribution due to the long-range 
pion cloud vanishes fast and hence /3mi{Q 2 ) increases. This 
nicely explains the positive slope of (3mi{Q 2 ) at Q 2 = 
and the maximum at Q 2 0.05 GeV 2 as indicated by the 
experimental data. 



FIG. 22 Results for P LL - P T t /s (left panel) and P LT (right panel). 
Dashed lines: dispersive ttN contributions. Dotted lines: asymptotic 
contributions with Aq, = 0.71 GeV and Kg = 0.51 GeV. S olid lines: 
full resu lts. The data are f rom | Qlmos de Leon et al\ i200lh (open tri - 
angles), iBourgeois et al\ teOOo) (full triangles), iRoche et al\ d2000h 
(squares), and lLaveissiere et al\ J2004I) as obtained by the LEX (open 
circles) and DR (full circles) analysis. The inner error bars describe 
the statistic al error, the outer erro r bars include systematical errors. 
Figure from lDrechsel et al\ ( 120031) . updated by B. Pasquini. 

According to Eqs. d97| i and ((98), the unpolarized VCS ex- 
periment gives a ccess to only 3 combin ations of the 6 GPs. 
As was shown bv l VanderhaegherJ d 1 997b . it takes experiments 
with polarized lepton beams and polarized targets or recoil 
nucleons to measure the remaining 3 GPs. These double- 
polarization observables require measuring the cross sections 
for a definite electron helicity h and recoil (or target) proton 
spin orientation parallel and opposite to a specified axis. As 
shown for the unpolarized cross section by Eq. d96l ), also the 
polarized squared amplitude has a low-energy expansion, and 
again the GPs are obtained from the term O (1). This term 
contains the structure functions Pfrp( q), P^jq), and P T \,(q ), 
which are related to the spin GPs by ( Vanderhaeghen, 1997) 

PI = ^. GmP (L1,L1)1 _ ^M±G E p(Ml,Ml)l (99) 

2q Q 

p'lt = ~IqG m P {L1 - L1)1 + ^G £ P^ 1 - M1 > 1 , 
2 Qq 

P ™ = ^ GM ( P(L1 ' L1)1 -\/l* P(M1 ' i2)1 )- 

While P[ T and P^ T can be accessed by in-plane kinematics, 
P L y requires an out-of -plane measurement. 

In Fig. [24]we compare the results of DR and HBChPT for 
the double-polarization observables, with longitudinally po- 
larized electrons and recoil proton polarization either along 
the virtual photon direction (z -direction) or in the reaction 
plane and perpendicular to the virtual photon (x-direction). 
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FIG. 23 Left panel: the electric GP olei (Q 2 ) obtained from the DR 
formalism with A Q = 0.7 GeV (solid line) and A Q = 1.79 GeV 
(dashed-dotted line). The solid line is the sum of the asymp- 
totic (dashed line) and dispersive (dotted) contributions. Right 
panel: the magnetic GP /3mi(Q 2 ) obtained from the DR formal- 
ism with A = 0.51 GeV (solid line) and A = 0.63 GeV 
(dashed-dotted line). The dashed-dotted line is the sum of the 
asymptotic (dashed line) and dispersive (dotted l ine) contributions. 
The data are from lOfmos de Leon et al.\ (l200ll) (open triangles), 
iBourgeois et a/.N2006h (full triangles). iRoche et a/.ld2000l) (squares), 
and lLaveissiere et al\ ([2004) as obtained by the LEX (open circles) 
and DR (full circles) analysis. The inner error bars describe the sta- 
tistical error , the outer error bars a lso include systematical errors. 
Figure from Laveissiere et al. ( 2 0041) . 



The large but well-known asymmetries from the Bethe-Heitler 
and Born terms have been subtracted in thi s figure in or- 
der to highlight the differences between DR dPasauini et all 
2001 ah and HBChPT at C(p 3 ) dHemmert et all |2000|) . We 
note that the latter approach yields significantly larger ef- 
fects due to higher predicted values for the spin GPs. Al- 
though these double polarization observables are tough to 
measure, a first test experiment is underway at MAMI. Con- 
trary to the scalar polarizabilities, the spin-flip GPs are still 
unknown territory. In Fig. [25] we compare the dispersive re- 
sults for the spin-flip GPs with the predictions of the non - 
relativistic constituent qu ark model (jPasquini et all 2001b). 
the HB ChPT to Q(v 3 ) dHemmert et all l2000t 1 1997al) and 
C(p 4 ) dK^o" and Vanderhaeahen, 2002), and the linear a- 
model (M etz and Drechse]Lll997l) . We refrain from comment- 
ing on the theoretical predictions which clearly open a wide 
range of values for the spin polarizabilities. An absolute 
must for further progress are dedicated experiments with a 
large sensitivity to the spin-dependent GPs. Such experiments 
are (I) unpolarized VCS with variation of the transverse pho- 
ton polarization e in order to separate the response functions 
Pll and Ptt and (II) double-polarization experiments as dis- 
cussed above. 



C. Polarizability of mesons 

Mesons are systems of a quark and an anti-quark and, there- 
fore, theoretically simpler to describe than baryons. Let us set 
the scene with the classical picture of two charges bound in 
a quark-antiquark potential in the presence of an additional 
static electri c field. For this system we can derive the follow- 
ing relation dWalchetil2006l) : 




V«leg) 



Vdeg) 



,7T + 

*E1 



(100) 



with C = 1/6 the effective charge of the system and a q q m 5 
derived from the heavy quark potential, of course somewhat 
outside its applicability. Furthermore, R is a characteristic 
dimension of the system, for example the equivalent charge 
radius. With the pion rms radius of Eq. d49l , we obtain 

R = VW3yp)7 = 0.86 fm, and as a result d^i = 3.2, 
here and in the following in units of 10~ 4 fm 3 . Comparing 
these numbers with the results of subsection IIV.AI we find 
that the pion is a "dielectric medium" with e w 1.001, i.e., 
even more rigid than the nucleon. 



FIG. 24 Predicted double-polarization asymmetries as function of 
the photon scattering angle # 77 and for the following fixed kine- 
matic values: q = 600 MeV, q' = 111.5 MeV, e = 0.62, and 
$ = 0°. In order to highlight the model dependence, the (known) 
Bethe-Heitler and Born contributions to the asymmetry have been 
subtracted. Solid lines: Results of dispersion relations for A a = 
1 GeV and A = 0.6 GeV dPrechsel e?a/ll2003l), dott ed lines: pre- 
dictions of HBChPT at Q(y 3 ) dHemmert et all 1 20001) . and dashed 
lines: HBChPT at 0(p 4 ) dKao et fl/.Ll2004h . See text for further ex- 
planations. 



The pion polarizabilities have been calculated in ChPT at 
the two-loop order, 0(p 6 ). Contrary to the situation of the 
nucleon, no "matter fields" with their own mass scale are 
present, and therefore the calcul ations can be performed in the 
origin al formulation of ChPT dGasser and Leutwylerl Il984l 
1985). This makes the following predictions for the polar- 
izabil ities a very significant test of this theory dGasser et all 
2006): 
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FIG. 25 The spin-flip GPs (without the 7r°-pole contribu- 
tion) from s everal calculat i ons. Solid lines: dispersive nN 
contributio n dPasguini et all l2001al) . short-dashed lines: C(p 3 ) 
HBChPT faemmert et all l2OO0t) . long-dashed lines: 0(p 4 ) 
HBChPT jKao et all |2004) . dashed-dotted lines: linear a model 
jMetz a nd Drechs ef^ 1 19961) dotted lines: non-relativistic CQM 
( Pasquini et al. , 2001b). For visibility, the tiny CQM results for 
P (Li Li)i md p(Mi,L2)i ^ multiplied by a f actor 100 . Figure 

from lDrechsel et all (2003) updated by B. Pasquini. 



of the experiment is dAhrens et al. 1 120051) 
= 11.6 ±1.5, 
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±0.5 n 



(103) 



which is at variance with the prediction of iGasser et all 
(2006) by two standard deviations. In view of the theoretical 
uncertainties from the fact that the Compton scattering is off 
a bound pion, the deviation from theory is an open problem. 
In particular we point out that the model error in Eq. d 1 03b is 
estimated by comparing the analysis with 2 specific models. 
This does not exclude that a wider range of models will lead 
to a larger model error. Considering the fact that the scattering 
is off a "constituent" pion in the nucleon, we may attribute 
the deviation to binding effects. For example, as suggested 
by Eq. ( 1 100b . an increase of the bound pion radius by 20 % 
would give a hand-waving explanation for the experimental 
data. Because the pion polarizability is extremely important 
for our understanding of QCD in the confinement region, it 
is prerequisite to check the given arguments by a full-fledged 
ChPT calculation of the reaction 15(7, 7'7r + n) . 

The second method to determine the polarizability, the Pri- 
makov effect, has been studied at Serpukhov with the result 
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Ml 



13.6 ± 2.8 stat ± 2.4, 



yst ? 



(104) 



in agreement with the value from MAMI. Recently, also the 
COMPASS Collaboration at CERN has investigated this reac- 
tion. However, at this time the analysis is still in a too prelim- 
inary stage to include the result. Unfortunately, the reactions 
e + e~ <^ 77 ^ 7r + 7r~ have led to even more contradictory 
results in the range 4.4 < a^ !1 < 52.6, as listed in the work of 
Gass er et al. (2006). In conclusion one has to wait for an im- 
proved analysis and possibly also independent experimental 
efforts before final conclusions can be drawn. 



The very small value predicted by Eq. ( 11011 ), that is Baldin's 
sum rule applied to the pion, makes a measurement of this 
observable close to impossible. The experiments are therefore 
analyzed with the constraint a 7 ^ ~ —ft mi- 

Unfortun ately, the exp e rimen tal si tuation is rather co ntra- 
dictory, see lAhrens et all d2005l) and IGasser et al. 1 (120061) for 
recent reviews of the data and further references to the exper- 
iments. There exist basically three different methods to mea- 
sure ctEi'- (I) the reactions e + e~ 77 ^ 7r + 7r~, (II) the Pri- 
makov effect of scattering a relativistic pion in the Coulomb 
field of a heavy nucleus, and (III) the radiative pion photo- 
production, p(j, r y , Tr + n), which contains Compton scattering 
on a (bound) pion as a subprocess. The latter reaction was 
recently investigated at the Mainz Microtron MAMI, by use 
of a kinematically optimized set-up consisting of the back- 
ward photon detector TAPS, a forward tt + detector, and a neu- 
tron detector realized by a large scintillator wall of dimensions 
3x3x0.5 m 3 . The largest error of this measurement is due to 
the systematic error of the neutron efficiency. The final result 



V. EXCITATION SPECTRUM OF THE NUCLEON 



A. Threshold production of mesons 



As outlined in subsection III. Dl the threshold photoproduc- 
tion of mesons provides a significant test of our theoretical 
understanding, because only few partial waves contribute and, 
therefore, all relevant multipoles can be directly determined 
by the experiment. The case of neutral pion photoproduc- 
tion on the proton, 7(g) + p(pi) — ► ir a (k) + p(p2) is of par- 
ticular interest. For the s-wave threshold multipole Eq + of 
this reaction, several authors had derived a low e nergy theo- 
rem (LET) based on current algebra and PCAC dDe Baenstl 
119701: IVainshtein an d Zakharo^ [l97l). According to the the- 
orem, the leading terms of the threshold multipole were di- 
rectly determined by the Born diagrams, evaluated with the 
pseudovector pion-nucleon interaction. However, this predic- 
tion had to be revised in the light of surprising experimental 
evidence. The reason for the discrepancy between the the- 
orem and the data was first explained in the framework of 
ChPT by pion-loop corrections. An expansion in the mass 
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FIG. 26 The real (Re) and imaginary (Im) parts of the s-wave ampli- 
tude E p+ for 7T° photopr oduc tion at threshold en ergies. The MAMI 
data of iBeck et al\ jl990j) and Fuchs et al. (1996) are represented by 
circles, the SAL data of Ber gstrom et al . ( 1996) by tr iangles . Dashed 
lines: predictions of ChPT at C(p 3 ) teernard et a/.Lll9 96a bl. solid 
lines: results from dispersion relations faanstein et a/.lll997h. The 
solid h oriz ontal line at about -2.2 shows the prediction of lDe Baenstl 
( 1970) and lVainshtein and Zakharovl Jl972h . Figure by courtesy of 
R. Beck. 
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= TOtt/M pa 1/7 yielded the result dBernard et all 
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(i - (i 



IM 



87rm T y •- 2 r ~ 16/2 ' 

(105) 

where g^N is the pion-nucleon coupling constant and 
f v w 93 MeV the pion decay constant. We observe that 
Eq + (tt°p) is proportional to fi, which suppresses this reaction 
relative to charged pion production. The first and the second 
ter m on the r hs of Eq. (1105|> are identical to the "LET " 
of iDe Baenstl dl970l) and IVainshtein and Zakharovl dl972l) , 
which however has to be corrected by the third term on the 
rhs. Although this loop correction is formally of higher order 
in fi, its numerical value is of the same size as the leading 
term. 

The energy dependence of E 0+ (ir°p) is shown in F ig. |26] 
The discrepancy between prediction s of IDe Baensl dl970h 
and IVainshtein and Zakharovl d 19721) and the experimental 
data obtained at the Mainz Microtron MAMI and at SAL 
(Saskatoon) is apparent. Furthermore, the real part of the am- 
plitude shows a characteristic "Wigner cusp" at the threshold 
for charged pion production, which lies about 5 MeV above 
the 7r° threshold. This cusp in the real part is related to the 
sharp rise of the imaginary part at the second threshold. The 
physical picture behind the large loop correction is based 
on (I) the high production rate of charged pions and (II) the 
charge-exchange scattering between the nucleon and the slow 
7r + in the intermediate state, which leaves a ir° in the final 
state. However, the direct experimental determination of the 



FIG. 27 The angular distribution of the photon asymmetry S for 
the reac tion ppy, 7r°)p . Da shed line: results of dispersion re- 
lations (jHanstein et all 1 1997b . dotted line: prediction of ChPT 
teernard et all Il996bl) . full li ne: empirical fit to the data with 
Eqs. (25) and (|26j. Figure from lSchmidt et al.\ d200ll) . 



imaginary part will require double-polarization experiments 
with linearly polarized photons and polarized targets. The 
excellent agreement between ChPT and the data for Eq + is 
somewhat flawed by the fact that higher order diagrams are 
sizeable, that is, the perturbative series converges slowly and 
low-energy constants appearing at the higher orders reduce 
the predictive power. 

For a more quantitative presentation of the results, the Eq + 
amplitude was param eterized as the sum of a direct and a 
charge-exchange term dBernstein ef aZ.L 1 19971) . 

E 0+ (ir°p) = A (q) + ia^+^o A + k n + (106) 
= A (q)+ipk n+ , (107) 

with Ao(q) — a® + a\{q — q t h T ) and A + describing the 
neutral and charged pion production in the absence of the 
charge exchange reaction, a^+^o the scattering length for 
charge exchange, and fc 7r + the momentum of the charged 
pion appearing in the intermediate state. This leaves the 3 fit 
parameters a > a i, an d P m oder to determine the s wave. As 
discussed in subsection III.D1 the unpolarized cross section 
gives information on only 2 of the 3 p-wave amplitudes. A 
complete experiment therefore requires measuring a further 
observable, e.g., the photon asymmetry £ of Eq. d26l i. The 
result obtained at MAMI is depicted in Fig. [27] 



The data of Schmidt et al. (2001) for the s- and p-wave am- 
plitudes are compared to the results of ChPT and dispersion 
relations in Table |U The predictions of ChPT for the p waves 
are in good agreement with the data, whereas P3 comes 
out much too small in dispersion theory. Since the latter 
approach is mainly based on input from the imaginary parts 
of the multipoles M\ + and Mi_ in the resonance region, this 
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FIG. 28 The total cross section for the p(e, e'p)ir° reaction as 
function of Q 2 for several values of the excitation energy AW = 
W — M — m„o . The da ta at Q 2 = 0.05 GeV 2 a nd 0.10 GeV 2 are 
from lMerkel et al\ d2002h and lDistler eTaZl dl998h. respectively. The 
solid line represents the prediction of ChPT dBernard et aZl Il996ch 
and the dashed line th e result of the phenomenological m odel MAID 
dDrechselef"a/U l999). Figure from Merkel et al. (2002). 



failure may indicate that even the structure of the low-lying 
resonances A(1232) and TV* (1440) is not yet completely 
unraveled. 





(Schmidt et al., 2001) 


ChPT 


DR 


£ 0+ (7T U ) 
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(Ivan den Brink et all 1 1995b and MAMI dDistler et all 1 1998b 
for Q 2 — 0.10 GeV 2 provided another confirmation of ChPT 
although at the expense of 2 new low-energy constants, 
which were fitted to the data. In order to further check this 
agreement, data were also taken at the lower momentum 
transfer Q 2 = 0.05 GeV 2 dMerkel et ail 12002b . The total 
cross section obtained by these measurements is compared 
with the predictions of ChPT in Fig. [28] which shows the 
total cross sections near threshold as function of Q 2 . The 
comparison was made on the basis of the total cross sections 
in order to eliminate all possible systematic and model errors 
connected with the separation in longitudinal and transverse 
parts. It is apparent that the Q 2 dependence of the data can 
not not be fully described by theory. We consider this an 
important issue that deserves further investigations. 

We c onclude this subse ction by presenting some recent re- 
sults of lWeis et al\ {2007). Whereas former experiments were 
only sensitive to the real part of the amplitudes, these au- 
thors also determined the fifth structure function (XT') given 
in Eq. ([28). This function can only be measured with po- 
larized electrons and out-of-plane, i.e., for finite values of 
the pion azimuthal angle with regard to the electron scatter- 
ing plane. Furthermore, its multipole decomposition is of the 
formlm (Lg + M\+ + - • • ), i.e., this function contains informa- 
tion on the phase of the s -wave ampl i tude. With the shorthand 
notation dui/dVl*^ = q, . IWeis et al\ d2007l) separated the par- 
tial cross sections <7rj = ot ± ott, and <jlt as well as 
the beam asymmetry A' LT corresponding to <j' lt . The result 
is displayed in Fig. [29] We observe that only the dynami cal 
Dubna-Mainz-Taipei (DMT) model dKamalov et aZ.ll200ll) is 
able to fully describe the experiment, in particular its predic- 
tion for the helicity asymmetry is right on top of the data. Such 
dynamical models start from a description of the pion-nucleon 
scattering phases by a quasi-potential, which serves as input 
for an integral equation to account for multiple scattering. In 
this sense the model contains the loop corrections to an arbi- 
trary number of rescattering processes, and is therefore per- 
fectly unitary, albeit on a phenomenological basis that may 
violate gauge invariance to some extent. 



TABLE I Experimental results of ISchmidt et al\ d200lh for 

Eo+(iv p) at the ty° and 7r + fhresholds in units of 10~ 3 /m 7r +, 

four combinations of the (reduced) P-wave amplitudes in units B. Nucleon resonances and meson production 

k ■ 10~ 3 /rnt+, and the parameter f3 (with statistical and systematic 



errors, in order) compared to the predicti ons of ChPT dBernard et ail 
ll996al lbT) and dispersion relations (DR) dHanstein et all 1 1997b . For 
the definition of the reduced P-wave amplitudes see section HTDl 



The great success of ChPT for photoproduction at threshold 
was a strong motivation to extend the experimental program 
to electroproduction. Because the virtual photon has an 
additional longitudinal component, 3 more partial waves 
appear to leading order: the longitudinal s-wave amplitude 
Lq + and the p-wave amplitudes Li + and Li_, describing 
the excitation of the A(1232)) and A^*(1440), respectively. 
Moreover, all the amplitudes are functions of Q 2 , that is, 
they probe the spatial distribution of pion production on 
the nucleon. The first investigations performed at NIKHEF 



The total photoabsorption cross section ot for the proton 
is displayed in Fig. [30] It clearly exhibits 3 broad resonance 
structures on top of a strong background. These structures 
correspond, in order, to magnetic dipole (Ml) excitation 
of the A(1232) resonance, electric dipole (El) strength 
near the resonances A r *(1520) and A r *(1535), and electric 
quadrupole (E2) strength near the Af*(1680). The figure 
also shows the contributions of the most important channels: 
the one-pion channels dominate up to v w 500 MeV, the 
two-pion branching ratio becomes comparable in the second 
resonance region at v w 700 MeV, and the large 77 branch 
of the resonance A^*(1535) is hidden in the background. 
Because the nucleon resonances lie above the one-pion 
threshold, any separation in a continuous background and 
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FIG. 29 The separate d cross sections ctq, gl t, and cttt as well a s the beam helicity asymmetry A' L T at fl* = 90°. Solid line: HBChPT 
^Bernard et ali 1 1996d) . dashed line: MAID jDrechsel et~aH Il999h . dashed-dotted line: DMT model iKamalov eFqll l200lh . Figure from 
IWeis et al.\ l200H) . 



discrete resonances is necessarily model-dependent. In 
particular, because background and resonance contributions 
interfere, a careful analysis of the partial waves and their 
relative phases is mandatory. Figure [3J] depicts the partial 
cross sections for the different decay channels investigated so 
far. The rapid increase of the two-pion contribution between 
400 and 600 MeV is clearly seen, from whereon it provides 
more than half of the total absorption. Another interesting 
feature is the dominance of charged pion production both 
below and above the A(1232) resonance. Finally, one notes 
the small T) decay branch, which corresponds to s-wave T) 
production mediated by the 5n(1535). 



In the past, the primary method to unravel the nucleon 
resonance spectrum were experiments with strong interac- 
tions, in particular pion-nucleon scattering. These da ta have 
been s ystem a tically studied by the Karlsruhe group dHohlen, 
Il983t iKochl 1 19851) and the GWU Collabo ration using the 
code SAID (|Arndt et al. L 120021 120061 12004 . A summary of 
the known spectroscopic information on nu cleon resonance s 
is given by the Particle Data Group (PDG) dYao et aZ.ll2006l) . 
It is the objective of these studies to determine the relevant 
characteristics of the resonances, their pole positions, widths, 
decay channels, and branching ratios. In a first step, the 
full data base is fitted within the framework of a partial 
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FIG. 30 The total photoabsorption cross section <jt on the proton as 
function of the photon lab energy v . The various lines represent the 
MAID results jDrechsel et al ., 1999) for the total cross section (solid 
line), one-pion channels (dashed line), more-pion channels (dashed- 
dotted line), and rj c hannel (dotted line) . Full circles: total cross 
section from MAMI ( MacC ormick et all 1 19961) . open circles: data 
from Daresbury ( Armstron g et q/.LIl972h . open triangles: two-pion 
production ( Brag hieri et q/.L l l995). Figure by courtesy of J. Ahrens. 
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FIG. 31 The partial and total cross sections for the absorption of 
photons on the proton as a function of the photon lab energy v 
obtained at MAMI. Open circles: total cross section as in Fig. [30l 
iMac Cormick et all 1 19961). s olid a nd open triangles: nir + decay 
channel dMacCormick et all 1 19961) . open diamonds: p-K decay 
channel I MacCormick et all 1 19961) . asterisks: two- pion production 
llBraghieri et all. 1 19951) . open squares: prj channel l lKrusche et all 
1 19951) . Figure by courtesy of J. Ahrens. 



wave analysis. If a particular partial wave shows some rapid 
increase of the scattering phase over a limited energy region, 
the fit is then repeated with a form containing both a smooth 
background and a resonance form, mostly of the Breit-Wigner 
shape. Because of the strong decay channels and large res- 
onance widths of typically 100 MeV and more, the ideal 
resonance form is only realized for the first resonance, the 
A(1232): the pion-nucleon scattering phase #33(1^) goes 
through 90°, the real part of the multipole vanishes, and the 
imaginary part has a maximum near W — Mr w 1232 MeV. 
Because of inelastic channels, overlapping resonances, and 
energy-dependent backgrounds, these conditions are not 
fulfilled by the higher resonances. A mere "bump" in a partial 
wave is not necessarily a resonance, it may also originate 
from the opening of a new channel, which usually produces 
an asymmetric resonance shape. The "speed-plot" technique 
is particularly useful to probe the resonance structure. It 
requires the derivative of the partial wave amplitude with 
regard to the energy W, which is then compared to the 
corresponding derivative of an ideal Breit-Wigner resonance. 
In this way one determines the pole position and the residue 
of the multipole in the complex energy plane, which are 
unique characteristics of a resonance. 

Because the cross sections for electroexcitation are sup- 
pressed by the fine structure constant a em , systematic studies 
of the resonance structure were only possible after the advent 
of the new electron accelerators providing high duty-cycle 
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FIG. 32 Resonance contributions to the total photoabsorption cross 
section of the proton. The calculations are from the unitary iso- 
bar model MAID2003 torechsel etaH Il999h . Figure by courtesy 
of R. Beck and L. Tiator. 



and flux density of photons and electrons. The virtues of these 
investigations are obvious: the selectivity through transverse 
electric and magnetic as well as longitudinal photon fields 
allows for the separation of multipoles and investigations 
of the spin degrees of freedom, the study of specific final 
states with different mesons by coincidence experiments, 
and the possibility to resolve the spatial distribution of 
charge and current densities by varying the virtuality of the 
photon. Figure [32] sets the scene by showing the individual 
contributions of the main nucleon resonances to the total 
photoabsorption cross section. The resonances in the figure 
are labeled in the spectroscopic notation, see subsection III.EI 
As an example, the A(1232) = ^33(1232) has the pion in a 
p wave and both isospin and spin are 3/2. As is evident from 
Fig. [32] the A(1232) is the dominant feature in the resonance 
spectrum. It contributes more than 400 fib at the maximum of 
the absorption cross section. Already in the second resonance 
region, there are several overlapping resonances, the small 
Roper resonance Pu (1440) with the same quantum numbers 
as the nucleon, as well as the relatively strong £>i3(1520) and 
the weaker Sii(1535), both mainly excited by electric dipole 
radiation, all on top of a large background (not shown in the 
figure!) and on the tails of the neighboring resonances. In the 
third resonance region, the figure shows only the dominant 
Fi5(1680) resonance with a concentration of the electric 
quadrupole strength. The number of known (and unknown) 
resonances incre ases in the higher part of the spectrum, see 
lYao et all (2006). However, somewhere above W « 1.5 GeV 
the notion of "resonances" becomes problematic and has to be 
replaced by a continuum which, of course, can be expanded 
in partial waves as discussed in subsection lll.EI The problem 
is clearly not on the experimental side, but in the modeling of 
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the spectrum, which for W > 1.5 GeV should be based on 
coupled channel calculations. As has been me ntioned before , 
the most valid technique is the "speed plot' ' dH5MeE[l983), 
which gives information about the resonance position in the 
complex energy plane. However, the distinction of "reso- 
nances" and "background" becomes increasingly difficult, 
and looses its meaning for energies somewhere above 2 GeV. 
The "missing resonances" issue of the constituent quark 
model (CO M) that predic t s more resonances than hav e been 
observed dlsgur and Karll Il979t iMetsch et all 120031) . may 
disappear in v iew of the "resonance-background" problem 
(lThomall2005l) . 

As in the case of the elastic form factors discussed in 
section [Till the Q 2 dependence of the multipole transitions 
provides information on the spatial distribution of these 
observables. Let us discuss this issue for the A(1232). 
The 3 (real) transition form factors G* M , G* E , and Gq are 
related to the (complex) partial wave amplitudes according 
to Eq. (f30l >. The magnetic dipole transition M%+ dominates, 
the electric and Coulomb quadrupole transitions Ei + and 
5i+ are much smaller. A finite value of the quadrupole 
moment requires that the wave functions of either the nucleon 
or the A(1232), or most likely both are deformed. In the 
CQM, such a deformation follows from the tensor force 
contained in the color-hyperfine interaction among the 



containe d in the co lor-hypernne interaction among the 
quarks dCapstick and Karl Il990t |Pe Ruiula etaH Il975t 



Isgur et ai, 119821). However, typical CQM calcu lations 



dCapstick and Keisteii 1 19951: |Pe Sanctis et all l2005bl) under- 
estimate the electric and Coulomb quadrupole amplitudes 
E\+ and Si+. In models with pionic degrees of freedom, 
the deformation arises naturally from the spin-dependent 
coupling of the pion to the quarks. The pions have been 
introduced in several a pproaches , such as chiral bag models 
dBermuth et all 1 19881 iLu et all Il997[ IVento et all Il980h 
dynam ical models (IKamalov and Y ang, 1999; Sato and Lee, 
200 lb. and effective field theor i es (iGail and H emmert, 120061; 
Pascalutsa and Vanderhaeghe"nl 120051 120061) . Although 
differing considerably in the details, all these models describe 
the experimental data reasonably well. In particular the chiral 
effective field theories (EFTs) are based on a systematic 
expansion in terms of the external momenta, the pion mass, 
and the NA mass splitting. Contrary to the dynamical models 
they are gauge and Lorentz invariant, however, being based 
on a perturbative expansion, the unitarity condition is only 
approximately fulfilled. 

Of course, the deformation can not be measured for the 
nucleon ground state, because the intrinsic static quadrupole 
moment Qq and the observed quadrupole moment are related 
by Qf s = [3(M) 2 - J(J + 1)]/[J(2J - 1)] Q . This rule 
allows one to observe the quadrupole moment of the A(1232) 
but, of course, the highly unstable A(1232) can not serve 
as a static target. Therefore, the deformation issue can only 
be accessed through the electromagnetic transitions from the 
nucleon to the A(1232). In the following we discuss the 
ratios of the multipole transitions at the resonance position, 
Rem and Rsm as defined by Eq. d33l . We repeat that 



the resonance position of the A(1232) is uniquely defined 
within the validity of the Fermi-Watson theorem, that is, (I) 
the real parts of all 3 NA amplitudes vanish at W — Ma, 
because all the amplitudes carry the same phase, namely the 
pion-nucleon phase shift in the partial wave with I = 1 and 
/ = J= lS 33 (W). 



Many experimental and theoretical investigations have been 
devoted to the N A transition with real photons, for a summary 
of this work see IBeckl (120061) . The most precise value for the 
ratio of the multipoles, 



R EM = (-2.4 ±0.16 stat ± 0.24 



systj 



(108) 



was obtained with linearly polarized photons dBeck et all 
2000: LeukelLl200ll) . These authors also studied neutral and 
charged pion production on the proton, which is necessary to 
isolate the isospin § amplitude relevant for the NA transition. 
At finite Q 2 the same physics questions were addressed by 
pion electroproduction, e + p — > e' + p' + 7r°, measuring 
the scattered electron in coincidence with the recoil proton 
detected in a high-resolution magnetic spectrometer. The 
small solid angle of these instruments does not limit the ac- 
curacy since the protons are focused by the relativistic boost 
along q. Such exper iments were first performed at MIT/Bates 
( Mertz et ~all l200ll) and then extended by ISparveris et al\ 
(2005) to out-of-plane angles with the OOPS spectrometer. 
This work was continued in the framework of the Al Collab- 
oration at MAMI/Mainz j Elsner et all l2006t ISparveris et all 
120061; IStave et all 120061) . In particular, IStave et all d2006l) 



also measured the structure function a' LT by use of polarized 
electrons and out-of-plane proton detection. The experimental 
results for the A multipoles in the low-Q 2 region are com- 
pared to several model calculations in Fig. [33] We observe 
that the leading multipole M i+ is described quite well by the 
models. It is also also worthwhile mentioning that the dy- 
namic models ascribe a third of the magnetic dipole strength 
to the pion cloud. The predictions scatter much more with 
regard to the quadrupole strength as shown by the ratios Rem 
and Rsm- Concer ning the origin of "deformation", both the 
dynam ic models dKamalov and Yangl 1999 : Sato and Lee , 
| 200ll) and effective field theory dGail and Hemmed, 12006 ; 
iPascalutsa and Vanderhaegheri 120051 |2006) agree that the 
multipoles E\ + and S\+ are essentially due to the pion cloud. 
For a detailed comparison of these models in the A region, 
see lDrechselandTiatorl d2007l) . 



In many of the mentioned contributions, the authors 
have pointed out a considerable model-dependence of the 
analysis. It i s ther efore a substantial progress in this field that 
iKellv et al\ d2007l) performed a series of double-polarization 
experiments at Q 2 = 1 GeV 2 near the A region. Altogether 
they extracted 16 of the 18 independent response functions, 
most of them for the first time. As mentioned before, the 
experimental cross section is obtained by summing bilinear 
products of multipole amplitudes, Ylee> M-\M-t> . Whereas 
the unpolarized response functions or cross sections, in 
shorthand ut, <Jl, &tt, and o^t are obtained from the real 
parts of these products, many of the polarized cross sections 
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FIG. 33 The ratios EMR=R E m and CMR=R E m together w i th the amplitude M, 3 ^ 2 of the NA(1 232) transition as function of Q 2 
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are given by the imaginary parts, the first example being 
the "fifth structure function" a' LT measured with polarized 
beams. Because the NA multipoles carry the same phase, 
their product can only contribute to responses built from 
the real parts. On the other hand, the responses containing 
the imaginary parts yield information on the interference 
be tween the NA and background multipoles. The findings 
of iKellv et all d2007l) can be summarized as follows: (I) 
response functions governed by real parts are in general 
agreement with recent model calculations, (II) response func- 
tions determined by imaginary parts may differ substantially 
from the experiment and among the calculations, (III) the 
multipole analysis yields better results than a (truncated) 
Legendre series, and (IV) the model builders should go 
back to the drawing board to get better control of the non- 
resonant background, in particular at the larger virtualities Q 2 . 

As mentioned above, the typical CQM calculations un- 
derestimate the A multipoles, in particular the electric and 
longitudinal ones. Therefore, any successful description of 



these observables needs a pion cloud, at least if one insists 
on a reasonable size of the quark bag. Such descriptions are 
the chiral bag models, the dynamical models, and effective 
field theories. Taking all facts together, one is again forced 
to accept the dominant role that pions play for the structure 
of the nucleon. Of special interest is the recent work to solve 
QCD on the space-time lattice. Since these calculations can 
not yet be performed at the small (current) quark masses 
corresponding to the physical pion mass, one uses very 
large quark masses leading to pion masses > 300 MeV. 
The results are then extrapolated to the physical pion mass 
by extrapolating functions, ideally as derived from chiral 
effective field theories. Such procedure is not undisputed 
because the chiral expansion is hardly valid at pion masses 
much larger than the physical mass. However, the chiral 
extrapolation shows that unexpected phenomena occur when 
the pion mass is lowered from a few hundred MeV to its 
physical value: near the pion mass value for which the 
A resonance can decay, the chiral extrapolation becomes 
non-analytic, which leads to a kink and strong curvature in 
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FIG. 34 The helicity amplitude A 1/2 for the Sii(1535) in the 
p(e, e'p)ri re a ction. The data are represented by fi lled circles: 
Denizl i et al\ d2007h . crosses: I Armstrong et al\ J 19991). and ope n 
circles: earlier publications as specified by iDenizli et al.\ d2007h . 
The predictions are based on qu ark models and given by the solid 
line dCapstick and Keisterl 1 1995 ) and t he dotted line I Aiel lo et ail 
Il998h . Figure from lDenizli et all J2007h . 



the extrapolation formula (Pascalutsa et al. , 2007). 

As a further instructive example we present some results 
for the electroproduction of r\ mesons, which are mainly 
produced by the decay of the 5n(1535) resonance with 
J = I = §, t = 0. As is evident from Fig. [3U the 5u(1535) 
is buried under the total cross section. However it is clearly 
seen in the -q channel, because this resonance has an 7/ branch- 
ing ratio of 45-60% com pared to a few per cent for other 
excitations of the nucleon dVrana et aZl boOO). Eta photopro- 
duction experiments at threshold show a strong increase of 
the cross section in the range of the resonance and an s-wave 
angular distribution i n agreement with the given assignment 
dKrusche et all 1 1995b . However, the shape of the resonance 
is very asymmetric, and also t he speed-plot analysis does 
not yield satisfactory solutions dHohletl 1 19931) . Because of 
this unusual behavior several alternative interpretations have 
been given, for exa mple in terms of a KS molecular state 
dKaiser et al.l\l99% . Also the helicity amplitude A1/2 of this 
resonance has an unusually soft form factor, as displayed in 
Fig. [34] for the jp(e, e'p)j ] reaction measured at the Jefferson 
Lab dDenizli et a/.ll2007l) . It is striking how flat this transition 
form factor stays compared to the typical dipole form for 
other form factors of the nucleon. We further notice that the 
constituent quark model calculations, relativistic or not, can 
not fully explain the slow falloff with Q 2 . It is also tempting 
to identify the structure near Q 2 fa 0.2 GeV 2 with a meson 
cloud effect as discussed in section lHI.AI 

Another small but interesting resonance in the second 
resonance region is the Roper resonance Pi 1 (1440) with the 
same quantum numbers as the nucleon. The CQM describes 
this resonance by a radial excitation of the nucleon, and 
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FIG. 35 The helicity amplitudes A 1/2 and S 1/2 of the Pn(1440) 
as function of Q 2 . The triangles are from the analysis of 
lAznaurvan et al\ d2005l) based on data of the CLAS Collaboration 
at JLab and obtained from dispersion relations (open triangles) and 
the unitary isobar model (filled triangles ) . The solid line is the global 
solution of MAID2007 dDrechsel etali |2007|) . a fit to the full data 
basis, and the filled circles are local fits to the data in th e particular 
energ y region. The d ata point at Q 2 = is from the PDG dYao et all 
l2006h . Figures from dDrechsel efflZll2007l) . 



therefore it should be sensitive to the radial form of the bag 
potential. However, its ma ss is mu c h lowe r than expected 
in simple quark models. iLi et all d 19921) suggested that 
the Roper could also be a quark-gluon hybrid state, which, 
however, can not be excited by the longitudinal current or 
Coulomb field. Recent data on the electroproduction of 
this resonance are shown in Fig. [35] and compared with the 
single-energy and global solutions of MAID2007. We note 
that the amplitude Aim (magnetic dipole transition) has a 
zero crossing at small momentum transfer. Furthermore, the 
amplitude S1/2 (Coulomb monopole transition) rises to quite 
large values, which rules out the quark-gluon hybrid model. 

In section [VI] we also discuss the helicity structure of the 
second and third resonance regions with regard to several sum 
rules for real and virtual photons. However, with increasing 
excitation energy it becomes more and more difficult to iso- 
late individual resonances by inclusive cross sections. It will 
therefore take a full-fledged program with polarized beams 
and targets as well as recoil polarization to analyze the higher 
mass region and to identify resonance structures on top of the 
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VI. SUM RULES 

The persisting problem of gaining a quantitative under- 
standing of nucleon structure is one of the reasons why we 
are interested in sum rules. Being based on quite general prin- 
ciples like causality, unitarity, Lorentz and gauge invariance, 
sum rules should be valid for every model or theory respect- 
ing these principles and having a "reasonable" high-energy 
behavior. Therefore, the agreement or disagreement between 
theoretical predictions and the sum rule provides invaluable 
information on the quality of the approximations involved and 
whether or not the relevant degrees of freedom have been in- 
cluded. Specifically, if we compare a sum rule value with 
accurate experimental data up to a certain maximum energy, 
we learn whether the physics responsible for the sum rule 
is provided by the phenomena up to that energy, or whether 
possibly new degrees of freedom come into the game. As 
an example, Baldin's sum rule relates the forward scalar po- 
larizability to an energy-weighted integral over the total pho- 
toabsorption cross section and thus allows for an independent 
check of the results from Compton scattering. Moreover, re- 
cent double-polarization experiments have determined the he- 
licity structure of this cross section for the proton. From these 
data, the forward spin polarizability has been obtained, and 
the Gerasimov-Drell-Hearn (GDH) sum rule has been verified 
within the experimental error bars of less than 10 %. Data 
have also been taken for the neutron and are under evalua- 
tion. With some caveat in mind, these sum rules can be gen- 
eralized to the scattering of virtual photons. The integrands 
of the respective integrals are related to the electroexcitation 
cross sections and, in the limit of deep inelastic scattering 
(DIS), to the nucleon structure functions. The resulting gen- 
eralized integrals and polarizabilities depend on the photon's 
4-momentum and therefore contain information on the spatial 
distribution of these observables. 



A. Sum rules for real photons 

1 . Forward dispersion relations and sum rules 

The forward scattering amplitudes / and g of Eqs. doTl ) - 
d63l can be determined by scattering circularly polarized pho- 
tons (helicity A = ±1) off nucleon targets that are polar- 
ized along or opposite to the photon momentum q as shown 
schematically in Fig. [36] If the spins are parallel, the he- 
licity of the intermediate hadronic state takes the value |. 
Since this requires a total spin J > |, the transition can 
only take place on a correlated 3-quark system. For oppo- 
site spins, on the other hand, the helicity is conserved and the 
scattering can also take place on an individual quark. Denot- 
ing the Compton scattering amplitudes for these two exper- 
iments by T 3 / 2 and T 1/2 , we find f(v) = §(Ti/ 2 + T 3 / 2 ) 
and g(v) — j(2i/ 2 ~ ^3/2)- Furthermore, the total absorp- 
tion cross section is given by the spin average over the helicity 



FIG. 36 Spins and helicities in the reaction 7 + TV — > N * . The open 
arrows denote the projections of the spin on the photon momentum, 
S z , the helicities h are the projections on the respective particle mo- 
mentum, and the photon has right-handed helicity, A = +1. 
Left: N(S Z = 1/2, h = -1/2) -* N*(S Z = h = 3/2), 
right: N(S Z = -1/2, h = 1/2) -> N*{S Z = h = 1/2). 



cross sections, 



°T = 2 ( a l/2 + C3/2) 



(109) 



and the helicity-dependent cross section by the helicity differ- 
ence, 



(7 J 1 J 1 



1 , 

2 (ai/2 



C3/2J 



(HO) 



Based on unitarity and causality, the optical theorem relates 
the absorption cross section to the imaginary parts of the re- 
spective forward scattering amplitudes, 



Im f(u) 



8tt 



(0-1/2 (f) + 03/2M) = j- &t{v) , 



lmg{v) = — (oi/ 2 (i/) - 0-3/2 (v)) = —(ttt{v). (Ill) 

Because of the smallness of the fine structure constant a em , 
all higher order electromagnetic processes are at the per cent 
level. We may therefore neglect the absorption below the 
threshold for pion production, vq = m^l + m T /2M) w 
150 MeV, thereby assuming that the scattering amplitude is 
real in this region. In the next step, one has to study the high- 
energy behavior of the absorption cross sections. As predicted 
by Regge theory and also seen by the data, the total absorption 
cross section increases at the highest energies reached by the 
experiment. Although this increase is not expected to continue 
forever, we can not expect that the unsubtracted dispersion in- 
tegral converges, and therefore we subtract f(v) at v = and 
identify /(0) with the classical Thomson amplitude. By use 
of the crossing relation and the optical theorem, the subtracted 
dispersion relation takes the form 



/(") 



v 2 r 



o T {v') 



dv' 



(112) 



For the odd function g{y) we assume the existence of an un- 
subtracted dispersion relation, 



qi/ 2 (*/)-g 3 / 2 (z/) 

— ; v dv . (113) 



If these dispersion integrals exist, they can be expanded as a 
Taylor series in v 2 which converges for \ v\ < Uq. Comparing 
these power series to the low en ergy theorems of Eqs. d62l and 
(1631 . we obtain the sum rule of lBaldinl(ll960t) . 
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(114) 
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the su m rule o f|G erasimo vl dl965Lll966h and lDrell and Heard 
(1966), 



2M 2 



dV — iGDH , 



(115) 



and the forward spin polarizability (Gell-Ma nn et al. i tun, 

1 r cr 3/2 (i/) - ff 1/2 (i/) , , 



4tt 2 



■dv' . (116) 



2. Photoabsorption cross sections for the proton 

The total photoabsorption cross section ctt in the resonance 
region of the proton is shown by Fig. [30] This figure displays 
3 resonance peaks on top of a large background. Above 
the resonance region, ot is slowly decreasing towards a 
minimum of about 115 [ib at W « 10 GeV. At the highest 
energies, W « 200 GeV (corresponding to v ~ 2 • 10 4 GeV), 
the experiments sho w a slow increa s e with energy of th e 
form cr T - W - 2 (lAidefa/l Il995t iDerrick et ali 11994 . 
in accordance with Regge para metrizations thro ugh a soft 
pomeron exchange mechanism (Cudell et al. , 2000). Given 
this information, the rhs of Eq. (II 141) can be constructed, with 
the most recent numerical result given by Eq. d64l ). In this 
way Baldin's sum rule provides a rather precise value for 
the sum of the 2 scalar polarizabilities, which serves as an 
important constraint for the analysis of Compton scattering. 

During the past years also the helicity difference <jtt has 
been measured. The pioneering experiment was carried out 
by the GDH Collaboration at MAMI for photon en ergies 
between 200 and 800 MeV ( Ahre ns et al .1 12QQ0L 1200 ll) . and 
then extended into the e nergy range up to 3 GeV at ELSA 
(lDutzefadl200ll2005h . These data allow us to verify the 
GDH sum rule for the proton within an accuracy of less than 
10 %. Because the integral for the forward spin polarizability 
converges much better, it is essentially saturated by the 
MAMI data at 800 MeV. As shown in Fig. [37] the helicity 
difference fluctuates much more strongly than the total cross 
section <tt- The threshold region is dominated by s-wave 
pion production, i.e., intermediate states with spin | that can 
only contribute to the cross section Ci/2- In the region of 
the A(1232) with spin J = |, both helicity cross sections 
contribute, but since the transition is essentially Ml, we find 
the ratio 03/2/ '01/2 ~ 3, and therefore the helicity difference 
becomes large and positive. The figure also shows that (73/2 
dominates in the second and third resonance regions. It was 
in fact an early success of the quark model to understand this 
feature by a cancelation of the spin and convection currents 
for <7i/2. The data at the higher energies indicate a fourth 
resonance region (1800 Me V< W < 2000 MeV) followed by 
a continuing decrease of Act with a cross-over to negative val- 
ues at v > 2.0 GeV, as predicte d by an extrapolation of data 
from deep inelastic scattering dBianchi and Thomasl 1 19991; 
ISimula etal. , 2002). At high v, above the resonance region, 
one usually invokes Reg ge phenomenology t o argue that the 
integral converges dBass and Brisudovalll999h . In particular, 
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FIG. 37 The helicity difference Act = a 3 / 2 — °"i/2 f° r the pro- 
ton as functio n of the photon energy v. The experimental data are 
from MAMI dAhrens et all |2000| . l200lh (full circles) and ELSA 
dDutz et all 12001 l2005h" (open c i rcles) . The various lines repre- 
sent MAID results dDrechsel et a/.LI 19990 for the total helicity differ- 
ence (solid line), one-pion channels (dashed line), more-pion chan- 
nels (dashed-dotted line), and rj channel (dotted line). Figure from 
iDrechsel and Tiatorl ( 120041) . 



one obtains for the isovector channel u\i% — (J3/2 — * ^ ai ' _1 
at large v, with —0.5 < ay < being the intercept of 
the ai(1260) meson Regge trajectory. For the isoscalar 
channel, Regge theory predicts a similar energy behavior 
with as — —0.5, which is the intercept of the isoscalar 
/i(1285) and /i(1420) Regge trajectories. However, these 
ideas have still to be tested experimentally. We observe that 
the large background of non-resonant photoproduction in <jt 
(Fig. [3Qb has almost disappeared in the helicity difference 
Act (Fig. [37] ). i.e., the background is "helicity blind". As a 
result the two helicity cross sections for real photons remain 
large and nearly equal up to the highest energies, at values 
of <ti/2 ~ CT3/2 ~ 120 /ib. We conclude that the real photon 
is essentially absorbed by coherent processes, which require 
interactions among the constituents such as gluon exchange 
between two quarks. This behavior differs from DIS, which 
refers to incoherent scattering off the constituents. As a 
consequence the ratio a'3/2/0'1/2 tends to zero with increasing 
virtuality Q 2 , because the absorption on an individual quark 
leads only to final states with helicity i. 

As shown in Fig. [38] the GDH Collaboration has also mea- 
sured the helicity difference for the deuteron dAhrens et ali 
l2006h . The upper panel of this figure displays the total cross 
section, which yields a considerably smaller resonance peak 
than predicted fo r free nucleons (MAID) and also by the dy- 
namical model of lArenhovel et al. (2004). However, the mea- 
sured helicity difference (lower panel) agrees quite well with 
both models. The experiment has b een continued t o energies 
up to W = 1.8 GeV at ELSA dDutz et ali 120051) . yielding 
positive values of typically 50 /ib above the resonance re- 
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FIG. 38 The total photoabsorption cross section or (upper panel) 
and the helicity difference Act = cr 3 / 2 — <?\i% ( lower panel) 
for th e deuteron. The experimental data are from lAhrens et al\ 
j2006l) (full circ l es). iMacCormick et al\ d 19971) (open circles), and 
lArmstrong et al\ dl972l) (asterisks). The theoretical predicti ons are 
from MAID03 (d ashed lines) and lArenhovel et al\ d2004l) (solid 
lines). Figure from Ahrens et al. ( 20Qg)^ 



gion. This contrasts the proton result, which turns out neg- 
ative in this region. Both experimental findings agree with 
the prediction of Regge theory that the asymptotic tail of Act 
should be positive for the neutron and negative f or the proton 
dBianchi and Thomasl [l999t ISimula et azll2002l) . 



3. The GDH sum rule 



The GDH sum rule, Eq. dl 151 ). relates the anomalous 
magnetic moment (amm) of a particle to an energy-weighted 
integral over the helicity-dependent photoabsorption cross 
sections. This relation bears out that a finite value of the 
amm requires the existence of an excitation spectrum, and 
that both phenomena are different aspects of a particle with 
intrinsic degrees of freedom. A further property of composite 
objects is their spatial extension in terms of size and shape, 
which reveals itself through the form factors measured by 
elastic lepton scattering. I n conclusio n , the discovery of 
the proton's large amm by IStern et al. I (119331) marked the 
beginning of hadronic physics. The experiment indicated 
that the proton was a microcosm in itself, and in this sense 
the findings of Stern and collaborators were revolutionary. 
Today, seventy-five years later, we are still struggling to 
describe the structure of the strongly interacting particles 
in a quantitative way. At this point the reader may well 
ask why should the GDH sum rule exist and what do we 
learn from it. In fact it was pointed out many years ago that 
the GDH sum rule holds at leading order in perturbation 



theory for the stand ard model of electroweak interactions 
(lAltarelli ef aZ.LIl972l) . Later on this result was generalized to 
any 2 — > 2 process in supersy mmetric and other field theories 
dBrodskv and Schmidn, 1 1995b . The essential criterion is that 
these theories start from point-like particles, and then the 
GDH sum rule should hold order by order in the coupling 
constant. As an exam ple, the GDH sum rule has been proven 
in QED up to C(e 6 ) dDicus and Vegal |2001|) . In passing we 
note that the GDH sum ru le was also inves tigated in quantum 
gravity to one-loop order ( Goldbergl |2000|) . The result was a 
violation of the sum rule, which may however be due to our 
ignorance of quantum gravity in the strong coupling (high 
energy) limit. 

The amm of a particle is defined by the following relation 
between the total magnetic moment fl and the spin S: 



A I 



(Q + k)S : 



(117) 



with eQ the charge and M the mass of the particle. We also 
recall that the ratio between the magnetic moment jl and the 
orbital angular momentum L of a uniformly charged rotating 
body is eQ/2M, whereas Eq. (Ill yields eQ/M as ratio 
between the "normal" magnetic moment fl and the spin S, 
because of the gyromagnetic ratio g=2 predicted by Dirac's 
equation for a sp in h particle. Contrary to the conjecture by 
iBelinfantd dl953l) that g=l/S, the "natural" value of the gy- 
romagnetic ratio is g=2 for every point particle, independent 
of its spin. This is necessary if one insists on a well-behaved 
high-energ y scattering amplitude and a reliable perturbative 
expansion ( Ferrara et aUll992tlWeinbera,ll970l) . in the sense 
that any deviation from this value must be related to finite 
size effects. Such spatially extended phenomena, however, 
do not affec t the high-energy limi t of Compton scattering. 
In particular Brodsky and Pri mackl d 19691) verified the GDH 
sum rule for a composite system of any spin on the basis of 
the spin-dependent interaction currents associated with the 
cm motion. To further illustrate this point let us consider the 
small amm of the electron, which can be evaluated in QED 
to 10 decimal places. Since the associated photon-electron 
loops are spread over a spatial volume of about 10 6 fm 3 , a 
high-energy photon of a few hundred MeV or a wavelength 
A < 1 fm will decouple from such a large volume. Therefore 
the amm does not affect the high-energy limit and, along the 
same lines, the amm of the electron does not keep us back 
from using the electron as a n ideal point particle to study the 
form factors of the nucleon dDrechsel a nd Bermut5 ll99lh . 



The lhs of the GDH sum rule, Eq. dTT5l yields i£ DH = 
205 fib and Iq D h — 233 fib for proton and neutron, respec- 
tively. However, the first estimates based on the then existing 
photoproduction data led to 261 fib for the proton and 183 fib 
for the neutron dKarrinen,ll973l) . Over the following years the 
predictions moved even further away from the sum rule values 
in spite of an improving data basis, simply because these data 
were not sensitive to the helicity difference of the inclusive 
cross sections. Many explanations for an apparent violation of 
the sum rule followed, but in view of the new experimental ev- 
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idence we may safely discard these ideas. TablelCTsummarizes 
our present knowledge on the GDH integral and the forward 
spin polarizability of the proton. The threshold contribution 
for v < 0.2 GeV is evalu ated by the MAID mu ltipole analysis 
of pion photoproduction dDrechsel et al.l [l999l) . with an error 
estim ated by comparing to the SAID analysis dArndt et al.l 
|2002|) . The resonance region up to v = 2.9 GeV is deter- 



mined by th e experimental data taken a t MAMI ( Ahren s et al.l 
l2000i 1200 ll) and ELSA dDutz et al.l I2003L [2005b, and the 
asymp totic contribution is based on the Regge ana lysis of 
DIS dBianchi and Thomasl fl999l ISimula et all 120021) . Sum- 
ming up these contributions, the GDH sum rule value is ob- 
tained within the given error bars. Because of the different 
energy weighting, the forward spin polarizability converges 
much better and is therefore completely determined by the ex- 
isting data. Because the helicity dependent cross section <ttt 



energy [GeV] 


7 GDH [Mb] 


7o P [ID" 4 fml 


<0.2 


-28.5 ±2 


0.95 ±0.05 


0.2-0.8 


226 ± 5 ± 12 


-1.87 ±0.08 ±0.10 


0.8-2.9 


27.5 ±2.0 ±1.2 


-0.03 


>2.9 


-14 ± 2 


±0.01 


total 


211 ± 15 


-0.94 ±0.15 


sum rule 


204 



TABLE II The contribution of various energy regions to the GDH 
integral Iq DH and the forward spin polarizability 7q of the proton 
(see text for explanation). 

is strongly energy dependent, the GDH integral is very sen- 
sitive to experimental errors. It is therefore quite satisfying 
that all the decay channels were separately identified in the 
range of 200 MeV < v < 800 MeV. The one-pion channel 
opens at w 150 MeV and dominates the cross section up 
to v rj 500 MeV, except for small contributions due to radia- 
tive decay and the onset of two-pion production. The helicity- 
dependent cross section for t he one-pion channel has t he fol- 
lowing multipole expansion dDrechsel and Tiatorlll992l) : 

Ott = 4t3 j|£ 0+ | 2 -|M 1+ | 2 ±6Re(£ 1+ Ai 1+ ) 



-3|£xh 
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-6Re(£2_M 2 _) ±3|M 2 -| 2 ± 



As we have seen, the threshold region was not covered by the 
experiment. The dominant multipole in this region is E$ + , 
which corresponds to an electric dipole (El) transition lead- 
ing to the production of (mostly charged) pions in an s wave. 
This multipole is well described by pio n photoproduction dat a 
at threshold (see s ection [VAl. ChPT dBernard et ali Il991al) . 
dispersion th eory dHanstein et a/ill 9971). and phenome nolog- 
ical analysis dArndt et al.l 120021: iDrechsel et al.l fl999). With 
increasing photon energy, the first resonance becomes more 
and more dominant, mainly because of the magnetic dipole 
transition to the A(1232). Although the associated elec- 
tric quadrupole transition is strongly suppressed by the ratio 
Rem — Ex+/Mx+, the GDH experimen t permits an inde- 
pendent measurement dAhrens et a/i 120041) . because the prod- 



uct El + Mi + appears with a factor 6 in Eq. ( II 18b . Alto- 
gether the MAMI data are in good agreem ent with the multi- 
pole analysis in the fi rst resonance region dArndt et al l l200l 
IDrechsel ef aZl fl999). However, even relatively small effects 
count, because this region provides the lion's share to the sum 
rule, that is about 175 fib between v = 250 and 450 MeV. 
In a similar way the N — > -D13 (1520) transition was studied 
dAhrens et al.l 120021) . The multipoles E 2 - and M 2 -, El and 
M2 transitions, respectively, are related to the helicity ampli- 
tudes of this resonance as follows: 



A 1520 P 
A l/2 &2 



-3M 2 



A 1520 



V3{E 2 



-Mo 
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The new data yield A 



1/2 



-38 ± 3 and A 



3/2 



147 ±10, to 



be compared with the listing of the PDG, — 24± 9 and 166 ± 5 
dYao et ali [2006b . all in units of 10" 3 GeV~ 1/2 . The given 
examples demonstrate that double-polarization experiments 
provide a very sensitive tool to study resonance properties. 

Although the threshold for two-pion production lies already 
in the A region, these channels become important only for 
v > 500 MeV. The channels nn + - K°, ptt + tt~, and pir° ir° 
were separately analyzed at MAMI dAhrens et al.l l2003albh . 
As an example, Fig. [39] shows the cross sections for the 
reaction 73?^ n7r + 7r°, which also exhibits a clear dominance 
of cr 3 /2 over <J\j2- The interesting and previously unexpected 
feature is the peaking of the respective cross section at 
v k, 700 MeV or W « 1480 MeV, definitely below the 
positions of the D13 (1520) and S\i (1535) resonances. This 
proves that two-pion production can not be simply explained 
by a resonance driven mechanism as was assumed in the 
earlier estimates for the sum rule. The rj channel provides 
another interesting contribution to the sum rule. This channel 
is dominated by the resonance S%i (1535), which has an 
exceptionally large branching ratio of about 50 % for 77 decay. 
Because this resonance has I = 0, it only contr ibutes to the 
helicity cross section er^ dAhrens et al.l l200"3bb . For further 
information on the helicity structure of the different channels 
we ref er the reader to the review bv lKrusche and Schadmandl 
(120031) . The work of the GDH Collaborati on at ELSA an d 
MAMI has recently been summarized by iHelbind d2006l) . 
Several other ongoing activities or proposals for future 
experiments have been reported at conferences by the LEGS 
Collaboration at BNL, several groups at the GRAAL and 
Spring-8 facilities, the CLAS Collaboration at JLab, and 
the E-159 collaboration at SLAC. Such further experiments 
will be invaluable as independent tests of the sum rules, 
in particular at the higher energies in order to probe the 
soft Regge physics in the spin-dependent forward Compton 
amplitude. 

Table [III] shows the amm for the nucleons and standard 
"neutron targets". The most striking observation is the tiny 
value of k for the deuteron, a loosely bound proton-neutron 
system, which has isospin 1 = and spin 5=1 and is 
essentially in a relative s state. The interplay between nuclear 
and subnuclear degrees of freedom in the deutero n has been 
studied in great detail by lArenhovel et al. I d2004l) . It turns 
out that the most important nuclear channel is deuteron 
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FIG. 39 The total cross section ctt and the helicity difference 
Act = ct 3 /2 — for the reaction 7 p — > nn + -K° . The theoreti- 
cal pr edictions are shown by solid lines dHirata e[ ai, 2003), dashed 
lines jNacher and Oseti|2002|), and dotted li nes jHolvoeAl200lh . The 
data are from MAMI I IAhrens et a/.ll2003al lbT). Figure completed by 
Lothar Tiator. 
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TABLE III The magnetic moment fi in units of the nuclear mag- 
neton /ijv, the amm k, and the GDH sum rule Jgdh in units of fib 
for electrons, protons, neutrons, deuterons, and s He nuclei as well as 
fully polarized hydrogen atoms and butanol molecules. 



disintegration, 7 + d — > p + n, which yields a maximum 
value of <jp — a a ~ —1800 fib at v w 2.3 MeV. This huge 
helicity asymmetry is due to the Ml transition 3 Si — ► So, 
which changes the magnetic moments of proton and neutron 
from parallel to antiparallel. Due to the energy denominator 
v in the GDH integral, it is precisely the small excitation 
energy of the weakly bound deuteron that provides this large 
negative contribution to the GDH integral. To the contrary, 
the N — > A transition aligns the quark spins and peaks at 
v w 330 MeV with a maximum value of op — a a ~ 1100 fib 
for free nucleons. As a result, the large negative contribution 
of deuteron break-up is canceled by large positive contribu- 
tions of the subnuclear degrees of freedom, and this happens 
to three decimal places. The other neutron target is 3 He, a 
system of two protons with spins paired off and an "active" 
neutron, essentially again in s states of relative motion. As a 
result we find fi 3Ho w fi n < 0, whereas Eq. ( II 171 ) predicts a 
"normal" moment of |/xjv- Therefore, the amm of 3 He has a 
large negative value, which leads to a large and positive GDH 
integral. 

The recent experiments to determine the GDH sum rule for 
the neutron, Iq D h> nave been performe d with a froz e n- spin 
deuterated butanol (C4D9OD) target (lAhrens et all 120061) 
and a frozen-spin 6 LiD target dDutz et al. L 120051) . Setting 
aside the problems on the molecular and atomic levels, these 
experiments provide polarized deuterons whose helicity- 
dependent response has been measured in the energy region 
200 MeV< v < 800 MeV. As is obvious from the above 



discussion, a quantitative extraction of Iq DH is necessarily 
model-depe ndent. Even the prese ntly most elaborate cal- 
culation of lArenhovel et all d2004l) misses the GDH sum 
rule for the deuteron by nearly 30 fib and overestimates 
the total photoabsorption in the region of the A(1232). It 
is therefore mandatory to also measure the different decay 
channels in order to constrain the theoretical analysis. Such 
experiments are in progress, and the first results for th e 
one-pion channels are already available dAhrens et all\200(h . 
To lowest order we may assume that ^gdh + ^gdh should be 
given by the GDH integral for the deuteron, if extended from 
pion threshold to infinity. On the basis of the present analysis 
we may conclude that more than 60% of this contribution 
is due to 7T° production and another third from two-pion 
channels with at least one charged pion. Furthermore, 
Fig. [38] shows a reasonable agreement between the data 
and the MAID model in the A(1232) resonance, in which 
region the model yields quite similar values for proton and 
neutron. The additional sum rule strength for the neutron 
should therefore come from energ ies above t he A (1232). 
And indeed, the ELSA experiment dDutz et all [2005) yields 
a contribution of about 34 fib between 815 and 1825 MeV, 
contrary to earlier estimates. Furthermore, the integrand 
remains positive at the highest energies, which could indicate 
a further contribution of about 40/ib on the basis of Regge 
models. In conclusion, the present experiments confirm 
the GDH sum rules for proton and neutron, however with 
a very large systematical error in the latter case. It is 
also likely that the isovector combination Iq DH — ^gdh turns 
out negative, as required by the sum rule prediction of — 28 fib. 



The interplay of nuclear and subnuclear degrees of free- 
dom is based on general principles like low-energy theorems 
and dispersion relations. A complete answer to the remain- 
ing questions calls for experiments covering both the nuclear 
and the subnuclear energy range. It is therefore very promis- 
ing that programs are being developed for energies between 
nuclear breakup and pion threshold at TUNL/HI7S (Duke) 
dWelleil 12003). both for the nuclear physics aspects by them- 
selves and as a test of many-body calculations that are in- 
evitably required for further studies of the GDH sum rule of 
the neutron. We conclude this subsection with a possibly aca- 
demic but still interesting question. Although we have talked 
about nucleons and nuclei as targets, we could just as well 
think about projects to measure the GDH sum rule for atoms 
or molecules. Table Hill shows the amm k and the GDH sum 
rule I for several such systems. The comparison of the differ- 
ent hierarchies is quite amusing. If one tried to reconstruct, 
for example, the amm of the hydrogen atom by a GDH in- 
tegral over its atomic spectrum in the eV or keV region, one 
would find physics "beyond": the physics of e + e~ pair pro- 
duction at the MeV scale (QED) and hadronic physics above 
pion threshold (QCD). 
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FIG. 40 The two-photon exchange diagram. The blob represents the 
possible intermediate states of the hadronic system. 



B. Sum rules for virtual photons 

Doubly-virtual Compton scattering (VVCS) offers a useful 
framework to stu dy generalized GDH integrals and polariz- 
abilities This process is based on the idea that an 

incident virtual photon with definite energy v and virtuality 
Q 2 > (space-like) hits a nucleon and excites this hadronic 
system, which eventually decays into a nucleon and an outgo- 
ing virtual photon with the same value of Q 2 . Although this 
reaction can not be realized experimentally, it is not merely 
a theoretical construct. As pointed out many years ago, the 
imaginary part of this amplitude can be obtained by elastic 
scatteri ng of transversely polar ized electrons off unpolarized 
targets (De Rui ula et a^.L 1 197 lb . The asymmetry with regard 
to changing the transverse polarization is parity conserving 
but time-reversal odd and therefore vanishes in the one-photon 
exchange approximation. It appears only at sub-leading order 
as the product of the (real) Born amplitude for one-photon ex- 
change and the imaginary part of the two-photon exchange 
diagram shown by Fig. [40] This imaginary part is related 
to the VVCS tensor, however the experiment can only de- 
termine a weighted integral over a range of virtualities Q 2 , 
whereas the VVCS tensor refers to a fixed value of Q 2 . A pi- 
oneering experiment t o measure this asy mmetry has been per- 



formed at MIT/Bates (Wells et al 



now available from MAMI (Maas, 



200 1 ), and more data ar e 
20051: iMaas eJali 120051) . 



These measurements show that even for moderate virtualities, 
Q 2 0.1 GeV 2 , the asymmetry is dominated by the excited 
states of the nucleon, notably the A(1232) resonance. The 
same physical effect can also be observed for other single-spin 
asymmetries, such as the transverse target asymmetry or the 
normal recoil polarization for an unpolarized electron beam. 



1 . VVCS and nucleon structure functions 



The absorption of a virtual photon on a nucleon N is de- 
scribed by inclusive electroproduction, e + N — > e' + X. The 
cross section for this reaction takes the form (iDrechsel et al.l 



l200Tl) . 



da 
dQ,2d,€2 
a 



T v a(v, Q 2 



(120) 



= cr T + ea L - hP x y / 2e(l - e)a LT 



hP z VT 



e 2 <jtt ■ 
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Comparing this form with Eqs. ( f2Tb and d28| ), we find that the 
transverse (cry) and longitudinal (ctl) absorption cross sec- 
tions are obtained by integration of the respective differential 
cross sections over the angles of the emitted pion, and in gen- 
eral by integrating over all angles and energies of the produced 
particles. The two (inclusive) spin-flip cross sections, <jtt 
and o can only be measured with polarized electrons (helic- 
ity h — ±1) and target polarization in the direction of the vir- 
tual photon momentum (P z ) and perpendicular to that direc- 
tion in the scattering plane (P x ), respectively, or equivalent re- 
coil polarizations. In order to avoid a possible misunderstand- 
ing, we note that these inclusive spin-flip cross sections are 
not related to the differential cross sections darr /<if2* and 
d<JLT I d£l%. of Eq. d28l ) but obtained from double-polarization 
cross sections not shown in that equation. The following dis- 
cussion concentrates on the spin-flip cross sections, which are 
related to the spin-dependent n ucleon structure functions g\ 
and gi by dDre chsel et aZ.ll2003l) 
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with 7 — Q/v and K the photon equivalent energy defined 
by Eq. (|29l l. The nucleon structure functions are usually ex- 
pressed as functions of the Bjorken variable, x = Q 2 / (2Mv), 
i-e., ffi.,2 = 9i,2(x, Q 2 ). 

The VVCS amplitude for forward scattering of virtual pho- 
tons generalizes Eq. doTt by introducing an additional longi- 
tudinal polarization vector q, 



T(v, Q 2 , 9 = 0) 



s'*-ef T (u, Q 2 ) + f L (v, Q 2 ) 
io -{s 1 * xe)g TT {v, Q 2 ) (124) 
i(e'* -£*) ■ (a x q)g LT {v, Q 2 ) . 



Because the total amplitude is crossing-even, the amplitude 
9tt{v) is an odd function of v whereas gvriii) is even. In 
order to set up dispersion relations, we have to construct the 
imaginary parts of the amplitudes with contributions from 
both elastic and inelastic scattering. The elastic contributions 
are obtained from the direct and crossed Born diagrams with 
nucleons in the intermediate states and expressed by the nu- 
cleon form factors of section [TLBl 
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with v B = Q 2 /2M. We note that the amplitudes of Eq. ( 1125b 
have been split in a real contribution and a complex term con- 
taining the nucleon poles at v — ±1/3 =p ie. The inelastic 
contributions are regular functions in the complex i^-plane ex- 
cept for cuts from —00 to —vq and to +00. The optical 
theorem relates the inelastic contributions to the partial cross 
sections of inclusive electroproduction, 

\mg TT {v,Q 2 ) = ^a TT { V ,Q 2 ), (126) 

47T 

lmg LT (v,Q 2 ) = ^-a LT {v,Q 2 ), (127) 

where the products K ijtt /lt &e independent of the choice 
of K, because they are directly proportional to the measured 
cross section. Comparing the above equations with the results 
for real photons in section [VI. A| we find distinct differences. 
In particular, the transition from real to virtual photons is not 
straightforward, becau se the limits Q 2 — > and v — » can 
not be interchanged djil ll993l) . This happens because the Born 
amplitudes of VVCS have poles, which also provide imagi- 
nary contributions to the amplitudes. In more physical terms, 
the crucial difference between the (space-like) virtual and the 
real photon is that the former can be absorbed by a charged 
particle whereas the latter can only be absorbed at zero fre- 
quency, v = 0. At this point, however, the real photon ampli- 
tudes can be expanded in a Taylor series whose leading terms 
are determined by the Born terms. In particular, g?p T (i/, 0) re- 
produces exactly the leading term of the spin-flip amplitude 
g{v) of Eq. d63l l. The GDH sum rule is then obtained by 
equating this power series to an expansion of the dispersion 
integral over the imaginary parts from the inelastic processes 
for v > vq. The virtual photon case differs in two aspects: 
The imaginary parts stem from both elastic and inelastic pro- 
cesses, and the real parts of the amplitudes have two poles in 
the v plane. 



2. Dispersion relations and sum rules 

We next turn to the sum rules for the spin dependent VVCS 
amplitudes. Subtracting the pole terms from the full ampli- 
tude, we obtain the following dispersion relation: 



where 



Re g%*(u, Q 2 ) = Re ( g TT (u, Q 2 ) - Q 2 ) ) (128) 
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Because this amplitude is regular below the first threshold, 
v = Vq, it can be expanded in a Taylor series at v = 0. The 
result is 



disp/ /->2\ 2a em 2 
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M 2 f°° K(v,Q 2 )a TT (v,Q 2 ) 



di/,(130) 



(n 2, 1 r K{ V ,Q 2 )aTT(v,Q 2 ) , 
1tt(Q ) = j v -4 ^ -(131) 

Comparing these expressions with Eq. d63l ). we find Itt (0) = 
— k 2 /4 and 7tt(0) = 70- The corresponding equations for 
the crossing-even amplitude g^x are 



Ilt(Q 2 ) = —J 
Tre 2 J u , 

5lt{Q 2 ) 
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di/,(132) 
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Both functions are finite in the real photon limit, because the 
factor oltIQ in the integrand is finite for Q 2 — * 0. Further- 
more, it follows from Eq. ( 1123b that ult can be replaced by 
gi + 52, the sum of the spin-dependent structure fun ctions. 
In the limit of large Q 2 , [Wandzura and Wilczekl (1 1977b have 
shown that g\ + 92 can be expressed in terms of the twist-2 
spin structure function gi if the dynamical (twist-3) quark- 
gluon correlations are neglected, 



gi(x,Q 2 ) +g 2 (x,Q 2 



l dy 9lM2) . (134) 
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In order to compare with the notation of DIS, we define the in- 
elastic contributions to the first moments of the spin structure 
functions, 



rr 2 '(Q 2 ) 
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Q 2 )dx=^I h2 (Q 2 ), (135) 



where the integration runs over the Bjorken variable from 
x = (or v — > 00) to x — x (or v = v ). Because the above 
integrals include the excited spectrum only, one has to add the 
elastic contribution in order to obtain the first moment of the 
structure functions. At small values of Q 2 the full structure 
functions are dominated by the elastic contributions, which 
can be constructed once the nucleon form factors are known. 
With increasing resolution Q 2 the coherent response of the 
many-body system "nucleon" decreases, whereas incoherent 
scattering processes on individual constituents become more 
and more important. In particular, the elastic contributions to 
the nucleon structure functions vanish like Q~ w . 

There exist two venerable sum rules for the moments of the 
structure functions. The first one was originally derived from 
current algebra and, therefore, is also a prediction of QCD. It 
deals with the isovector combi nation of the first spin structure 
function in the limit Q 2 -> 00 (Bjo rkenl[l966lll970h . 
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with cja the axial-vector coupling constant and a s the run- 
ning coupling constant of the strong interaction. The other 
sum rule is a prediction for the sec ond spin structure function 
dBurkhardt and Cottin gham, 1970h. 



g»(x, Q 2 )dx = 0. 



(137) 



The Burkhardt-Cottingham (BC) sum rule relies on a "su- 
perconvergence relation" for the associated VVCS amplitude 
such that the dispersion integral exists not only for the odd 
amplitude Si{v) but also for the even amplitude vSi(v). As a 
result the sum of the elastic and inelastic contributions should 
vanish. This allows us to cast Eq. ( 1135b in the form 



1 



I?(Q 2 ) = ~F»(Q 2 ) (FS(Q 2 ) 



F$(Q 2 )) . (138) 



Provided that the assumed convergence criterion is indeed 
given, the BC sum rule connects the spin structure of the ex- 
citatio n spectrum with g round state properties for each value 
of Q 2 . iTsai et al\ d 19751) have proven that the BC sum rule is 
fulfilled for QED to lowest order in the fine structure constant 
a em . Along the same lines, this sum rule is also fulfilled i n 
perturbative QCD to first order in a s dAltarelli et ail Il994l) . 
Furthermore, it is consistent with the Wandzura-Wilczek rela- 
tion, as can be proven by integrating Eq. ( 11341 ) over all values 
of the Bjorken parameter x. However, the BC sum rule is pre- 
dicted to be valid at any Q 2 , whereas the Wandzura-Wilczek 
relation neglects dynamical (twist-3) quark-gluon correlations 
and higher order terms. The different integrals discussed 
above are connected by Eqs. (1122b and ( 1123b . In particular, 
their values at the real photon point can be expressed by the 
charge and the amm of the nucleon, 



/f(0) = I?(0) = \ KN (e N + KN ) : 



(139) 
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3. The helicity structure of the cross sections 

The helicity-dependent cross sections utt{v,Q 2 ) and 
&lt(v, Q 2 ) determine the GDH-like integrals and the 
polarizabilities defined above. For momentum transfers 
Q 2 i5 0-5 GeV 2 , the bulk contribution to these cross sec- 
tions stems from one-pion production, which is reasonably 
well known over the resonance region. The threshold 
region is dominated by s-wave production (Eq + , Sq+) 
accompanied by much smaller contributions of the p waves 
(Mi±, Ei±, Si±). Low-energy theorems, the predictions of 
ChPT, and several new precision experiments have provided 
a solid basis for the multipole decomposition in that region. 
The data basis is also quite reliable in the first resonance 
region. Although the leading Mi + multipole drops with Q 2 
somewhat faster than the dipole form factor, it dominates 
that region up to large momentum transfers. In the higher 
resonance regions, the multipole decomposition is known 
only semi-quantitatively. In particular, there is as yet little 
reliable information on glt, except that it is generally small. 
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FIG. 41 The helicity difference Act = <7 3 / 2 — o"i/2 for the pro- 
ton as function of th e cm energy W for different values of Q 2 
torechsel et all 1 19991) . The figure shows the total helicity differ- 
ence (solid line) as well as the contributions from one-pion (dashed 
line), more- pion (dashed-dotted line), a nd r\ (dotted line) production. 
Figure from Drechsel and Tiator (2004]). 



This fact is not really consoling in the context of the sum 
rules, because the integral Tlt of Eq. ( 1132b does not converge 
well. However, great improvements in the data basis are 
expected from the wealth of ongoing and planned polarization 
experiments. 

The helicity difference Act of the proton is displayed in 
Fig. EH as function of the cm energy W for several values of 
the virtuality. The figure shows negative values near thresh- 
old, because there the pions are produced in s waves lead- 
ing to <Ji/2 dominance. With increasing values of Q 2 , the 
s-wave production drops rapidly. The A(1232) yields large 
positive values because of the strong Ml transition, which 
aligns the quark spins (paramagnetism). For very small val- 
ues of Q 2 , also the second and third resonance regions con- 
tribute with a positive sign. However, this sign has changed 
already at Q 2 — 0.5 GeV 2 . Let us study this effect in more 
detail for the N*(1520) with multipoles E 2 - and M 2 _. Ac- 
cording to Eq. dl 181 ) this resonance yields a term Act ~ 
\E 2 - 1 2 + 6Re (E|M 2 _) - 3|M 2 -| 2 - This value is posi- 
tive at the real photon point where the electric dipole radia- 
tion (El) dominates over the magnetic quadrupole radiation 
(M2). Yet as the magnetic term increases with Q 2 faster than 
the electric one, the helicity difference becomes negative for 
Q 2 w 0.3 GeV 2 . The latter finding is in agreement with per- 
turbative QCD, which predicts the dominance of helicity i 
states at sufficiently large momentum transfer. Figure HTl also 
shows an overall decrease of the resonance structures with in- 
creasing values of Q 2 , because the coherent resonance effects 
are of long range and therefore strongly damped by form fac- 
tors. Finally, for momentum transfer beyond 4 GeV 2 , the res- 
onance structures become small fluctuations on top of a broad 
background, the low-energy tail of DIS. 
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FIG. 42 The Q 2 dependence of the integral If defined by Eq. J135t , 
The open circles show the resonance contribution (W < 2 GeV), 
the s olid symbols also include the DIS contribution. Th e data are 
from l Fatemi et al. ( 20031), CLA S Collaboration (circles ), lAbe et al\ 
Jl998h . SLAC (diamonds), and lAirapetian et al\ ( f2003h . HERMES 
(triangles). Full line: MAID including all channels up to W = 
2 GeV, dashe d line: one-pio n channel only, dotted: 0(p 4 ) prediction 
of HBChPT dJi et aZ.Ll200Ch. dashed-dotted: relativistic baryon ChPT 
teernard et a l., 2002b, 2003), dash-dot-dotted: interpolating formula 
of lAnselmino et a/!Ul989l). asterisk : sum rule value at Q 2 — 0. Fig- 
ure from Drechsel and Tiatqi] i2004l) . 



4. Recent data for GDH-like integrals 

The Bjorken sum rule of Eq. ( 11361 > has been confirmed by 
a series of experimen ts. A fit to all the available DIS data 
dAnthonv et allfloodt) yields the following asymptotic values: 
T{ = 0.118 ±0.004 ±0.007, T" = -0.058 ± 0.005 ± 0.008, 
and hence T{ - J?? = 0.176 ± 0.003 ± 0.007, in good 
agreement with the sum rule prediction of 0. 182 ± 0.005. The 
small value of rf , on the other hand, led to the "spin crisis" 
of the 1980's and taught us that less than half of the nucleon's 
spin is carried by the quarks. With regard to the second spin 
structure function, the BC sum rule predicts that ^(Q 2 ) 
vanishes identically for all Q 2 , and therefore the inelastic and 
elastic contributions should have the same absolute value, 
namely O (Q -10 ) in the scaling limit. 

Figure |42] shows the Q 2 dependence of If. The rapid 
increase from large negative values near the real photon point 
to positive value s in the DIS region is particularly striking. 
The JLab data of F atemi et al.\ (120031) clearly confirm the sign 
change of If at Q 2 0.3 GeV 2 . These data are in good 
agreement with the MAID estimate which covers the same 
energy region. However, with increasing momentum transfer 
the DIS contributions at the higher energies become more 



FIG. 43 The Q 2 de pendence of the integr al I" defined by Eq. (US) . 
The data are from Umarian et al\ ( l2002h . JLab E94-010 Collabo- 
ration (circles). For furt her notation see Fig. [42] Figure from 
iDrechsel and TiatoJ J2004h . 



and more important. Altogether we see a rather dramatic 
transition from resonance-dominated coherent processes at 
low Q 2 to incoherent partonic contributions at large Q 2 . This 
physics is driven by (I) the strong damping of the long-range 
coherent effects by form factors, and (II) the change from | 
to \ helicity dominance in the second and third resonance 
regions. The neutron integral as displayed in Fig. |431 shows 
a similar rapid increase with Q 2 , except that If approaches 
zero right away. The MAID prediction for the one-pion 
channel is in reasona ble agreement with the resonance data 
(lAmarian et all 120021) except for the region of very small 
momentum transfer. This disagreement may have its origin 
in uncorrected binding effects of the "neutron target" in the 
MAID multipoles and/or the data analysis. In Fig. [42] we also 
dis play the pred ictions from heavy baryon ChPT to 0(p 4 ) 
of \hetal.\ (120001) and relativi stic baryon ChPT to 0(p 4 ) 
of Bernard et al\ d2002bl 1200 3l) . The figure shows that the 
chiral expansion can be only applied in a very limited range 
°f Q 2 ^ 0.05 GeV 2 . The value at the real photon point is, 
of course, obtained by inserting the amm as a low-energy 
constant. However, the slope and the curvature of Ii(Q 2 ) 
come about by a complicated interplay of s-wave pion and 
resonance production. Whereas the former process is well 
described by the pion loops of ChPT, the transition form 
factors and widths of the resonances require a dynamical 
treatment of both the A(1232) and the higher resonance 
region. However, ChPT is able to describ e the di f ferenc e 
If — If over a much larger Q 2 range (Burkertj, BOOll) . 
because the contributions of the A(1232) and other isospin | 
resonances drop out in the isovector combination. This opens 



49 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



,1111,1 


i ■ . . . i . . . . 


: \ 




: \ 














_ 

















~.—0~r~ 



0.00 0.25 0.50 0.75 
Q 2 (GeV 2 ) 



1.00 




0.00 0.25 0.50 0.75 
Q 2 (GeV 2 ) 



1.00 



FIG. 44 The Q 2 dependence of the integral IS defined by Eq. fi35\ . 
The ne utron data we r e obta ined by the JLab E94-010 Collabo- 
ration ( Amaria n et all l2004bh . Dash-dot-dotted line: Burkhardt- 
Cottingham sum rule of Eg. dl35b. For further notation see Fig. [42] 
Figure from lDrechsel and Tiatoa ( I2004f) . 



FIG. 45 The Q 2 dependence of the neutron integral /£ T defined by 
Eq. J132t . O pen circles: the resonance co ntribution (W < 2 GeV) 
measured by I Amarian et al\ (2002, 2004b) at JLab, full circles: res- 
onance contribution plus est imate for PIS region. See Fig. I42I for 
further notation. Figure from Drechsel and Tiator ( 20o3)- 



the possibility to bridge the gap between the low and high Q 2 
regimes, at least for this particular observable. 

Figure l44l com pares the neutron d ata of the JLab E94-010 
Collaboration dAmarian et all I2002I) and the MAID predic- 
tion for 7™(Q 2 ) with the BC sum rule. The MAID result 
overshoots the BC sum rule at small Q 2 but agrees with the 
data for Q 2 > 0.1 GeV 2 . At the much larger momentum 
transfer of Q 2 = 5 GeV 2 , the SLAC E155 Collaboration has 
recently evaluated the BC integral in the measured region of 
0.02 < x < 0.8. The results indicate a small deviation from 
the sum rule, but this could well be compensated by contribu- 
tions from the unmeasured region. The integrals I^ T and 7£ T 
have been derived from the 3 He data of t he JLab E94-0 10 Col- 
laboration (Am arian et al. , 2002, 2004b) and corrected for nu- 
clear effects according to lCiofi degli Atti and Scopettald 19971) . 
The data for T^ T show qualitatively the same behavior as the 
discussed integral I". The integral 7£ T is displayed in Fig. [45] 
This observable deserves particular attention, because it sam- 
ples the information from the longitudinal-transverse cross 
section. As indicated by Eq. ( 1132| >, the convergence of 7£ T 
requires that ct^t drop faster than 1/v at large v. Because the 
longitudinal-transverse interference involves a helicity flip, 
this is likely to happen at sufficiently large v. However, there 
is little experimental information on cflt over the whole en- 
ergy region, and therefore the phenomenological description 
is on shaky ground. The zero of I£ T at Q 2 — is particularly 
interesting, because this requires a complete cancelation of 
resonance and P IS contributions. The agreement between the 
new JLab data dAmarian et all l2004bl) and MAID in the res- 
onance region (W < 2 GeV) is quite satisfactory, except for 
the real photon point where both the experimental and the the- 
oretical error bars increase. Furthermore, the contribution of 
DIS is known to be large and negative over the full Q 2 region, 
which brings the integral much closer to zero at Q 2 = 0. Con- 



cerning the ChPT calculations dBernard et al. 1 12003b IJi et all 
2000), the zero value at the photon point is of course taken for 
granted, but the steep slope is a prediction. 

5. Generalized polarizabilities 

In section [VI. A| we have discussed the delicate cancelation 
between negative and positive contributions to the GDH-like 
integrals, and in particular the rapid change of the integrals as 
functions of momentum transfer. For the generalized polariz- 
abilities, the integrands are weighted by an additional factor 
v~ 2 or x 2 , which enhances the importance of the threshold 
region relative to resonance excitations, and suppresses the 
contributions of the DIS continuum above W = 2 GeV. Fig- 
ures @6] and @7] display the transverse-transverse (jtt) and 
longitudinal-transverse (Slt) polarizabilities as function of 
Q 2 . A s in the case of I^ T , the data for 7^ T dAmarian et all 
2004a) show considerably more strength at small Q 2 than pre- 
dicted for the one-pion contribution. However, the agreement 
f° r Q 2 si 0.4 GeV 2 is again quite satisfactory. In view of 
the additional weight factor towards the low-energy region, 
this behavior is another indication that the "neutron problem" 
near the real photon point should be related to low-energy and 
long-range phenomena. A comparison with the predictions 
of ChPT shows that 7^ T is a particularly sensitive observable 
because of the cancelation between s-wave pion production 
and A resonance excitation, as is apparent from Eq. ( II 181 ). 
In fact, the additional weight factor v~ 2 increases this cance- 
lation considerably relative to the integral Itt- It is there- 
fore no big surprise that ChPT cannot describe jtt with- 
out including the A resonance. And indeed, the 0(p 3 ) and 
0(p ) approximations of HBChPT change from positive to 
negative values, and also the newly developed Lorentz invari- 
ant version of ChPT misses the real photon point. This be- 
havior changes for the longitudinal-transverse polarizability 
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FIG. 46 The Q 2 dependence of the generalized neutron polarizabil- 
ity 7tt defined by Eq. ( 1131b . The open circ l es are t he data of the 
JLab E 94-010 Collaboration dAmarian e?oZll2004al). See Fig. |42] 
for further notation. Figure from Drec hsel and Tiatog ( l2004l) . 
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FIG. 47 The Q 2 dependence of the generalized neutron polarizabil- 
ity S2t defined by Eq. d 1 33fc . Note that <5£ T has been multiplied by 
a factor of 10 4 Q 6 /(16a e mA/p) in order to compensate for its rapid 
decrease with increasing Q 2 . The open circles are t he data of the 
JLab E 94-010 Collaboration jAmarian et all l2004al). See Fig. |42] 
for further notation. Figure from lDrechsel and Tiatorl | |2004|) . 



shown in Fig. |47j because Slt is dominated by the s-wave 
term Sq + Eq + in the multipole expansion. The contribution of 
the A resonance is proportional to S* + E±±, which is much 
suppressed because both amplitudes are small. As a conse- 
quence S lt decreases rapidly as function of Q 2 without show- 
ing any pronounced resonance structures. In other words, the 
longitudinal-transverse polarization takes place in the outer 
regions of the nucleon and is mostly due to the pion cloud. 
This notion is well supported by the fact that the ChPT pre- 
diction for 5lt is much better than for jtt- 



VII. CONCLUSION 

In this review we have concentrated on the bulk prop- 
erties of hadrons seen at low momentum transfer, such as 
shape, polarizability, and low-energy excitation spectrum. 
We have presented many new precision data implying the 
compositeness of the hadrons based on the interaction among 
the constituents, quarks and gluons, as prescribed by the 
QCD Lagrangian. Experiments at high momentum transfer 
have established this theory as the fundamental theory of the 
strong interactions, with previously unexpected features like 
asymptotic freedom and chiral symmetry. The great success 
of QCD at high Q 2 has become possible by the decrease 
of the strong coupling constant with increasing Q 2 , which 
allows for a perturbative treatment of QCD in this region. 
This is definitely not possible in the Q 2 region below 1 GeV 2 , 
which defines the realm of non-perturbative QCD. Therefore, 
QCD has not yet passed its final tests in the low-energy 
domain, in which new phenomena show up, as for instance 
the confinement of the many-body system "hadron" and 
the spontaneous breaking of the chiral symmetry. In fact, 
QCD encounters very fundamental conceptional problems 
at low Q 2 : One has to deal with the relativistic many -body 
aspect, the theoretically not yet understood confinement, 
and the fundamental strong interaction at the same time. 
The simplifications of the parton model typical for the high 
Q 2 physics do not work in the non-perturbative region, and 
neither is it possible to describe the bound many-body system 
of the light u and d quarks by basically non-relativistic 
physics as is done for systems of heavy quarks. 

This review has two punch lines. On the experimental side 
one has to insist on measurements with the utmost precision. 
As it became evident only a new generation of cw accelerators 
together with modern detectors and data acquisition allowed 
for progress after almost two decades of stagnation. How- 
ever, since the most interesting observables like spin observ- 
ables are still difficult to get at, the experiments have to be 
well chosen. They have to be "significant" for contributing to 
the understanding of hadrons in the framework of QCD. For 
this purpose a profound theoretical guidance is needed. As 
promising examples of such significant experiments we men- 
tion: 

• A more complete and even more precise study of meson 
threshold production, in particular through a study of 
spin observables (see section[VA} and the production of 
strange quarks, for example by kaon electroproduction. 

• Real and virtual Compton scattering in order to dis- 
entangle all the spin polarizabilties of the nucleon by 
double-polarization experiments (see section HVl). 

• A fresh approach to investigating the excitation spec- 
trum of the nucleon, with more attention given to the 
resonance-"background" separation (see section IV.BI ). 
Once more, such dedicated experiments require an in- 
tense study of the spin observables. Such investigations 
are also the basis for further progress in understanding 
the sum rules discussed in section [VTl 
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• The nucleon form factors remain a challenge at both 
low and high Q 2 (see section [TTTb . The improvement 
of the data base at low Q 2 would however take a very 
strong effort with the existing facilities. It is in place 
here to regret the closure of the MIT/Bates storage ring 
which offered the best means to perform such experi- 
ments at low Q 2 . 

On the theoretical side the challenges are not smaller. In 
this review we have met several effective models such as con- 
stituent quark models, dynamical baryon models, and effec- 
tive meson field models. All are effective in the sense that they 
describe the data more or less well and provide an insight into 
the physics. They also provide the necessary guidance for the 
experiments. However, these models are only "inspired" by 
QCD but not fundamentally based on this theory. As of today 
only two schemes exist that are both directly based on QCD 
and able to describe the low-energy region: chiral perturbation 
theory (ChPT) and lattice gauge theory (LGT). Ironically the 
first one works best for very small quark masses, whereas the 
second one requires large quark masses, at least with the com- 
puter power of today. The extrapolation between these two 
approximations is therefore still on shaky ground. However, 
the ongoing efforts promise decisive turns in the near future: 

• The possibilities for calculating loops beyond the lead- 
ing order in ChPT is improving continously, and the 
necessary low-energy constants get better and better 
constrained by experiment and, maybe, in the near fu- 
ture also by LGT. 

• Lattice gauge theory is progressing through more ef- 
fective algorithms and the growing performance of the 
computers. It appears that the used quark masses will 
soon approach sufficiently low values to make contact 
with ChPT by a chiral extrapolation of the numerical 
results from LGT. This would indeed mark an essential 
breakthrough in low-energy hadron physics. 

As became evident in this review, the agreement between 
experiment and theory is not yet satisfactory in quite a few 
cases, which leaves major challenges in hadron physics for 
both sides. 

We conclude that the study of hadrons at low energy and 
momentum transfer is a unique possibility to unravel the struc- 
ture of a relativistic many-body system in the realm of QCD, 
one of our most fundamentally bas ed quantum fie ld theories, 
as is highlighted by the quote from lWilczekl d2000l) : "QCD is 
our most perfect physical theory." 
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